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Douc as, Long Beach Division, engineers 
used CEC Dynamic Recording Equipment 
to test the new “Single-Point Fueling” 
system for the C-124 Globemaster. The 
12 fuel tanks, 6 in each wing, hold ap- 
proximately 11,000 gallons—equal to a 
full tank car. Fueling normally requires 4 
to 5 hours. The new system cuts this time 
to 18 minutes by pumping through a single 
entry at 600 gpm, filling all 12 tanks at 
once. 

By using three CEC 5-114 Recording 
Oscillographs and five sets of System “DD” 
amplification (shown at left), Douglas was 
able to simultaneously record all tank in- 
let pressures, tank internal pressures, 
manifold pressures, tank shut-off time, 
pressure drop through the main inlet noz- 
zle, and the fuel level in each tank. In all, 
the findings of 29 pressure pickups and 12 
capacitance fuel gages were recorded. 


Dynamic Recording Systems 


such as the one shown here are designed and manufactured 
by Consolidated. Variations in the arrangement of the 
equipment are infinite. Applications are widely varied 
throughout industry and the sciences. A typical recording 
system includes pickups, amplifiers or bridge balances, and 
a recording oscillograph. Write for Bulletin CEC 1500B. 
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PROPULSION 
FOR A MISSILE 


The art of propelling a missile has progressed a long 
way since the era of the rock-throwing Roman catapult. But 7 


the design of a modern missile, like that of the old 


stone catapult, is best done by those with missile experience. — 


Engineers at Fairchild’s Guided Missiles Division 
are among the most experienced in their field. 
Beginning with one of the Armed Services’ very 
first missiles projects, Guided Missiles Division 
engineers have played an important role in the 
design and development of complete modern 
missile weapons systems. Fairchild missile 
projects have included both rocket and 


turbo-jet powered missiles. 


Fairchild’s broad experience encompasses all 


phases of missiles weapons systems, including 
_ propulsion, airframe, guidance and such intricate 


associated equipment as ground and 


shipboard radar. 


Pe 


ENGINE AND AIRPLANE CORPORATION 


FAIRCHILD 
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—the Douglas F3D Skyknight 


Out of Korea come new reports of the 
Douglas F3D Skyknight in action, down- 
ing Migs for the United States Marine 
Corps during advanced night and foul 
weather operations. 

Designed for the U. S. Navy, the all- 
weather Skyknight flies at near-sonic 


Be a Naval flier—write to 
Nav Cad Washington 25, D.C. 


speeds, operates from aircraft carriers as 
well as small advanced airfields. A side- 
by-side seating arrangement of pilot and 
radar operator results in closer combat 
teamwork—permits Skyknight’s modern 
radar search and fire control equipment 
to be operated with maximum efficiency 


when against marauding enemy planes, 

Performance of F3D Skyknight in 
action is another example of Douglas 
leadership in aviation. Planes that can 
be produced in quantity to fly faster and 
farther with a bigger payload are a basic 
rule of Douglas design. 


ARS JourNAL 


impt 
We 

since 
max! 
weig 
for a 
desis 
weig 
auxi 
sis sl 


one 
2 
‘ie 
Radar eyes see in darkness, storm, or 
to lock this twin-jet fighter on its pre 
= 
A 
fluic 
valu 
and 
ity 
dise’ 
i 
Depend on DOUGLAS First in Aviati Socr 
ia epen on irst in Aviation peri 


\L 


SEPTEMBER- OCTOBER 1953 @ Ki [3 Editor-in-Chief 
VOLUME 23 NUMBER 5 MARTIN SUMMERFIELD 
socierTy 


= 
Coma of and Fluids for 


Rocket Propulsion | 


i J 
> KRAFFT A. EHRICKE! 
oe Guided Missile Development Group, Redstone Arsenal, Huntsville, Ala. 

A series of chemical propellants and atomic working If the project net weight becomes larger than the assigned 
fluids, covering a wide range of density-specific impulse net weight—for a given gross weight—it is technically not 
values, are compared with respect to vehicle performance feasible to use this particular fluid, unless the cut-off velocity 
and design parameters. The changing influence of dens- and/or payload weight are lowered. 
ity and specific impulse with increasing performance is Comparing several propulsion fluids, one can define an opti- 
discussed. Ranges for favorable application of different mum fluid which requires a minimum of fluid weight and net 
fluids are indicated and their respective limits of applica- weight per unit payload weight. Minimization of the net 
bility shown for ascent in a gravity field as well as for space weight should be understood as reduction of the weight of 
flight. those items which are influenced by the fluid selection; that 

-* ——~ a is, mainly the power plant weight, determined by propulsion 
oe oe . system and tank volume. 

=P ) Introduction? Obviously, this fluid is most desirable from a general per- 

: : : formance viewpoint. Its rating for a particular project de- _ 
CHEMICAL propellant or a working fluid for thermonu- pends on additional considerations of great practical impor- 

clear propulsion determines the mechanical performance tance, such as economy, availability, handling, storage, mr 

of a rocket vehicle through its characteristic values of density portation, effects on power plant development, and others. 
and specific impulse. The interesting parameters describing These viewpoints have been excluded intentionally in the 
the mechanical performance are cut-off velocity, attained aforesaid definition, in order to maintain its general validity. 
upon termination of powered flight along a given trajectory, This paper deals with the assigned net weight, or its dimen- 
and payload weight carried. : oll ob ase sionless parameter, in relation to propulsion fluid and me- 

The gross weight of a rocket vehicle can be subdivided into chanical performance with the purpose of determining the ap- 
payload weight, propulsion fluid weight, and the residual proximate range of favorable applicability of individual fluids. ; 
weight which shall be called net weight, to be discussed below. The conclusions arrived at should be understood in the light of 

If, for a certain gross weight, the payload weight and cut-off the above remark that, for practical reasons, a different fluid — 
velocity along a specific trajectory are given, then, by select- may be found preferable. It is believed, however, that the aos 
ing a propulsion fluid, the fluid weight and consequently also results may be found useful in that they show which pro- 
the net weight are fixed. Comparison between different pulsion fluids deserve greatest attention from the viewpoint of 
fluids will yield different fluid weights, depending on specific research and development in various performance ranges up 
impulse and density, and therefore also different net weights. to interplanetary flight. ans 
We shall designate this net weight as assigned net weight, : Se 
since it results from the given gross weight and represents the z . bus 
maximum weight available for the items comprising the net General Analysis 
weight. Whether or not the assigned net weight is sufficient For powered flight in a gravity field, vertical ascent only 
for actual construction and flight depends on a more detailed will be considered for the sake of simplicity. 7 
design and project analysis with given fluid and purpose of From the fundamental equation of rocket motion, the cut- ‘ 
project, since the fluid density influences the construction off velocity is then given by 
weight and the project purpose has a bearing on the weight of . 
auxiliary devices. The net weight resulting from this analy- cols [1] 
sis shall be designated as project net weight. mi 
where denotes the cut-off velocity, g the gravitational ac- 

= celeration at the surface, J;, the specific impulse, m and m, 

Received September 23, 1952. a _ the over-all mass of the vehicle at take-off and cut-off, re- 

Chief, Gasdynamics Section. Member ARS. spectively, and ¢, designates the burning time, given by the © 


? At the time of completion of this study, a paper b a 
Tsien was published in the JouRNAL OF THE AMERICAN atin relation 
Socrery, Vol. 22, No. 4, pp. 200-203, which compares the Top 


performance of different power plants by a method similar to the [2] 
one derived in this paper. : c 
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 pellant weight over gross weight of the vehicle, and C,, is the 
absolute acceleration at take-off in terms of g, defined by 


| 


Fis thrust. 
The equation for the cut-off velocity can then be written 
in the form 


The mass ratio is connected with the propellant density by 
the relation 


‘iiheiae Yp denotes the mean specific weight of the propellant, 
and ¢ is the propellant capacity factor of the vehicle, defined 
___ as the ratio of available propellant container volume per unit 
empty weight 

W designates the propellant weight. 
> The loading factor can then be expressed in terms of pro- 

pellant specific weight and capacity factor 


Wo 


The cut-off velocity, as given in Equation [4], depends on 
y specific impulse, mass ratio, loading factor, and initial accelera- 


tion. Since the last parameter is of minor interest for this 
discussion, a given value shall be assumed, say Cy, = 2. Then, 
using Equations [5] through [7], the cut-off velocity becomes 


-= makes v, a function of the interesting parameters Isp, 


= glep [im (1 + ype) 20 


Yp, and ¢ only. 

__In order to be applicable for performance calculations, the 
propellant capacity factor has been defined in Equation [6] 
in terms of the empty weight of the vehicle. This includes 
not only the construction weight, but also the payload weight. 
Therefore, the capacity factor is not strictly a function of the 
design quality of the vehicle, but depends also on the payload 
weight carried. The capacity factor can be connected with a 


weight parameter by the following consideration. Let the 
net weight of the vehicle be defined by 


where W;, is the payload weight. The term Wy comprises the 
residual weight of the vehicle; that is, construction weight 
W., weight of the auxiliary fluids W4r for energizing the prime 
mover of the feed system, and additional weight for equipment 
and installations which cannot be considered as payload and 
therefore are added to the construction weight to give the net 
weight. (The term net weight has been preferred to the term 
empty weight because the latter includes the payload 
weight.) Although the auxiliary fluid contributes to the over- 
all thrust, its effect is comparatively small and shall be ne- 
glected in this case for the sake of simplicity. 

Using net weight and propellant weight, a net weight factor 
can be defined as 


Using the net weight factor, a simple relation can be de- 
rived for the payload factor, defined as ratio of payload weight 
over gross weight of the vehicle or of the stage under consider- 
ation 


where A denotes the loading factor, that is, the ratio of pro- 


Substituting the capacity factor in Equation [8] by expres- 
sion [12], and making the cut-off velocity a dimensionless per- 


(11 


Yrh 
1 — fy (1 —fw)(1 + yp) 


From Equation [11] the capacity factor follows as function 
of net weight factor, payload factor, and specific weight of the 


propulsion fluid 


- 


- FF) | 


o= 


formance parameter by dividing it by g/sp, yield 
vy, 1 
, = 
0 


In Equation [8] a value of Cy, = 2 was inserted. 


In the preceding analysis, Equations [8] and [13] define 
the relationship between the three sets of parameters describ- 
ing the propulsion fluid (yp, Isp), the vehicle performance 
(vi, Wi/W,), and the vehicle design (¢, fy). Consequently, 
for a given performance a set of design parameters can be as- 
signed to each propulsion fluid considered. These assigned 
values have a considerable practical bearing on the relative 
merits of the fluid and even on its applicability under given 
conditions. 

The capacity factor indicates the vehicle size required for 
the particular fluid on account of its density and specific im- 
pulse. The net weight factor, in turn, is a function of capa- 
city factor and fluid density if the payload factor is given. 
This is shown by Equation [11] which takes the influence of 
vehicle performance into account, since it contains the loading 
factor. Inspection of Equation [11] shows that for a given 
payload factor the assigned net weight factor becomes largest 
for the fluid which requires the smallest loading factor. In- 
deed, from the viewpoint of technical feasibility this fluid is 
most favorable, since a large assigned net weight factor facili- 
tates design and construction of the vehicle. 

However, the seeming advantage of this fluid over another 
fluid with smaller assigned net weight factor can be utilized in 
practice only if the increase in assigned net weight factor is 
larger than the increase in project net weight factor. This is 
not necessarily the case. Equation [10], which involves 
vehicle data only, shows that the net weight factor must in- 
crease when only the propellant weight is reduced. If, in ad- 
dition, the absolute net weight becomes larger, for instance, if 
the alternate fluid is lighter (heavier feed system, plumbing, 
insulation, etc.), then the resulting increase of the project net 
weight factor is even greater. Comparison of the new values 
obtained from Equations [10] and [11] indicates in the first 
approximation whether or not the theoretical superiority of 
the alternate fluid can be made effective in practice In the 
same manner, different mixture ratios of a propellant com- 
bination can be compared. 

These considerations apply also when a group of propulsion 
fluids is compared in the light of increasing vehicle perform- 
ance. In this case the assigned net weight factors of all 
fluids will decrease as the performance is increased. In cer- 
tain performance ranges they will not be within the limits of 
technical feasibility. These limits will be different for the 
various fluids. Their exact determination again requires @ 
detailed design analysis. However, knowledge of the cor- 

' related capacity factors is useful for a preliminary estimate. 
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The Capacity Factor 


For this reason, it is of interest to study the change of the 
capacity factor with density and specific impulse for a given 
cut-off velocity. For this purpose, a number of propellants 
and working fluids have been considered with the aim of cover- 
ing a wide variety of values of specific impulse and density. 
The data as used here are taken from a previous publication 
in which the performance of various fluids was calculated for 
shifting equilibrium as function of mixture ratio and pressure 
ratio (1).3 

The performance data of working fluids and additional per- 
tinent data are summarized in Table 1. The calculations 
were based on a temperature before expansion of 3000° K and 
on a chamber pressure of 20 atm. Two exit pressures were 
chosen, representing operation at low and high altitude, re- 
spectively. For mercury, a higher chamber pressure of 40 
atm was taken and the expansion limited to 1 atm, in order to 
avoid condensation. The mercury values should be consid- 
ered aS approximations only, since thermodynamic data at 
higher temperatures were not available. For heating up to 
648° K a specific heat of 6.52 cal/mole was estimated (2). For 
mercury vapor, the specific heat of a monatomic gas was 
taken. This is certainly not correct in the lower vapor tem- 
perature range where real-gas corrections would lead to higher 
values. However, the performance of mercury is low and its 
application for propulsive purposes would hardly be consid- 
ered, so that an approximation may be sufficient. Dissocia- 
tion of diatomic gases was taken into account, except in the 
ease of nitrogen. 

The propellants considered are oxygen-hydrogen, a rep- 
resentative of light combinations, ozone-methane, oxygen- 
ammonia, oxygen-hydrazine, designating medium propellants, 
and hydrazine with the heavy oxidizers, nitric acid (RFNA), 
tetranitromethane, and chlorine trifluoride. Except in the 
combination with ozone, only hydrogen and hydrogen-nitro- 
gen fuels were used, since they yield a higher specific impulse 
than most carbonaceous fuels and, in the case of hydrazine, 
show a favorably high density. Methane is especially useful 
in combination with ozone, since its boiling point lies 50 
centigrades below that of ozone, but still above the freezing 
point of the latter. The main hazard in using ozone seems to 
be its explosive vapor rather than the liquid which appears to 
be stable if kept free from impurities which might act as cata- 
lysts (3). Ar:anging the ozone tank inside the methane 
tank would reduce the vapor pressure of ozone to 12 mm of 
mercury, using the vapor pressure data given in reference (4), 
when the vapor pressure of methane is 1 atm (5). In com- 
bination with-methane the use of ozone might become practi- 
cable and therefore this oxidizer has been included in the com- 
parison. 

Fig. 1 presents the maximum theoretical specific impulse for 
different expansion ratios of propellants and working fluids, 
plotted against their mean density or density (working fluids), 
in liquid state. The mean density of the propellants in- 
teases slightly with decreasing exit pressure, since the mix- 
lure ratio for maximum specific impulse shifts toward the lean 
side. 

The mixture ratio corresponding to maximum specific 
impulse is, in most cases, not the optimum mixture ratio for 
maximum flight performance, except when the densities of 
both propellant components are similar—as in the case of 
xygen-hydrazine. For propellants with greatly different 
component densities, the mixture ratio for maximum flight 
performance is considerably closer to stoichiometric—the 
oxidizer being the heavier component—or even in the range of 
widizer excess. It is not the purpose of this paper to deter- 
mine these optimum mixture ratios. However, for reasons of 
comparison, a particularly notable case, that of oxygen-hy- 
drogen, is included. The performance of this propellant is 
chown in Fig. 1 for the mixture ratios yielding maximum spe- 


* Numbers in parentheses refer to References on page 296. 
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CHARACTERISTIC AND PERFORMANCE DATA OF WORKING FLUIDS FOR THERMONUCLEAR PILE PROPULSION 
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cific impulse at the respective pressure ratios and for stoichio- 
metric ratio, r = 7.94, where r designates the weight ratio of 


e = 20/latm 
° R/Re 20/07 
' 20/01 
800) 20/0.01 
x 40/1 
RFNA + 85% Yo NOz 
e Oe | 
700 Nel He 
He NHs Cig | O2 | A Hg 
| 
| 
i | 
400 
L | 
300g" 
3, 
a | 3 
w 
z 
IN x 
—{_I|scace 
ie ] if | | 
MEAN SPECIFIC GRAVITY ] 
08 09 10 12 13 14 «15 14 
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oxidizer to fuel. This mixture ratio has a higher density, but 
it is not the correct optimum value which lies in the region of 
large oxygen excess and, therefore, has little practical interest, 

The exhaust velocities g/;p and the specific weights of all 
fluids are listed in Table 2 (see page 294), together with com- 
ponent densities and mixture ratios for the propellant combina- 
tions. As nitric acid, RFNA with 15 per cent NO» was taken, 

The specific impulse of most working fluids is seen to be in- 
inferior to that of chemical propellants. The difference is so 
large in many cases that the gas temperature of the working 
fluids would have to be raised substantially, in order to obtain 
comparable values. Notable exceptions are, of course, liy- 
drogen and helium, but also ammonia and hydrazine. 

Figs. 2 and 3 show the cut-off velocity obtained as function 
of the capacity factor for the different propellants and working 
fluids. Most of the propellants show little difference with the 
exception of O. - NH; and O» - He, both mixture ratios, which 
require higher capacity factors for equal performance. Since 
the dimensions of a rocket vehicle are mainly dependent on 
its capacity factor, these graphs indicate the relative velicle 
size required for different propulsion fluids as a function of 
flight performance. Among the propellants only oxygen- 
ammonia and oxygen-hydrogen yield a particular increase in 
vehicle dimensions, The capacity factor of the A-4 (\-2) 
rocket, using oxygen and diluted alcohol, is about 2.1 cu m per 
metric ton for a cut-off speed of about 5000 fps. For the 
same velocity, O2 - H, for maximum specific impulse requires 
3.5 cu m/t and (Os: - He)s:, the stoichiometric mixture ratio, 
2.7 cu m/t, which means that the vehicle would be about 70 
per cent or 30 per cent larger, in spite of superior specific im- 
pulse. The heavy propellants which require the smallest 
capacity factor at low velocities are replaced gradually, with 
increasing performance, by combinations of medium density 
but higher specific impulse. Among the working fluids shown 
in Fig. 3, the dominant role of density for low-speed vehicles 


i6 
ea 
LZ 
= 
| 
| ex) 
‘aa 
< 
LZ 
| 
S 
16 32 48 64 
0 2 4 feum/t] 6 


CAPACITY FACTOR, @ 
FIG. 3 INFLUENCE OF CAPACITY FACTOR ON CUT-OFF VELOCITY 
FOR DIFFERENT WORKING FLUIDS (EXPANSION RATIO 20/1 ATM) 


ARS JourNAL 


—= 


[ft/sec] 


id 


IDEAL CUT-OFF VELOCITY, Vv, 


a 


FIG 


a 
= 
» 
ste 
2 
he 
cel 
Q 
V4 ar’ 
‘wal 
acc 
: 
3 
~ 
3 
3 
ae 
— 
E 
~ 
A 
8 
6 
I 
i 


but 
1 of 
est. 
“all 
m- 
cen, 
$0 
‘Ing 
in 
hy- 


rion 
‘Ing 
the 
lich 
nee 
on 
icle 
1 of 
ren- 
in 
per 
the 
ires 
tio, 
70 
im- 
lest 
vith 
sity 
Wn 
cles 


3 


IDEAL CUT-OFF VELOCITY, 1073 


S 
T 


and the growing importance of specific impulse are illustrated 
even more clearly in the case of mercury and argon. Figs. 2 
and 3 are based on Equation [8]. The specific impulse cor- 
responds to expansion from 20 to 0.7 atm for the propellants, 
and 20 to 1 atm for the working fluids. 

A wider range of capacity factors is presented in Figs. 4 and 
5 which are based on gravity-free powered flight, yielding the 
so-cilled ideal velocity at cut-off. This applies to horizontal 
propulsion in vacuo as, for instance, in the upper stage of 
satellite vehicles and propulsion of space ships. Due to the 
special construction of space ships, higher capacity factors 
will become feasible. The specific impulse and density used 
for computing Figs. 4 and 5 are based on 20 to 0.1 atm expan- 
sion ratio. Equation [8] has again been applied, putting the 
second term in the brackets equal to zero, since Cy, formally 
becomes infinite. The rapidly growing importance of ozone- 
methane is notable, while the denser propellants now require a 
larger capacity factor for the same performance. The largest 
factors still are needed for the oxygen-hydrogen mixtures. 
Among the working fluids, hydrazine and ammonia are out- 
standing, followed by water. The curves for hydrogen and 
helium show the steepest gradient, indicating the rapidly 
growing importance of high specific impulse. 

The changes in Figs. 4 and 5, as compared to vertical as- 
cent, are due to the fact that gravitational losses in Figs. 2 and 
3 affect the fluids to a different degree, because constant ini- 
tial acceleration rather than constant burning time has been 
assumed. This yields reduced burning time for the light 
fluids at small capacity factors and for the heavier fluids at 
large capacity factors. However, constant burning time 
would have involved a considerable variation of the initial 
acceleration for the different fluids. This hardly would be 
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practicable in actual operation. Due to the differences in 
density, the same capacity factor corresponds to a different 
loading factor (Equation [7]). Therefore, the amount to be 
subtracted as gravitational loss at vertical ascent (the term 
A/Cy, in Equation [4]) is different. With increasing loading 
factor, the natural logarithm of the mass ratio grows more rap- 
idly than the term A/C,, in Equation [4]. Consequently, 
vehicles with a large loading factor, using dense propellants, 
are comparatively less affected. For gravity-free powered 
flight, in turn, light fluids become comparatively more efficient 
as far as small capacity factors are concerned. 

Comparison between Figs. 2 and 3 and between Figs. 4 and 
5 shows that most working fluids require larger capacity fac- 
tors than the propellants for equal cut-off velocity. Excep- 
tions are hydrazine and ammonia. 

Fig. 6 shows the connection between capacity factor and 
mass ratio for the fluids under consideration. The little cir- 
cles in this figure terminate the lines corresponding to mixture 
ratios for pressure ratio 20 to 0.7 atm, used in Fig. 2. The 
loading factors which correspond to the mass ratios indicated 
in Fig. 6 can be found immediately, since 


(mo/m,) — 1 
mo/m, 


A comparison between Fig. 6 and the preceding performance 
diagrams reveals the significance of different density and spe- 
cific impulse in terms of capacity factor and mass ratio. Take, 
for instance, CIF;-N2H, and O.-He, both mixture ratios. 
Fig. 2 shows that for a cut-off velocity of 5000 fps the re- 
quired capacity factors are, in cubic meters per metric ton 
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The corresponding massratiosare 


ee 2.4 2.14 1.93 


In other words, oxygen-hydrogen requires slightly lower mass 


ratios, but much larger capacity factors. 

Comparing this result with a high performance case, as- 
suming, for instance, an ideal velocity of 26,000 fps, the re- 
quired capacity factors and corresponding mass ratios are, ac- 
cording to Figs. 4 and 6 
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FIG. 6 MASS RATIO AS FUNCTION OF CAPACITY FACTOR 


This ideal velocity corresponds approximately to the value re- 
quired for cireumnavigation of the moon, starting from a 
satellite at 2500 miles altitude and establishing a circular 
orbit at 30 miles above the surface of the moon. 

These two examples illustrate the increasing superiority of 
lighter propellants of high specific impulse with increasing 
performance. The difference in mass ratio, hence in loading 
factor widens. Simultaneously the difference between the 
respective capacity factors is reduced. Still, the capacity 
factors indicate that the hydrogen vehicle would become con- 
siderably more bulky in either case, particularly in the first 
example. This would raise practical difficulties in terms of de- 
sign and transportation. It also means that, in the first ex- 
ample, the oxygen-hydrogen vehicle gross weight will hardly 
be smaller than that of the ClF;-H.H, vehicle due to the very 
large difference in capacity factors. With increasing cut-off 


velocity the differences in gross weight will be reduced gradu- 


ally, and eventually the light-propellant vehicle can be ex- 
pected to have the lower gross weight for a given payload. 


However, even then this vehicle will have the higher con- 


struction weight. 


The Net Weight Factor 


bes 5 As pointed out before, the net weight comprises the weight 
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of construction, of the auxiliary fluids and additional weight, 
Waa, depending on the project purpose. Therefore, the net 
weight, as defined in Equation [10] is, for a given propulsion 
fluid, not exclusively a function of the construction weight. 

For a given propellant, the net weight factor depends on 
size and purpose of the vehicle. The factor will decrease with 
increasing vehicle size, because of the more favorable ratio 
between volume and surface area. For equal size, the pur- 
pose of the vehicle has a decisive bearing on the net weight 
factor. A ballistic rocket, for instance, can be expected to 
have a smaller factor than a glider rocket. The A-4 (V-2), a 
ballistic rocket, has a net weight factor of about 0.28. lor 
its winged version, the A-4b which was fired at Peenemiinde 
in 1944-1945, the net weight factor was about 0.31. Both 
vehicles were operated with liquid oxygen and diluted alcohol. 
Smaller net weight factors can be expected for future booster 
rockets in multistage vehicles. The large planned 
booster rocket A-10 with a thrust of about 200,000 kg and the 
A-4b or a similar vehicle as payload stage (13,400 kg) had a net 
weight factor of not more than 0.24 or 0.25. This version 
was supposed to use a combination of nitric acid (RFNA) «and 
a heavy crude oil, the mean specific gravity being 1.33 as com- 
pared to about 1.0 for oxygen-alcohol. 

If the vehicle size is given, the net weight factor decreases 
rapidly with increasing propellant density. As an example, 
the Peenemiinde rocket project A-8 may be mentioned, which 
was laid out for nitric acid and crude oil. It had a gross 
weight of about 22,700 kg at 1000 kg payload. The pro- 
pellant tank volume was 13 cu m (459 cu ft) and its net weight 
was about 4420 kg. This yields a net weight factor of roughly 
0.2, mainly due to the large propellant weight of 17,300 kg. 
If the tanks were filled with oxygen-alcohol, keeping the net 
weight unchanged, the net weight factor would be raised to 
0.25. If oxygen-hydrogen, at a mean specific gravity of 0.27, 
would be the propellant, the net weight factor would become 
0.55, keeping the net weight constant. Using oxygen-hy- 
drogen in stoichiometric mixture (mean specific gravity 0.424), 
the net weight factor would become 0.43, again at constant 
net weight. These are the effects of changing propellants 
without consideration of correlated changes in absolute net 
weight. A net weight increase of the order of 30 and 20 per 
cent for the lighter and heavier oxygen-hydrogen combina- 
tion, respectively, would have to be expected, mainly due to 
increased propulsion system weight and additional insulation. 
This would raise the net weight to 5500 and 5100 kg, respec- 
tively, yielding net weight factors of 0.61 for the lighter and 
0.48 for the heavier mixture. These values are considerably 
higher than the one required for the heavier propellant, al- 
though the influence of mixture ratio in the use of oxygen-hy- 
drogen is quite strong. 

These examples indicate the magnitude of net weight fac- 
tors to be expected for present-day vehicles. 

For very large vehicles of the future, a substantial net 
weight factor reduction can be expected. As an example, the 
weight study for a three-stage vehicle may be mentioned 
where a winged third stage is supposed to carry 11,000-lb 
payload into a circular orbit at 138 miles altitude. The’ main 
data are summarized in reference (6). In this reference the 
given propellant weight includes the propellant used for gas 
turbine feeding (2 per cent). If this propellant is added to 
the construction weight to obtain the net weight, the net 
weight factors become about 0.13 for the first stage and 0.11 
for the second stage. The first stage was laid out for nitric 
acid (RFNA) and hydrazine, the second and third stage for 
oxygen and hydrazine. The absolute initial accelerations as- 
sumed were 1.5 g, 1.1 g, and 1.0 g, respectively, for the first, 
second, and third stage. The net weight factor for the third 
stage becomes about 0.4, although this vehicle is somewhat 
larger than the A-4b, having a factor of 0.31. The main rea- 
son for this increase in net weight is that the additional weight, 
due to pressurized cabin, cooling system, landing gear, etc., 
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PAYLOAD FACTOR AS A FUNCTION OF LOADING FACTOR FOR 
DIFFERENT NET WEIGHT FACTORS 


is considerably larger than in an unmanned rocket glider of 
the guided missile type. 

Equation [11], which is of general validity, shows that the 
payload factor decreases with increasing net weight factor. 
The exact amount of change depends on the loading factor. 
In Fig. 7 the payload factor is plotted as a function of the 
loading factor for different net weight factors. For a given 
cut-off velocity the loading factor may be decreased if the 
specific impulse can be increased, allowing for a higher pay- 
load factor and/or net weight factor. The plot shows that an 
increase in net weight factor by the amount mentioned before 
in connection with oxygen-hydrogen calls for a substantial 
decrease of the loading factor. That this is difficult to obtain 
with smaller vehicles was illustrated in the previous section. 
If the gross weight of the light fluid vehicle is lower than that 
of the heavy fluid vehicle, as is the case at very great vehicle 
size, then the payload factor must also be increased, if a given 
absolute payload weight is considered. This limits the 


permissible increase of the net weight factor. tenia 
Application 


From Fig. 2 through Fig. 5 the capacity factor required for 
using a given fluid can be found as function of the cut-off 
velocity for vertical ascent. Similar curves can be drawn for 
ascent along any desired trajectory. A generalized representa- 
tion, even for a given trajectory, is not possible due to the 
particular combination of density and specific impulse for 
each fluid. 

The assigned net weight factor can be presented in general 
form as function of the dimensionless velocity parameter de- 
fined in Equation [13]. This is done in Fig. 8 for the two 
cases of vertical ascent and gravity-free powered flight, taking 
the payload as parameter. With the aid of Table 2, the as- 
signed net weight factor for a given cut-off velocity and pay- 
load factor ean then be found immediately for different pro- 
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pulsion fluids. The corresponding capacity factors follow 
from Fig. 2 through Fig. 5. Of course, for final conclusions 
the assigned net weight factor must be compared with the 
project net weight factor to determine the technical feasibility 
for a given application. 

However, it has been seen in the previous sections that the 
net weight factor decreases with increasing capacity factor. 
Therefore, knowledge of assigned net weight factor and cor- 
related capacity factor permits at least estimating the relative 
merits and the range of favorable application of different 
fluids. This shall be illustrated by a few examples. 

A sounding rocket, ascending vertically at 2 g absolute ini- 
tial acceleration and with a payload factor of 0.05, shall attain 
a cut-off speed of 7000 fps. The propulsion fluids considered 
are listed in Table 3. The upper part of this table shows the 
required capacity factors from Figs. 2 and 3, and from Equa- 
tion [12], the velocity parameters v,/g/sp, using Table 2 (ex- 
pansion ratio 20 to 0.7 or 1 atm, respectively) and finally the 
assigned net weight factor according to Fig. 8. Again, the 
application of light fluids appears difficult, since they require 
capacity factors which are considerably larger than those of 
heavier fluids, while their assigned net weight factors are not 
increased in preportion to offset this effect. Moreover, for 
vehicles using the dense propellants, a design analysis can be 
expected to yield project net weight factors which are smaller 
than the assigned values given in Table 3. 

Comparing, for instance, (O2-H2),. with O3-CH,, it follows 
that the first propellant yields a vehicle about twice as large 
at half the mean fluid density of the ozone combination; but 
practically the same net weight factor is assigned to both ve- 
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; ie For liquefied gases the density at 1 atm vapor pressure was used. 


hicles. It is very obvious in this case that O;-CH, would be 
preferable from the viewpoint of vehicle performance. Also, 
in comparison to the other propellants, oxygen-hydrogen 
would be inferior in the present example, even if its assigned 
net weight factor would technically be feasible, which is doubt- 
ful. 

Comparing the propellants with the working fluids, it is 
rather obvious that, in this performance range, they cannot 
compete with the propellants. This applies to hydrogen in 
particular, since it appears very unlikely, in view of the sup- 
posedly heavy atomic propulsion system, that a net weight 
factor of 0.51 could be realized at all. Ammonia and hydra- 
zine would impose less severe problems from the design view- 
point, but also in their case the assigned net weight factors 
seem inadequate for actual construction. Among the heavier 
propellants, ozone-methane would be outstanding, since its 
density is but slightly lower than that of oxygen-hydrazine, 
while its assigned net weight factor is comparatively much 
larger. The difference between the oxygen ana nitric acid 
combination with hydrazine is small at this performance level. 
Preference’ of either one would be influenced by other factors. 

As a second example, assume a space ship capable of at- 
taining an ideal velocity of 20,000 fps with a payload factor of 


TABLE 2 CHARACTERISTIC DATA OF PROPULSION FLUIDS CONSIDERED 


Expansion Mean Exhaust —Component specific— 
ratio specific velocity gravity * 
Propulsion fluid (atm) gravity (fps) Oxid. Fuel 
1 Hydrogen 0.071 23,902 
28 , 687 
20/0.1 17,185 
3 0.68 12,272 
4 Oxygen-hydrogen 20/0.7 0.266 11,914 1.14 0.071 
(max. spec. impulse) 20/0.1 0.30 13,556 
4a Oxygen-hydrogen ibe 20/0.7 0.424 10,110 
(stoichiometric) 3 20/0.1 0.424 12,560 
5 Hydrazine : x 20/1 1.01 10,507 
1 12,487 
6 Ozone-methane ep 20/0.7 0.96 9,963 0.424 
20/0.1 0.97 11,946 
7 Oxygen-hydrazine 1.065 9,177 1.14 1.01 
i 1.07 10,819 
8 Oxygen-ammonia ahs 20/0.7 0.88 8,791 1.14 0.68 
20/01 0.90 10,401 
9 Tetranitromethane- we 20/0.7 1.32 8,501 1.65 1.01 
hydrazine 20/0.1 1.33 9,982 
10 RFNA-hydrazine 1.276 8,372 1.57 1.01 
(RFNA: 85HNO;- 20/0.1 1.28 9,821 
15NO.) | 
Chlorine trifluoride- 1.4 8,275 1.76 1.01 
hydrazine 20/0.1 4 9,467 
Water 1.0 7,190 
Nitrogen 81 6,46 
Ca 20/0.1 7,651 
20/0.1 7,293 
Argon al 20/1 1.59 4,852 
20/0.1 5,438 
Mercury 40/1 13.546 2,418 


0.10. The same fluids shall be compared. Their relevant 
data are listed in the lower portion of Table 3. The picture 
has changed greatly, as compared to the previous example. 
The increase in capacity factor for hydrogen is very small, 
since no gravity losses occur in this case. The net weight 
factor of the oxygen-hydrogen ships is now about one third 
and one fourth, respectively, of that of the hydrogen ship, 
while in the previous example they were better than one half 
of this value. The change is due to inferior specific impulse 
and to a larger payload factor. Nevertheless, the higher 
density of oxygen-hydrogen and the smaller propulsion sys- 
tem weight would make this propellant at least equivalent to 
hydrogen, let alone greater simplicity of design. The data 
also show that oxygen-hydrogen could compete with nitric- 
acid-hydrazine in this performance range. This was not the 
case in the previous example. However, it still would be in- 
ferior to ozone-methane and to oxygen-hydrazine, leaving 
these medium propellants in the highest bracket. As in the 
previous example, ozone-methane appears to be the best pro- 
pellant because of its favorable combination of specific im- 
pulse and density. Oxygen-hydrazine and oxygen-hydrogen 
have gained in importance while nitric-acid-hydrazine has 
lost its high rating. Among the working fluids, ammonia and 


TABLE 3 COMPARISON OF PROPULSION FLUIDS 
v1 (fps): 7000 (W/W. = 0.05) 
2 5.3 4.13 1.66 1.85 1.89 1.77 
0.571 0.59 0.691 0.666 0.705 0.765 0.836 
0.394 0.38 0.33 0.345 0.325 0.297  ~— 0.27 
(fps): (Wi/W, = 0.10) 
¢ (cu m/sec) 14.25 4.6 11.4 3.95 4.5 . i 
v:/gIsp (—) 0.695 1.41 1.48 oe 1.605 1.68 1.84 2.03 
fw (-) 0.447 0.16 0.14 0.113 0.11 0.095 0.065 < 0.05 
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hydrazine probably are superior to hydrogen, in spite of lower 
assigned net weight factors, because of much higher density. 
Generally the position of the working fluids with regard to the 
chemical propellants is greatly improved in comparison to the 
previous example. Already the data of the heavy working 
fluids appear competitive when compared with nitric-acid- 
hydrazine or even oxygen-hydrazine, disregarding develop- 
ment problems, More detailed specifications depend on an 
exhaustive comparison of atomic and chemical propulsion sys- 
tem weight requirements, but the general tendency is 
apparent. 

If interplanetary flight shall ever rise beyond the status of 
occasional scientific expeditions, it seems necessary, for in- 
stance, to exploit valuable deposits on other celestial bodies, 
preferably on satellites and asteroids which could be contacted 
directly by space ships. Their payload factor must be as high 
as possible for economic exploitation. If possible, no auxiliary 
ships should be needed, to avoid loss of costly engines, particu- 
larly if atomic piles are used. This means that still higher 
ideal velocities must be required for the individual space 
ship. For a trip to the inner moon of Mars, e.g., which 
requires about the lowest of all ideal velocities for inter- 
planetary flights, the over-all value would be about 41,000 fps 
when leaving from an artificial satellite. Assume, as a third 
example, this ideal velocity and a payload factor of 0.2. In- 
spection of Fig. 8 shows that even a hydrogen atomic pile 
ship is allowed only a net weight factor of about 0.055 for ac- 
complishing this performance. For all the other fluids con- 
sidered, the net weight factor would be zero already, long 
before this ideal velocity could be attained, even at smaller 
payload factors. This means, if such a ship can ever be built 
at all, the working fluid would have to be hydrogen. 

These three examples verify what was previously revealed 
tentatively in Figs. 2 through 5, namely, that high density is 
important and necessary for vehicles of comparatively low 
cut-off speeds. As the mechanical performance is increased, 
specific impulse gradually becomes more important and ulti- 
mately emerges as the decisive parameter. 

Solving Equation [11] for fy, the assigned net weight factor 
for a certain payload factor can be found if the loading factor 
or the capacity factor is determined for a given value of v/gI sp 
from Equation [4] or [8]. Thus the significance of Fig. 8 can 
be specified for any desired propulsion fluid, as it has been 
done in the afore-mentioned three examples for a few selected 
fluids. The result can be plotted as net weight factor versus 
cut-off velocity. Such a graph is shown in Fig. 9 for vertical 
ascent and for a payload factor of 0.10. The highest net 
weight factors are assigned to the lightest fluids, since they 
yield the smallest loading factor. Those portions of the 
curves which are represented by solid lines indicate the favora- 
ble range of application; that is, the range in which the as- 
signed net weight factor is likely to satisfy the technical re- 
quirements. Within the performance range covered in Fig. 
9, the upper four fluids have little chance for favorable ap- 
plication. An atomic rocket using hydrazine might become 
applicable in the upper range of cut-off speeds considered in 
Fig. 9. For speeds between 7000 and 10,000 fps, ozone would 
be a favorable propellant A similar range applies for oxy- 
gen-hydrazine, while the heavier propellants may be applied 
more favorably in a somewhat lower range of speeds. Due 
to lower specific impulse and density and because of the pre- 
sumably heavier weight of the atomic pile unit, working fluids 
such as water and nitrogen would require net weight factors of 
the order of 0.4 at least. This, in principle, would make these 
fluids eligible for a lower speed range between 2000 and 6000 
fps. Actually, of course, this range is occupied by technically 
more preferable propellants such as oxygen-alcohol] or acid- 
aniline, which are not considered in the graph. The specific 
impulse of nitrogen is similar to that of these two propellants 
which, however, have higher density and can be used in pro- 
pulsion systems already developed. Mercury theoretically 
would fall in a similar velocity range. Its net weight factor 


SEPTEMBER-OcTOBER 1953 


would be extremely small. Fig. 9 shows that, with the pos- 
sible exception of hydrazine as working fluid, atomic propul- 
sion systems, in fact, would be of no advantage in this per- 
formance range, since chemical propellants are available which 
do as good or better. If may be emphasized that the solid 
portions of the curves are only very rough indications of 
favorable ranges of application and do not imply accurate 
limitations. 

For velocities exceeding those considered in Fig. 9, either 
one of the upper five fluids would become eligible or, if heavier 
propellants are used, the vehicle must consist of more than 
one stage. 

Considering the propulsion fluids as density-impulse points 
(Fig. 1), one can visualize these points sliding down on the net 
weight factor scale as the performance is increased. All 
points are high up when the vehicle performance is small. 
They cannot rise above fy = 1, hence must be crowded at the 
upper end of the scale. Many fluids are available for seleec- 
tion to suit specific projects. In general, that fluid is most 
favorable, from the viewpoint of performance, which, for a 
given gross weight, yields the highest ratio of assigned net 
weight factor over project net weight factor. With increasing 
vehicle performance, the points slide down on the net weight 
factor scale and begin to spread, since the light fluids with 
high specific impulse drop less rapidly. Since there is a lower 
technical limit for the net weight factors connected with each 
fluid (namely, when the ratio of assigned over project net 
weight factor becomes less than one), the heavy propellants 
are the first to drop below this limit. They become princi- 
pally inapplicable when they pass this limit at zero payload 
factor. This selective effect gradually leaves a higher per- 
centage of light fluids for practical application as the perform- 
ance is increased. Eventually only hydrogen remains, and 
thereafter either the multistage principle must be applied or 
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FIG. 10 NEW WEIGHT FACTOR FOR DIFFERENT IDEAL VELOCITIES 


different methods of propulsion (for instance, high-speed 
ion jets) must be considered. 

The process described before is illustrated in Fig 10 for 
space flight conditions. The fluids are designated by num- 
bers in accordance with Table 2. In order to show the in- 
fluence of the payload factor, three values, 0.05, 0.10, and 
0.20, are represented. The assigned net weight factors cor- 
responding to zero payload factor are indicated for several 
fluids. These points determine the highest possible loading 
factor for the space ship using the particular fluid at the given 
ideal velocity. 

It can be seen that now, in contrast to Fig. 9, some working 
fluids, particularly hydrazine, ammonia, helium, and finally 
hydrogen become considerably superior to all propellants, 
including oxygen-hydrogen. While there is a wide variety of 
fluids available at 10,000 fps, their number is reduced con- 
siderably at 20,000 fps. At this ideal velocity, ammonia and 
hydrazine, at a payload factor of 0.1, possess the highest den- 
sity among the fluids shown’ in the net weight factor range 
between 0.1 and 0.3. Ozone-methane and oxygen-hydrogen 
appear as favorable propellants and also oxygen-hydrazine 
may still be applicable. Between 20,000 and 30,000 fps, a 
range which comprises the energy requirements for circum- 
navigation of the moon when departing from an artificial 
satellite, most propellants cease to apply with the possible ex- 
ception of ozone-methane and oxygen-hydrogen for small 
payload factors. At 30,000 fps, practically hydrogen and 
helium alone are available for payload factors above 0.05. 
At 40,000 fps, hydrogen alone can accommodate payload 
factors above 0.05. In other words, unless the payload factor 
is very small (Wi/Wo < 0.05), chemical propellants apply 
for ideal velocities up to about 20,000 fps to be attained by a 
single ship or stage. In the region between 20,000 and 
40,000 fps, working fluids definitely become superior to chemi- 
cal propellants. Particularly between 30,000 : and 40, 7 fps, 


only working fluids allow for some net weight at all if the pay- 
load factor isabove zero. Depending on the weight per unit 
thrust ultimately obtainable with thermonuclear propulsion 
systems, hydrazine and ammonia would be very promising 
working fluids also for terrestrial use in single or two-stage 
vehicles designed to attain cut-off speeds of the order of 
15,000 to 20,000 fps. 

Finally, it may be emphasized that the speed range between 
35,000 and 41,000 fps represents about the lower speed limit 
required for interplanetary voyages (Mars) when taking off 
from an artificial satellite. Thus it is very unlikely that, in 
the foreseeable future, a single space ship can reach other 
planets in our solar system. However, by appropriate selec- 
tion of propellants and especially by developing atomic pro- 
pulsion, the number of supply ships could be reduced to « 
minimum, therewith simplifying space flight substantially. 


Conclusions 


The procedure outlined herein permits a rapid determin:- 
tion of the relative merits and of the range of applicability of 
propulsion fluids in rockets for different cut-off speeds and 
payload factors. This is accomplished by evaluating the 
characteristic capacity factor and net weight factor assigned 
to each fluid at given cut-off speed and payload factor. Know- 
ing these assigned values, it can be determined by design 
analysis whether or not, for a given project, their realization 
is technically feasible. 

The method has been applied to a general investigation of 
chemical propellants and working fluids for atomic propulsion 
over a wide range of vehicle performance. It was found that 


up to 10,000 fps cut-off velocity, that is, for most terrestrial 


applications as well as for satellite operations with multistage 
vehicles, a large number of propulsion fluids are available. 
Among them, working fluids are inferior to propellants or at 
least do not offer practical advantages, unless the payload 
factor is very high. Among the propellants, heavy and me- 
dium combinations appear preferable. If higher speeds are 
to be attained with single or two-stage vehicles, thermo- 
nuclear propulsion systems, preferably using hydrazine or 
ammonia, are advantageous. 

For space flight, chemical propellants have competitive 
qualities, with emphasis on high specific impulse rather than 
density, up to ideal velocities of about 20,000 fps or somewhat 
higher, if payload factors of 0.05 or below are accepted. A- 
side from technical considerations, their range of applicability 
is extended, because, in space flight, the initial acceleration 
can be kept small so that high loading factors become feasible 
without yielding cut-off accelerations of physiologically pro- 
hibitive magnitude. The propellants applicable for these 
high ideal velocities are ozone-methane, oxygen-hydrogen, or 
similar combinations, and perhaps oxygen-hydrazine. Be- 
yond 20,000 fps, working fluids become increasingly superior, 
particularly at payload factors greater than 0.05. Above 30,- 
000 fps they are the only ones that appear technically applica- 
ble at all. Beyond ideal velocities of about 40,000 fps, ther- 
mal propulsion methods cease to apply unless the ship is re- 
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A method for establishing the relative desirability of 
competing rocket power plant components is developed 
for long-range, two-stage rocket-powered vehicles. The 
measure of the desirability proposed is minimum net cost 
wherein missiles of the same class are required to perform 
identical functions; e.g., same range carrying the same 
payload. The relationships developed are derived from 
simple missile performance parameters and lend them- 
selves to readily calculable results. 

Characteristics of a hypothetical two-stage, long-range 
rocket vehicle are presented, and for this missile exchange 
rate curves are given which show the order of compromise 
allowable between weight and efficiency for any power 
plant component or between competing power plants. 
Several power plant arrangements are analyzed for this 
missile, and the proposed technique for establishing desira- 
bility is demonstrated. 


Nomenclature 

AF = difference in propellant weight, lb 
AF,» = difference in first-stage propellant weight, |b 
I = effective specific impulse, sec 
m = mass ratio 
n = portion of AT that is jettisoned 7 
Ap = difference in component weight, Ib 
AP; = difference in jettisoned component weight, |b 
AR;;, , = difference in propellant flow rate required by compo- 

nent in second stage, lb/sec 7 
AT = difference in tank weight, lb eine 
ts_4 = second-stage burning time, sec = 
AV = difference in tank volume, cu ft 
Wo = gross weight at take-off, Ib tae. 
We = weight at burnout, lb & 
W. = weight at end of first stage, lb ., 
W; = weight at start of second stage, lb 
Wy.» = weight of propellant consumed during first stage, Ib 
a = ratio of AT’ to AV 
p = propellant density, lb/cu ft 


Introduction 


NE of the major recurrent problems in designing a 

rocket-powered vehicle is the selection of individual com- 
ponents to make up the power plant. In choosing a com- 
ponent, the missile designer gives consideration to availa- 
bility, reliability, weight, efficiency, and cost. The first two 
of these yardsticks are dependent on market conditions and 
state of the art of the component in question and are outside 
the scope of this paper. If we limit consideration to quantity- 
produced vehicles, it would seem that the most sensible step 
would be to select the component which can perform its func- 
tion satisfactorily for the minimum cost without affecting the 
basic missile function. Cost, in this sense, involves not only 
the price of the component but also the cost penalty of its 
weight. Some relationship is required which will describe the 
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rate of exchange between component weight and efficiency as 
a function of the specific weight and performance characteris- 
tics of any missile. This relationship should also include 
terms to evaluate the penalties incurred if other than even 
exchanges between weight and performance are considered. 
The penalty could be considered in terms of extra propellant 
and structure required to permit the missile to carry out its 
design function when operating with heavier or less efficient 
components. 

The general approach of the work presented in this paper 
is to consider the effect small changes in component weight 
and operating efficiency have on required propellant load 
when maintaining missile range and payload weight. It is 
assumed that for long-range vehicles, where the drag is small 
relative to the weight, a reasonable, though approximate 
measure of the range, for constant initial acceleration, is the 
parameter J In m, where J is the average effective specific im- 
pulse, and m is the mass ratio. On this basis, for any given 
missile to maintain constant range, J In m must be held con- 
stant. If competing components having different operating 
efficiencies are considered, they will produce a change in the 
effective impulse and, in order to keep the J In m parameter 
constant, the mass ratio will have to be adjusted. To some 
extent, it will be immediately adjusted by the change in com- 
ponent weight; but in order to balance the range parameter 
completely, some change in the weight of propellant carried 
will probably have to be made. If the costs of propellant and 
component per pound were known, the net costs of component 
and propellant differences can be compared to establish the 
cheapest system to carry the same payload the same distance. 
If the propellant weight change is greater than the component 
weight change, there will be a difference in the initial accel- 
eration. For the comparatively small component variations 
considered in this paper, the results are not greatly affected. 


Development of General Relationships 


Two basically similar missiles, denoted by the subscripts A 
and B, will be considered. These missiles will be identical in 
all respects except for a power plant component which will 
differ in both weight and efficiency, where component effi- 
ciency is measured in terms of the over-all specific impulse. 
Both missiles will be two-staged and be required to carry the 
same payload over the same range. 

To maintain the same range 


I,4in ma = const = Ig 
which can be rewritten as 


By definition 
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where Wo, = gross weight at take-off 
W. = weight at end of first stage =e 
W; = weight at beginning of second stage “Fs 
We = weight at burnout 


The difference between W: and W; is the jettisoned 


weight, W;, so Equation [3] can be written as 


Combining Equations [2] and [4] 


=F mp = = oB [1 

In this case, the pertinent differences between missiles A 
and B are: (1) Some change in component weight AP and a 
corresponding change in jettisoned weight, AP;; (2) some 
change in propellant weight, AF, required to operate the 
extra component weight and also to maintain the J In m 
parameter constant; and (3) some change in tank and struc- 
tural weight required to carry AF. 

To evaluate this latter quantity, let AZ’ represent the 
change in tank weight and structure incurred by AF, and as- 
sume that AT is proportional to the change in tank volume, 
AV. Then 


where p = propellant density. Let gee 


On this basis, missile B weights can be expressed in terms of 
missile A weights in the following manner 


Wos = Wor + AP + AF + [3] 

Ws, = Wa, + (AP — AP;) + — AF......... [9] 
Wi, = Wi, + + — (10) 

Wey = Woy — (Wig, + 


= Wo, +4P + (1 + *) (Wry, + AFi-2).. [11] 


Substituting Equations [8], [9], [10], and [11] into [4] 


yields 
Wo, + AP + aF (1 +) 
mg = 
We, + AP — AP; +“ aF 
4 p 
Wi, + AP; + “* aF 
=> AF AF3-, = AF , [13b] 


where AR, ,_, represents the difference in propellant consump- 
tion rate of the two missiles, due to the component change, and 
t;_, equals the second-stage burning time. In Equation [13b] 
it has been assumed that all of the extra propellant in excess of 
that required to operate AP has been consumed in the first 
stage. Although this is not exact, it greatly simplifies the 
analysis and the resultant error is quite small. 

Substituting Equation [13b] into [12] and equating the 


298 


tion. 


AF 
p 


Wo, 
A A 
na 


Wa, 
Wi, + AP; +7, 


We, + aP + ar ( 


Ws, + AP — AP; + 


1- 


Wo, + AP AF — Wy, + ARys-e 


[Wo 


na 


Vi, + + aF 


4. AF + ARys-4 ts-4 


(Wo, + AP)) . [15] 
1+ 
= 
ef p 
[21] 


Substituting Equations [16] through [22] into [15], the re- 
sulting relationship can be rewritten as 


AF? + Bo) + Bo) CoBo?| + AF [Eo (1 Bo) 
G + Bo) + Ey By Do (1 + By) + By’ Ey CoAoBo 
CoBoGo] + [ Eo? + + Co Ao0Go] = 0. [23] 


This is in the form of a simple quadratic equation, the 
solution of which can be written as 


AF = 
where 


A; = Bo(1 + Bo) + Bo’ (1 + Bo) — CoBo? 


B, = Eo(1 + Bo) + Go(1 + Bo) + EoBo — Do(1 + Bo) + Bo’ Eo— 
CoByEo Cy BoGo rrr [26] 


[27] 


B, (Bi? 4A,C;)'/2 


C; = E,? + EvGo DoEo CoAnEo Co 


Equation [24] essentially provides an exchange rate be- 
tween component weight and efficiency as a function of mis- 
sile performance. With the differences in component and re- 
quired propellant weight of the proposed systems established, 
a cost comparison could be performed if dollars per pound of 
component and propellant were known. The most desirable 
system would be that whose combined cost of component and 
propellant to operate and carry the component is the least. 

Detailed examination of the terms in Equation [24] in- 
dicates that for small changes in AF the terms representing 
tankage and structure incurred by changes in propellant load 
can be neglected with only slight effect on the resultant solu- 
If this were done, Equation [12] could be rewritten as 
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(WPo, + AP). .|29] 


Wo, + AP + AF 


We, + (AP — AP)j) 
Wig + AP; 
Wo, + AP AR, 
and following the same procedure as before 
(Wa, + AP OP 
AF = Wj, + AP; 
Wo FAP = Wi, FAR, 


Because of its relative simplicity, 
used throughout the remainder of this paper. 


Equation [29] will be 
It can be seen, 


change rate between AF, AP, and J4/Iz is a function of the 
weight and performance characteristics of the missile in 
question. The “best”? component and/or arrangement of 
components is not necessarily the lightest and, similarly, not 
necessarily the most efficient. Rather, a compromise is re- 
quired, permitting the missile to perform its designated func- 
tion reliably, at the minimum total operating cost. As might 
be expected, Equations [24] and [29] do not, in themselves, 
provide the basis for an absolute judgment as to component 
selection. The information they do provide is in the form of 
a measure of desirability which, when coupled with knowledge 
of relative cost, availability, and reliability, will provide a 
basis for choosing the most suitable configuration. These — 
relationships may also be used to advantage when comparing — 
components which do not affect the propellant consumption — 
because of their efficiency characteristics but do present a 


from this relationship, that the actual magnitude of the ex-_ 
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problem in that they differ in weight and have varying por- 

tions of their weight jettisoned. In this case, J4/Iz is equal 

to unity, and for the various combinations of AP and AP; 

in question, the resultant AF’s could be established. Re- 
lationships such as these can be developed, with relatively 
small effort, for missiles having any number of stages. 

To illustrate the application of these relationships, con- 
sideration will be given to a hypothetical long range ballis- 
tic missile. Fig 1 is a sketch of this vehicle and Table | lists 

some of its pertinent characteristics. The two stages are 


TABLE | WEIGHTS AND DIMENSIONS OF SAMPLE MISSILE 


Weights 
Ib 


Missile 
Warhead (payload of 1000 Ib) 2,050 
Instrument compartment 1,350 
Center section 5,385 
Propulsion unit 2,775 
Aft structure 1,116 
Burnout weight 12,676 
Booster 
Tank section (5) 6, 805 
-_ Propulsion unit (5) 9,750 
structure (5) (including jettisoning mec — 
and structure) 4,168 
Empty weight 33 , 399 
Launching weight 269 , 000 
Total propellant weight 235 , 600 
Stage 1 propellant ee 188 , 480 
Stage 2 propellant 47,120 
Weight at end of stage 1 80,520 
Dimensions 
Over-all length, ft =: > 82 
Missile body diameter, in. t. = 60 
Booster length, ft : 64 
Booster body diameter (each of 5), in. 60 


parallel mounted, with five V-2 power plants, used in the 
boosters, set around the aft circumference of the missile. 
Detailed power plant data were obtained from Sutton (1)? and 
Kooy and Uytenbogaart (2). No attempt was made to design 
the best possible missile, and such detailed weight and perfor- 
mance changes as were made were done more to demonstrate 
the use of the techniques described in this paper than to ad- 
here closely to realistic values. The missile will have a ballis- 
tic trajectory, a range of approximately 1750 nautical miles, 
carry a payload of 1000 Ib, and have an over-all mass ratio of 
15.8. Propellant combination, sea-level impulse, and general 
arrangement will be identical to those used on the V-2. 
Substituting values from Table 1 into Equation [29], a 
curve sheet such as Fig. 2 can be prepared showing the varia- 
tion of AF with AP using J4/I, as a parameter. The inter- 
septs and slopes of these curves are very much dependent on 
the weights and performance of the missile in question, and it 
follows that the exchange rate level for any component will be 
established quantitatively by these basic characteristics. 
For the case of zero efficiency change, where /4/Iz is unity, a 
decrease in component weight is reflected by a corresponding 
decrease in the propellant load, corresponding with the mass 
ratio. As #, the ratio of AP; to AP, is increased, the pro- 
pellant weight-saving is seen to decrease. If component 
efficiency is decreased, producing an increase in I4/Iz, the 
general level of the exchange rate curves is raised accord- 
ingly and greater component weight reductions are re- 
quired to effect a comparable propellant weight-saving. 
Increases of range and performance requirements of the mis- 
sile also tend to raise the level of the exchange rate curves. 


2 Numbers in parentheses refer to the References listed on 
page 300. 
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Fig. 3 is a plot of AP and J4/Js for the case of zero change in 
propellant weight, and indicates the component weight- 
saving required to compensate for reductions in efficiency 
while maintaining constant propellant load. 

To illustrate the application of these curves, two selection 
problems that might arise will be investigated. The nu- 
merical values used in these examples were chosen primarily 
for their illustrative quality and are neither representative of 
any special design nor the result of detailed study aimed at 
modifying the power plant. 

Consider first, the possibility of modifying the pumping 
plant. The original system uses six turbopumps weighing 
350 lb each, five of which are jettisoned at the end of the first 
stage. As opposed to this, the proposed system will be 
lighter, each turbopump weighing only 300 lb, but will re- 
quire 0.50 per cent more propellant per sec. Comparing 
these data yields AP = —300 lb, AP; = — 250 lb, = 0.834, 
and J4/Ig = 1.005, where subscripts A and B represent the 
original and proposed systems, respectively. From Fig. 2, 
AF = +1500 lb, which represents a 1500-Ib increase in the 
quantity of propellant required if the new system is used. It 
can be seen from Fig. 3 that for the same component weight- 
saving, the minimum allowable I4/J, to maintain constant 
propellant weight is 1.0022. This is approximately equiva- 
lent to an increase in required propellant flow rate of 0.25 per 
cent. At first inspection, it would seem that the original sys- 
tem is best. However, considering production of a large 
number of vehicles, the actual selection should be made on a 
cost-per-missile basis, since it is apparent that with relatively 
minor modifications in structure, either configuration could be 
used without preventing the missile from completing its task. 
It will be arbitrarily assumed that machinery costs 55 dollars 
per pound, structure costs 25 dollars per pound, and propel- 
lant, 0.07 dollar per pound. From Equation [6] we know 
that AT ='a/p-AP,and for the purposes of this example we 
can assume a/p= 0.043. Using these values, the saving in 
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component cost is 16,500 dollars, while the cost of extra pro- 
pellant and structure is 1717 dollars. If the assumed costs 
per pound are correct, the proposed system would be an 
economic improvement. 

Another possible selection problem might be modification 
of the air bottle arrangement. The original system has an 
air bottle weight of 2070 Ib and jettisons 1385 lb at the end of 
the first stage. The possibility will be considered of replacing 
this arrangement with one that has an initial weight of 500 |b 
more but jettisons 2000 lb, thereby reducing the final burnout 
weight. Using the proposed arrangment as a base, AP = 
—500 Ib, AP; = —615 lb, @ = 1.23, and since no efficiency 
change is involved, J4/Ig is zero. From Fig. 2, AF = +1000 
Ib, indicating that the proposed system would permit a 1000- 
Ib propellant reduction. As before, the selection decision 
will be affected by the individual costs in question. 
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Nomenclature 


PT = dimensionless instantaneous value of the time lag re- 
duced by the characteristic time 6/2 = 7; + 7 
7, 7; = dimensionless steady-state value and instantaneous 


value of the part of the total time lag which is in- 
sensitive to the pressure oscillation 

7,7 = dimensionless steady-state value and instantaneous 
value of the other part of the total time lag which 
is sensitive to the pressure oscillation 


n = exponent in the pressure dependence of the time lag = 
pressure index of interaction 

m; = rate of injection of the propellant per unit cross- 
sectional area of the combustion chamber 

rity = rate of the generation of the hot gas from the combus- 


tion of the propellant per unit cross-sectional area 
of the combustion chamber 

L = combustion chamber length from injector end to the 
entrance of the de Laval nozzle = reference length 
scale for nondimensionalization 


z° = distance from injector end. along the combustion 
chamber 

t* = dimensional time 

u* = dimensional mean flow speed of the gas along the 
combustion chamber axis 

p*, p*, T7*, andc* = dimensional instantaneous values of 


pressure, density, temperature, and speed of sound 
in the burned hot gas 


T*o, and c*> = values of p*, p*, and c* at injector 


end 
= reference quantities for nondimensional- 

ization 

p*, o*, T*, andc* = steady-state values of p*, p*, 7'*, and c* 

6/2 = L/c* = characteristic time = time required for 
a sound wave to travel the entire length of the com- 
bustion chamber filled with stagnant burned gas 

t = t*/(L/c*>) = dimensionless time 

z = 2°/L = dimensionless length 

ui = u*/c%, = dimensionless velocity of the gas in 

_ unsteady and steady-state operation 

M, Af. = Mach number of the gas flow in unsteady and steady- 
state operation 

p, p, T, ande = dimensionless instantaneous values of 


pressure, density, temperature, and speed of sound 


p, 6, T,and@ = dimensionless steady-state values of p, 


p, T, and c. 

T’, and .c’ = dimensionless instantaneous perturba- 
tions over their respective steady state values 

a = +i = root of the characteristic equation with 
the dimensionless time as the independent variable 

rN = dimensionless amplification coefficient 

= dimensionless angular frequency 
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Combustion Instability Rocket Motor Motor 


2 = absolute value of the angular frequency = (0/2) 
7 = adiabatic index of the combustion gas ; 
= p’ exp (—at) 
6 = p’ exp (—at) 
= wu’ exp (—at) 
é = dimensionless distance of the concentrated combus- 


tion front from the injector end expressed as a 
fraction of the characteristic length L 
8 = reduced angular frequency of the oscillation 
= angular frequency divided by the velocity gradient 
in de Laval nozzle 


z = reduced velocity parameter = #(y + 1)/2 
v/t 
I(6, %) = [ ; A = the ratio of fractional variation of veloc- 
ity to fractional variation of density at combustion 
chamber exit or entrance to de Laval nozzle 
k,S = real and imaginary parts of J(8, 7) 
h,k = integers characterizing the modes of the oscillation 
D/Dt = substantial derivatives along the path of a propellant 


element 
O( ) = the order of magnitude of the quantity in the bracket 


Subscripts 

xort = partial differentiation with respect to x or t 
lor 2 = the quantities evaluated in the flow field 1 or 2 


1 Introduction 


OUGH combustion as a result of large pressure oscilla- 
tions in the combustion chamber of a liquid propellant 
rocket motor has been observed under different circumstances 
in two distinct ranges of frequencies: the low-frequency range 
of less than 100 eps, and the high-frequency range of several 
hundreds or several thousands cps. Such rough combustion 
not only gives fluctuating performance but also shortens the 
life of the rocket motor. An understanding of the basic 
mechanism of producing unstable pressure oscillations that 
lead to rough combustion is therefore of great practical im- 
portance. 

It has long been recognized that for unstable operations 
the oscillation of the chamber pressure and the oscillation of 
the rate of hot gas generation produced by the pressure oscil- 
lation must be properly out of time phase so that an increase 
of the rate of hot gas generation occurs at an overpressure 
period and further increases the overpressure in the combus- 
tion chamber. This time phase difference is originated from 
the fact that the propellant element does not burn imme- 
diately after being injected into the combustion chamber, but 
burns after a certain time interval, called the “time lag,” 
during which the fuel and the oxidizer particles mix properly 
and absorb the necessary amount of activation energy. 
In (1, 2, 3),4 the low-frequency oscillation has been analyzed 
based on the assumption that the time lag is constant and 
independent of the oscillations of the gas system in the com- 
bustion chamber. In these analyses, a pressure-sensitive 
feeding system which prov ides a varying injection rate under 
the pressure oscillation in the combustion chamber is assumed 
to be the self-exciting mechanism which creates the variation 


Numbers in parentheses refer to the References on page 313. 
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of the rate of hot gas generation. The senior author of the 
present paper has pointed out in (4) that the pressure sensi- 
tivity of the feeding system is not a necessary self-exciting 
mechanism for producing unstable pressure oscillations. By 
assuming a varying time lag 7 which depends on the chamber 


it is shown that unstable oscillations of beth the low-frequency 
range and the high-frequency range can be produced even if 
the injection rate is constant. This kind of combustion in- 
stability cannot be eliminated by properly designing the feed- 
ing system but is intrinsic in the nature of the combustion 
processes. 

The case of low-frequency intrinsic combustion instability 
is extensively studied by the senior author in (4), while the 
high-frequency case is only briefly discussed, using a simpli- 
fied model of a single concentrated combustion front near the 
injector end and for the particular value of n = 1/y where 
y is the adiabatic index of the burning gas. As n represents 
the extent of interaction between the pressure oscillation and 
the combustion processes, n must be an important parameter 
as a stability criterion. In practical cases, a large part of the 
combustion often takes place in a narrow region somewhere in 
the combustion chamber but usually not too close to the in- 
jector end. Therefore the analysis in (4) is extended in the 
present paper to the case with arbitrary location of the con- 
centrated combustion front and for arbitrary values of n. 
The effect of distributing the combustion along the combus- 
tion chamber axis presents a more difficult problem and is 
analyzed in (5). The problem of shifting part of the concen- 
trated combustion at the injector end to arbitrary axial loca- 
tion is briefly studied in the present paper with a view to some 
indication of the effect of distributing the combustion axially. 

The unsteady supercritical flow in a de Laval nozzle with 
linear steady-state velocity profile in the subsonic portion 
has been studied in (6) and (7) from which the boundary con- 
dition for high-frequency oscillations can be deduced. This 
boundary condition will be used in several specific examples. 
In the limiting case of a very short nozzle, this boundary con- 
dition is equivalent to the boundary condition of constant 
flow Mach number at the entrance of the nozzle which is used 
in (4). This short nozzle boundary condition is simple and 
admits analytical solution of the characteristic value prob- 
lem. The results with these two boundary conditions will 
be compared. 


2 Formulation of the Problem 


2.1 Assumptions and Simplified Models of Gas Flow Sys- 
tem 


The combustion chamber of a liquid propellant rocket 
motor is often a straight duct of constant cross-sectional 
area and is filled with the hot burned gas. The propellant 
elements injected into the combustion chamber are mostly 
in the form of atomized liquid droplets suspended in and 
carried along by the burned gas stream without occupying 
appr-ciable volume of the combustion chamber. The hot 
gas generated from combustion is recirculated actively to the 
region near the injector end where the hot gas supplies the 
activation energy to the unburned propellant elements. The 
recirculation and the mean flow patterns of the hot gas are ex- 
tremely complicated, depending largely on the design and 
arrangement of the injectors. For the present analysis, we 
shall consider the gas flow inside the combustion chamber to 
be a one-dimensional flow of the hot gas only, with the un- 
burned propellant elements suspended in and carried along 
by the hot gas. Oscillations in the transversal plane normal 
to the chamber axis is not being considered. The hot gases 
all over the combustion chamber are generated from the 
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combustion of the same propellant under slightly different 
pressure, and thus are essentially at the same stagnation tem- 
perature. The combustion process in a liquid propellant 
rocket motor does not primarily increase the specific energy of 
the gas system. The heat released by the combustion of a 
propellant element is used to raise the temperature of the 
combustion products of this element to the temperature 


of the burned gas system. Therefore, from the point of 
view of the flow of the burned gas, the combustion proc- 
ess is essentially a process of generating or introducing 
new mass of the burned gas into the flow system, while the 
specific energy of the gas system remains substantially con- 
stant. As the combustion chamber is a duct of const:nt 
area, the mean flow velocity of the burned gas must incre:ise 
whenever additional mass of burned gas is introduced into 
the flow system. The ‘new’ burned gas and the “old” 
burned gas are assumed to mix intimately and accelerate to- 
gether, carrying with them the suspended unburned propel- 
lant elements. This process of mixing and the process of 
accelerating the suspended particles give rise to an entropy 
variation in the gas flow field even if the specific energy re- 
lease due to combustion is assumed constant. This entropy 
variation is shown in (5) to be a higher-order small quantity if 
the square of the Mach numbers of the gas flow system is 
negligibly small. Thus the gas flow system can be consid- 
ered as isentropic to the proper order of approximation. A 
concentrated combustion front in such a gas flow system does 
not separate two regions of gases having significantly differ- 
ent thermodynamic states, but is only a sharp discontinuity of 
the mean flow velocity of the gas. The mean velocity dis- 
tribution along the combustion chamber axis indicates the 
distribution of combustion, while the pressure, the density, 
and the temperature of the gas in steady-state operation are 
essentially uniform throughout the combustion chamber 
within the proper order of approximation. 

The shape of the velocity profile or the distribution of com- 
bustion in a liquid propellant rocket motor varies consider- 
ably. It is often found that most of the combustion is con- 
centrated in a narrow region. Therefore, as a rough approxi- 
mation, we consider the combustion as a sharp discontinuous 
front. In (4), this concentrated combustion front is assumed 
to be near the injector end, which means the length of the com- 
bustion chamber is much longer than what is necessary. In 
the present paper we shall consider the concentrated combus- 
tion front to be located at arbitrary axial position, and the 
problem will be formulated for a model with two concentrated 
combustion fronts, one situated near the injector end, and 
the other at arbitrary axial position. Thus we have a two- 
step velocity profile which can be reduced to different simpler 
limiting cases of special interest. 

Consider an idealized liquid propellant rocket motor whose 
injectors provide at a constant rate two uniform streams of 
propellant elements. The propellant elements in the first 
stream have a common small value of total time lag so that 
these elements that are injected into the combustion chamber 
at the same instant will burn simultaneously at a place very 
close to the injector end. The propellant elements in the sec- 
ond stream are assumed to have a uniform total time lag which 
is much larger than that of the elements in the first stream 
and will burn at a distance £ from the injector end. In 
steady-state operation the hot gas generated from the first 
concentrated combustion front moves downstream with a ve- 
locity u*, in the region 1 bounded by the two concentrated 
combustion fronts. At the second concentrated combustion 
front a new mass of burned gas is introduced into the system 
and thoroughly mixed with: the burned gas from region |. 
They move downstream as a single unit with velocity u*, in 
the region 2 bounded by the second concentrated combustion 
front and the exit of the combustion chamber. Under the 
present simplified model there is no combustion taking place 
anywhere ‘else except at the two concentrated combustion 
fronts. Thus u*; and u*, are constants in the two regions 
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DRAWING OF ROCKETS 
COMBUSTION FRONTS 


SCHEMATIC 


WITH CONCENTRATED 


indicated by the subscripts 1 and 2 as shown in Fig. 1 (a). 
Both when u*; ~ u*, and when & > 0, as shown in Fig. 1 (b) 
and 1 (¢), we have the case of a single concentrated combus- 
tion front at the injector end. When u*, — 0, as shown in 
Fig. 1 (d), we have a single concentrated combustion front 
at arbitrary position ¢ The dimensionless velocities 1 
and uw, obtained by dividing u*, and “*, through the stag- 
nation sound speed c*p as a reference quantity are assumed 
to be so small that 2 and Ww, coinciding practically with the 
square of the Mach numbers of the gas flow, are negligible 
compared to unity. 


2.2. The Time Lag and the Burning Rate 


It is explained in (4) that during the total time lag, the pro- 
pellant elements undergo a series of complicated processes 
which ultimately lead to complete combustion of these pro- 
pellant elements. Some of these processes, such as the atomi- 
zation of the fuel and the oxidizer and the proper; mixing of 
such atomized particles, are rather insensitive to the pres- 
sure and temperature oscillations of the burned gas in the 
combustion chamber. Many other processes, such as the 
vaporization of the propellant elements and the activation 
through ordinary heat transfer or other means, are rather 
sensitive to both the pressure and the temperature oscillations 
of the burned gas. The total time lag 7, is therefore composed 
of a constant or insensitive part 7; and a varying part 7 which 
is sensitive to the oscillations of the burned gas in the combus- 
tion chamber. Without any precise knowledge of these proc- 
esses taking place during the period of the time lag, we have 
to assume some form of the dependence of 7 on the pressure 
and the temperature oscillations of the burned gas. For the 
case of small oscillations of the burned gas about the steady- 
state conditions, it will be assumed that the temperature and 
the pressure oscillations are correlated and that the effect of 
the temperature oscillations can be expressed in terms of the 
pressure oscillations. The average rate of the rates of the dif- 
ferent local processes is assumed to be proportional to a con- 
stant power n of the local gas pressure acting on the propel- 
lant element. Thus, the relation defining the pressure sen- 
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sitive time lag r of the element burning at the instant ¢ is 
given in (4) as 
fi, pale"), U1 dt! = const C........ [2.2.1] 
where p [z’(t’), t’] is the gas pressure acting on the propellant 
element at the position X’(¢t’) and at the instant t’. Both the 
constant value C of the integral and the constant pressure 
index n of interaction are characteristic constants of the pro- 
pellant under the operating steady-state chamber pressure. 
The index n is assumed to be constant throughout the time 
lag period and therefore represents the average extent of in- 
teraction between the combustion processes and the pressure 
oscillations. The effect of the temperature oscillations is also 
included in the index n. Our knowledge of the kinetics of 
the individual process is not sufficient for a theoretical pre- 
diction of such an over-all parameter. This constant n for 
given propellant can be determined through experiments only. 
Now we shall see how this pressure-sensitive time lag can 
lead to varying burning rate in the simplified model as ex- 
plained in the previous section. Consider for simplicity the 
case of a single concentrated combustion front. Let m, (t) be 
the total amount of the propellant injected from the beginning 
of the operation up to the instant t, and m,(t) be the total 
amount of the propellant burned up to the same instant t. 
At this instant ¢, the propellant injected into the combustion 
chamber during the interval t —7, to ¢ has not burned. 
Therefore 


The instantaneous burning rate is hence 
m,(t) = 


. dr\ | 
= (: = (1 — [2.2.3] 


Asweare considering the intrinsic type of combustion instability 
in the idealized liquid propellant rocket as described in section 
2.1, the injection rate m,(t —7,) is a constant equal to mi. 
Differentiate equation [2.2.1] with respect to ¢ to obtain 


1-(dr/dt) and substitute 1-(dz/dt) into equation [2.2.3]. The 
burning rate is thus given as 
m,(t) = m [2.2.4] 


Since the gas pressure varies with time and position, the burn- 
ing rate is sometimes larger and sometimes smaller than m,, 
which is equal to the burning rate in steady-state operation. 
This variation of burning rate produced by the pressure oscil- 
lations in the combustion chamber is the self-exciting mecha- 
nism of producing the intrinsic combustion instability. 

In the case of low-frequency oscillation, when the length of 
the combustion chamber is much smaller than the wave length 
of the pressure oscillation, the pressure in the combustion 
chamber is almost uniform at any instant. Thus p[z(t—7), 
t—r]is approximately equal to p[t—7]. In the case of high- 
frequency oscillations, when the length of the combustion 
chamber is of the same order of magnitude as the wave length 
or several times larger than the wave length of the pressure 
oscillations, it is not obvious that p[x(t— 7), t—7] can be re- 
placed by p[.x(t), £—7] as is done in (4). However, it is shown 
(5) that the spacewise variations of the pressure in the com- 
bustion chamber has a contribution which is a higher-order 
small quantity in the stability calculation compared with the 
contribution of the timewise variation of the chamber pres- 
sure at a given location if the pressure sensitive time lag r is 
much smaller than the pressure insensitive time lag 7,. In 
the simplified model which we are considering now, 7 is much 
less than 7;, and the spacewise distribution of combustion is 
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approximated. We shall hence put p|z(t—7r), t—7r] 
t—r], and the burning rate m,(t) is given as 


my(t) = 


where both pressures are evaluated at the same location x 
where the propellant element burns. It should be noticed 
that when the total time lag varies with the gas pressure, the 
position where the propellant elements burn also varies. We 
shall neglect the oscillation of the concentrated combustion 
front about its steady-state position and evaluate p(t) and 
p(t—r) at the steady-state position of the concentrated 
combustion front. 


2.3 Small Perturbation Equations for the Gas System 


In analyzing the low-frequency oscillations, the combustion 
chamber pressure is assumed to be uniform at any instant but 
fluctuates as a whole. This is justifiable because the period 
of the low-frequency oscillation is very large compared to the 
time required for a pressure wave to travel the length of the 
combustion chamber, so that before the periodic oscillation of 
the gas has produced any appreciable changes in the gas prop- 
erties, the pressure wave has traveled many times back and 
forth and has made the flow properties nearly uniform. Thus 
the consideration of mass balance is sufficient to formu- 
ulate the low-frequency oscillation problem. In the case of 
high-frequency oscillations, the characteristic time for the 
pressure wave to travel the entire length of the combustion 
chamber is comparable to the period of the oscillations. 
Therefore the wave propagation phenomena will have to be 
considered along with the combustion process. The gas 
dynamic equations of continuity, momentum, and energy will 
then be used. From previous discussion in section 2.1 the 
energy equation will be replaced by the equation of isentropic 


change of state which is justified within at oo order of 


approximation. Thus we have 


p*[u* + u*u*z+] = . [2.3.1] 
p*p* = const 


where all quantities are dimensional and subscript ¢* or x* 
means the partial derivative with respect to the corresponding 
variable. 

Use the following scheme to make al] these quantities di- 
mensionless 


with superscript * indicating that the quantity is dimensional, 
and subscript » indicating that the quantity is evaluated at 
the stagnation condition. C™%*) is the sound speed in the stag- 
nant gas and L is the length of the combustion chamber. 
The characteristic time 6/2 is defined as the time required for 
the sound wave to travel the combustion chamber length in a 
stagnant gas, i.e., 0/2 = L/C*). This characteristic time is 
one half of the characteristic time defined in (4) under the 
approximation u? < <1. Equations [2.3.1] when expressed 
in dimensionless quantities are ge (a4 
pi + (pu), = 0 


These equations govern the unsteady flow of the gas in region 
1 and in region 2, respectively, separated from each other by 
the velocity discontinuity at the concentrated combustion 
front. 


For the study of the small oscillations in the gas system, 


steady-state flow. Thus define 
=ptp’, 


where the mean quantities p and f are practically unity in 
dimensionless form under the approximation u?< <1. The 
mean velocity u is constant either in region 1 or in region 2. 
Both % and tw are assumed to be small such that u? is much 
less than unity. Introduce these perturbations into Equa- 
tion [2.3.3] and linearize the equations with respect to small 
and we have 


andu =u+u’ 


p=pte’, 


ll 


= rp’ 


These equations are recognized as equivalent to the simple 
wave equations, and admit solutions of the type uw’ = wu’(¢) 
and p’ = p (¢) with ¢ = t—az. Substituting these into 


Equation [2.3.4], one obtains 
| 


In order to have nonzero solutions of u’y and o's, we must 
have (1 —ua)? —a? = 0. Therefore only two values oi a 


(1 — wa)u’s — = 
—au's + (1 — tia)p’s 


— anda, = — - where 1/a, is 
1+u l1—u 


the speed of the downstream moving wave, and 1/a, is that of 
the upstream moving wave as observed from the combustion 
chamber wall. Thus the general solution of Equations 
[2.3.4] is 


are possible, i.e., a, = 


p’ = p’(t — an) + p’s(t — [2.3.5] 
= — apt) + — 


Put Equations [2.3.5] into Equations [2.3.4], separate the up- 
stream and the downstream moving waves, and integrate with 
the boundary condition that in steady state both p’ and wu’ 


must vanish. We find the relations 
Hene e Equations [2.3.5] become 
p’ = u'-(t — apt) — — apt)) 
a’ u'(t — apr) + u’,(t — 


Let us investigate the stability of periodic solutions of 
exponential type 


p's pS —u', 


[2.3.6] 


u’, = cr exp [a (t — are)] 2.3.7) 


where c, and ¢, are integration constants, and a = \ + iw 
with \ = amplification coefficient and w = angular frequency 
of the wave. For simplicity, let us also write the perturba- 
tions as 


= v(x) exp (at), p’ = exp (at), and p’ = ¢(z) exp (al) 


Then the solutions for the functions v(x) and 6(z) are 


= cr exp (—arar) + 
= c, exp (—a,ar) — cs exp (—a,az) 


These solutions apply to regions 1 and 2 respectively. In 
region 1 the ratio of the two integration constants c,,/c,, can 
be determined by the boundary condition at the injector end. 
In region 2, the corresponding ratio can be determined by the 
boundary condition at the combustion chamber exit. The 
two sets of solutions in the two regions will have to be matched 
at the second concentrated combustion front as required by 
the boundary conditions at such front. This matching of 
the two sets of solutions defines completely the complex 
quantity a = } + ww. If d is positive, the disturbance will 
grow exponentially with time and therefore is unstable. 
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If \ is negative, the disturbance will die out exponentially 
with time and is stable. The boundary between the stable 
and the unstable regions as a relation between the character- 
istic constants of the gas flow system, 7, n, ¢, u etc., is ob- 
tained if we put \ = 0. The determination of such stability 
boundary is one of the major objects of the present investiga- 
tion. 


24 The Boundary Condition 


The boundary condition at a concentrated combustion 
front will be investigated first. As has been explained in 
section 2.1, the concentrated combustion front is not a dis- 
continuity of pressure, density, and temperature but is only 
a discontinuity of velocity of the flow. Hence the boundary 
condition at a concentrated front consists of two parts: 

1 The steady-state values as well as the small perturbation 
values of the gas pressure and the gas density are continuous 
at every instant across the concentrated combustion front. 
That is, p’2 = p’: and p’2=p’;. These are equivalent to the 
conditions 


= a =a + 


The last equality in Equations [2.4.1] indicates that the oscil- 
lation frequency and the amplification rate are the same on 
both sides of the concentrated combustion front. 

2 The fractional increase of the difference of mass flow 
rates across the concentrated combustion front is equal to 
the fractional increase of the burning rate at the concentrated 
combustion front. The fractional increase of the burning 
rate (m» —mi)/m; can be obtained from Equation [2.2.3]. 
But in the analysis of the simplified model, we shall use Equa- 
tion [2.2.5]. Thus neglecting higher-order small quantities, 
we have 


ty — dr pr(t) 
pr(t — 7) 


The second boundary condition at the concentrated combus- 
tion front separating region 1 and region 2 is obtained by 
equating (m, — m,)/m, to the fractional increase of the dif- 
ference of the mass flow rates. Thus replacing the small per- 
turbations by their periodic form and canceling the common 
factor exp (at), we obtain 


dt 1 = n[p'(t) — — 


— + — m)(1 — + — exp (—a7) = 0 
[2.4.2] 


In Equation [2.4.2] we have replaced 7 by = in the coefficient 
of the small perturbation 6,, neglecting the difference 7 — 7 as 
a higher-order small quantity. 

If the concentrated combustion front is located at the in- 
jector end, the upstream side of the combustion front has no 
oscillation. Hence the boundary condition at x = 0 is ob- 
tained by putting the disturbance and the mean velocity of 
the upstream flow in Equation [2.4.2] to zero. 


(1 + exp (—a7) =(Q@... [2.4.3] 


Now we come to the boundary condition at the combustion 
chamber exit where the gas enters the converging section of 
the de Laval nozzle. The reflection of a one-dimensional 
pulse at the entrance of the nozzle is essentially three-di- 
mensional and will lead to a complicated problem somewhat 
like the problem of Mach reflection from a wedge. The result 
of such analysis, even if it could be obtained in reasonably 
simple form, is not suitable as a boundary condition for the 
one-dimensional flow in the combustion chamber. To be con- 
sistent the boundary condition must be obtained from one- 
dimensional consideration, neglecting completely the two- 
dimensional effects. The physical boundary condition at the 
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nozzle entrance is that there are no discontinuities of the gas 
flow properties at this station. In other words, the unsteady 
one-dimensional motion of the gas in the combustion chamber 
and that in the nozzle must be joined continuously at this 
entrance. In (6) the unsteady flow in a de Laval nozzle with 
linear steady-state velocity in the subsonic portion is deter- 
mined for arbitrary frequency, and the boundary condition 
at the entrance to the nozzle is presented in graphical 
form as the ratio of the fractional variations of velocity and 
density at the entrance. This boundary condition, when 
y = 1.20, is reproduced in Fig. 2 with I = (»/u)/(6/p) = 
R + iS plotted against the reduced frequency 8 = wlsu»/- 
(V2 (y+1)—u) where J.» is the length of the subsonic 
portion of the nozzle as a fraction of the combustion chamber 
length. 

If the length of the subsonic portion is very small, i.e., 
leu. & 0, the values of 8 corresponding to the fundamental 
or the first few modes of oscillation are very close to zero. 
Hence for very short nozzle, one has 


10, i) = (2.4.4) 
6/p/ B=0 2 


if the one-dimensional result holds good in such limiting case. 
This boundary condition has been shown in (6) and (7) to 
correspond to constant Mach number of the gas flow at the 
entrance of the nozzle. 


2.5 Final Formulation 


Having established the boundary conditions for the solu- 
tions in different flow regions, we can proceed to formulate 
the equation for the determination of the complex quantity 
a = + w fora given system. 

For region 1, the boundary condition at the injector end 
as given in Equation [2.4.3] can be used to determine the 
ratio of c,, and c,,. Call this ratio — A. 

1 — — yn) — tiyn exp (—a7) 


1+ — yn) + tyn exp (—a7) 


The solution in region 1 can hence be written as 


= cs, exp (—aa,x)[1 — A exp { — a(a, — a,)x}] 


= —Cs, exp (—aa,x)[1 + A exp { — a(a, — a,)x}] 


6 
I( Ba): =R+iS eq 
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FIG. 2 RECIPROCAL OF NOZZLE IMPEDANCE FOR OSCILLATORY 
FLOW OF FREQUENCY w = 78 
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In a similar manner the boundary condition at the combus- 
tion exit X = 1 can be used to determine the ratio 


C,,/Cs, = —B exp [a(a, — a,)]........ [2.5.3] 


When the nozzle is very short the boundary condition given 
as Equation [2.4.4] is used, and we have 

1+ (y — 
1 — (y — 1)a/2 


When the nozzle is long, B is a fraction of B and u with1(8, u) 
given in Fig. 2. 


1 + 1(6, 


BG, = ae 


In general, we can write the solutions in region 2 as 


= c,, exp (—aa,x) [1 — B exp { a(a, — a,)(1 — x)}] 
8(z) = — ¢,, exp (—aa,r)[1 + B exp {a(a, — a,)(1 — x)}] 
.. [2.5.6] 


Now the two sets of solutions [2.5.2] and [2.5.6] are to be 
matched at the second concentrated combustion front z = é. 
The oscillation of the second concentrated combustion front 
about the mean steady state position tis neglected. Equating 
6,(é) and 6,(¢) as given in Equations [2.5.2] and [2.5.6], we 
have 


Co 1 + A exp [—a(a, — a,)é] Ps 1 + A exp (—2 at) 
B exp[a(ar — a,)(1 — Bexp [2 a(1 — 


where a, — a, = 1/(1 + uw) + 1/(1 — &%) = 2/(1 — w) =2 
has been substituted. 

By introducing Equation [2.5.2], [2.5.6], and [2.5.7] into 
Equation [2.4.2] and dividing the resulting expression with 
6;, we have 


1 — Bexp [2 a(1 — £)] 4 1 — A exp [—2 aé] - 
1 + Bexp [2 a(1 — §)] 1 + A exp [—2 aé] Zz 


(tig — %)[((1 — yn) + yn exp (—a7)]......... [2.5.8] 


This is the final form of the equation for determining the com- 
plex quantity a = \ + iw of a combustion system with two 
steps of concentrated combustion, one at the injector end 
and the other at ¢. For the determination of the stability 
boundary, we put \ = 0 or a = ww, then separate the real and 
the imaginary parts in Equation [2.5.8] to get two real equa- 
tions from which we can eliminate w. The eliminant is the 
equation defining the stability boundary. We see that in 
Equation [2.5.8] both A and B are, in general, complex quan- 
tities. The separation of the real and the imaginary parts of 
Equation [2.5.8] is quite laborious and the resulting real 
equations cannot be handled easily. 

For the case of a single concentrated combustion front 
at arbitrary position ¢ (Fig. 1 (d)), we put z, = 0 in Equation 
[2.5.8]. Thus 


1 — Bexp [2 a(1 — &)] = 1 — exp (—2 aé) 3 
1+ Bexp [2e(1— 1+ exp(—2aé) 


ti,[((1 — yn) + yn exp (—a7)]............ [2.5.9] 


When the nozzle is short, B(0, 2) is 
a real constant given by Equation [2.5.4]. Thus a is the 
only complex quantity in Equation [2.5.9]. The solution of 
the problem is therefore greatly simplified algebraically. 
When the nozzle is long the value of B is given by Equation 
|2.5.5] where J is a complex function given in Fig. 2. Then 
the determination of the stability boundary will have to be 
done graphically. If ¢ = 0, Equation [2.5.9] becomes 
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where A becomes unity. 


2 
2 cos? w 
yn sin wr = tan w 


(1 — yn)tiz + yntie exp (—a7) = 1 — B exp (2 a) 


1 + B exp (2 a) 


This case corresponds to a single concentrated combustion 
front at the injector end. This Equation [2.5.10] can also be 
obtained from Equation [2.5.8] by taking the limit of either 
t = 0 or & = & (Figs. 1 (b) and 1 (c)). When yn = 1, 
Equation [2.5.10] reduces to a form identical with Equation 


[12.15] in (4). 


Solution with Short Nozzle 

3.1 Single Concentrated Combustion Front at Injector 

End 


Rewrite Equation [2.5.10] and drop the subscript 2. 


1 — B exp (2 a) 


1 + B exp (2 a) I 


ynii exp (—az) — (yn — 1)a = 


where 


Putting a = ww for the neutral oscillations and separating 

the real and the imaginary parts of Equation [3.1.1] we have 
1 — B? 

1 + B? + 2B cos2w 


ynii cos we = —(1 — + 
[8.1.2] 


” 2B sin 2 w “= 3 
1+ B2+ 2B cos2w 
with the approximation % << 1 we can write 


Sin wi 


(y — 


1+ B? 1+ [(y — 1)a#/2]* 
Thus Equation [3.1.2] beeome 


yn COS wF == —(1 — yn) — 


When vn is of the order of unity, tan w is of the order of 4; 
therefore we can neglect tan? w as compared to unity. Thus 
the first relation in Equation [3.1.4] immediately gives 


yn 
= 
ind 
1\2?7 
Sin = 1- (1-44 [3.1.6] 
To 2 yn 
Thus 
ij 1\2]/2 
tanw = 1 — (: - 2+") (3.1.7 


can we have real solutions of w7 and w for neutral oscillations. 
Therefore the index n must be bigger than the minimum 
value (y + 1)/47 if the combustion system is to have unstable 
pressure oscillations. As the values of y for most of the com- 
bustion products of the common rocket propellants are usually 
close to 1.2 or 1.3, the minimum value of n is about 0.45 and 
is the same for all modes of the high-frequency oscillations 
when the constant Mach number boundary condition is used. 
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1.1] 


From Equation [3.1.7] we obtain the frequency of the neu- 
tral oscillations as 


1\2]'/2 
2 yn 


where k = 1, 2,3,..... The fundamental mode of the high- | 
frequency oscillation is obtained when k = 1; solution with - 
k = 0 is identified to be the low-frequency solution in (4) and 
is therefore discarded in the present investigation. Equa- 
tion [3.1.9] shows that the frequencies of the neutral oscilla- 
tions are close to the integral multiples of x which in dimen- 
sional form are the natural organ-pipe frequencies of the gas 
system. The fractional deviations of the neutral frequencies 
from the corresponding natural organ-pipe frequencies are of 
the order of magnitude of the Mach number of the gas flow. 

Now we shall demonstrate on which side of the stability 
boundary is the oscillation unstable. Taking B = 1 and 
differentiating Equation [3.1.1] with respect to 7, we have on 
the stability boundary 


da dy + 
da 
ynt Cos? wi 


or 


dw 


COS? wr 


dd Sin w7 


[3.1.10] 
dw COS? w — COS wi 

When n < (y + 1)/2y7, then cos wr < 0, thus the sign of 
d\ ‘dw is the same as the sign of —sin #7 or the sign of —tan 
w. When dd/dw is positive on the stability boundary, the 
region which is reached by increasing w from the boundary is 
the unstable region. (From the usual argument of continu- 
ity, this statement would be expected to hold good for n > 
(y + 1)/2y without making an effort for detailed proof.) 
Therefore, the unstable range of the frequencies is 


1\2]'/2 
ty ana (1 - <w< 
2 yn 
1\2 1/9 
ke + “ (1 “s (3.1.11) 
2 yn 


It is therefore concluded that all the frequencies of the un- 
stable oscillations are very close to the natural organ-pipe 
frequencies. 

The corresponding ranges of the values of 7 for such un- 
stable oscillations are given as 


1 
ka + (: - 
2 yn 
= — y+1 1) 
(2h + + sin ( 


ke mat (1 
where k = 1, 2, 3,... indicates the successive modes of the 
high-frequency oscillations, and h = 0, 1, 2... indicates the 
successive higher ranges of the values of 7 for unstable oscil- 
lations of a given mode. The unstable range is shaded as 
given in Fig. 3 (a) for %@ = 0.213. We see that the unstable 
ranges of the values of 7 for a given set of values of h and k 
increases with increasing n. 

When n equals the minimum value (y + 1)/47, Equations 
(3.1.11] and [3.1.12] give 


wr = (2h + el 
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The oscillation is neutral and the range of the unstable values 
of 7 vanishes at this minimum value of n. 


3.2 Single Concentrated Combustion Front at Arbitrary 
_ Axial Position 


- Rearrange Equation [2.5.9] and drop the subscript 2. a) 
exp (—az) = —(1 — yn)a + 
— B exp (2 
[1 + exp (2 at)][1 + B exp {2 a(1 — 


. [3.2.1] 


Separate the real and the imaginary parts of Equation [3.2.1] 
for the neutral oscillation with a = iw la 


1 + B? + 2B cos [(1 — £)2 a] 


2B sin (2 — £)w + (1 + B?) sin tw 
cos gw{1 + B? + 2B cos [(1 — £)2]} | 


ynt cos wr = — (1 — yn)ai + 


[3.2.2] 


sin wi = 


Using Equations [3.1.3] we have : 


ynt cos wt = — yn)i — 
7 2 cos? (1 — £)w 
[3.2.3] 


> 


=, 


sin w 


cos cos (1 — £)w 

By squaring and adding the two equations in [3.2.3], one ob- 
tains an equation from which we can solve for the frequencies 
of the neutral oscillations. For qualitative discussion of the 
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(1 — yn)(y 


me 


results, we write this equation as a quadratic of 2 cos? (1 — é) 
w and solved for 2 cos? (1 — ¢) w. Thus Aleset 


9 2(1 — 8) - 1 
cos? 2 yn — 1 wi? 
(1 — yn)(y — {| sre 1 


cos? fw 2 yn — 1 


2yn—-—1 


2 yn — 1 
The minus sign before the square root is dropped because 
cos? (1 — £) wisnon-negative. Introducing the following two 
inequalities 

sin? w 1 
cos? gw 2 yn 


2 > 2cos? (1 — t)w > 0, and 


into Equation [3.2.4] we obtain the following mocenaty con- 
dition for any possible real solutions of w 


(1 — ya)(y — 1) = 4 


2yn-1 2yn-—1 


or 


Y 
which checks with the minimum value of n as given in Equa- 
tion [3.1.8]. 

From Equation [3.2.4] we have another necessary condi- 
tion for real solutions of w if only the inequality 2 > 2 cos? 
(1 — £) w is introduced into Equation [3.2.4] 


sin? w < +n- 1] - 


As n is of the order of unity, we must have 
ae 


[3.2.6] 


cos? . 


sinw = + #:0(1) 
or 


Since the unstable regions of w for the case of ¢ =0 is defined in 
Equation [3.1.11], we conclude that the frequencies of the 
neutral and the small unstable oscillations for all possible 
values of ¢ are always close to the natural organ-pipe frequen- 
cies kx. 

The right-hand side of Equation [3.2.6] vanishes when 

= (y + 1)/4 y or when cos t»# = 0. When either of the 
two conditions is satisfied, sin w must be zero. The results 
with n = (y + 1)/4 y are the same as the result obtained 
from Equation [3.2.5]. The case cos & = 0 leads to an im- 
portant restriction on the regions of ¢ that admit real solu- 
tions of w. When both sin w and cos gw are zero, cos (1 —é)w 
must vanish; then Equations [3.2.3] cannot admit any real 
solutions for w7. 

From Equation [3.2.7] we know w is approximately kz; 
therefore the zeros of cos & are approximately ¢ = 1/2k, 
3/2k... solong ast <1. These positions correspond to the 
nodes of the pressure oscillations of the kth mode. Around 
each of these nodes there is a range of positions ¢ of the com- 
bustion front which is always stable. The critical positions 
=, which define such stable range around each node are those 
where the oscillations are neutral with w = km and zero un- 
stable ranges of 7. Equation [3.2.4] then gives ¢, as 


1 
1 =1— — cos"! 


— 1) >0. [3.2.8 
2yn— 1 )> 


For the fundamental mode, k = 1, there are two values of 
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é, lying between 0 and 1 symmetric with respect to 1/3. For 
the kth mode, we have 2k values of ¢, between 0 and 1, defining 
k stable regions about the & nodes of the kth mode of oscill:- 
tions. The results are plotted in Fig. 4 with n vs &, for k = 
1,2, and 3, and y = 1.20. It is observed that ifn = (y + 1 
4y, the critical values of é, are 0, 2/2k, 4/2k, etc., and the en- 
tire region of ¢ from 0 to 1 are stable except at those critical 
positions where the oscillations are neutral; and that if 
n <(y + 1)/4 ¥, no real values of ¢, exist. This result agrees 
with Equations [3.1.8] and [3.2.5], and in addition reveals 
the fact that if the combustion is concentrated at the injector 
end ¢ = 0, the combustion system is most liable to have un- 
stable high-frequency oscillations. The positions 0, 2/2/, 
4/2k, ete., are approximately the antinodes of the kth mode of 
pressure oscillations where the amplitude of the pressure os- 
cillations is the largest. 

Equation [3.2.4] can be used to solve the frequency of the 
neutral oscillation for arbitrary values of ¢& But with Equa- 
tion [3.2.7], a rapidly converging iteration procedure can be 
used to solve the critical values of 7 and w from equations 
[3.2.3] and [8.2.4]. By replacing (1 — £) w® and gw 
with (1 — &) kw and &kz, we can calculate sin w from Equation 
[3.2.4] and get wow = ke + 1%. Using w in the place 
of w, we can determine w® = kx + I1-% and repeat the 
process until the necessary accuracy isacquired. Then Equa- 
tion [3.2.3] gives the values of w7. For the cases calculated, 
two iterations are sufficient. Calculation is carried out for 
k = 1, 2, and 3 and h = 0, 1, and 2 with n = 1/y = 0.833, 
and u = 0.20. The results are plotted as given in Fig. 5. 
These curves are not symmetric. The shaded regions for a 
given mode of oscillation represent the values of 7 and é that 
make this mode unstable. We see that there are only a few 
spots where all the three modes are stable. 

From this calculation, it is clear that the position of the 
concentrated combustion front along the combustion chamber 
axis has considerable importance. If the combustion is 
mostly concentrated in a narrow region whose width is only a 
small fraction of the distance between two adjacent positions 
é, of a given mode of oscillation, the stability behavior of this 
mode can be analyzed with the simplified model of concen- 
trated combustion. If the width of the combustion zone is 
larger than the distance between two adjacent é, of a given 
mode, it is hardly possible that the simplified model could be 
used satisfactorily. As the distance between two adjacent 
positions ¢, decreases rather fast when k increases, the approxi- 
mate simplified model becomes less satisfactory for the higher 
modes of oscillation. Consequently, if most of the combustion 
in a liquid propellant rocket motor is concentrated in a narrow 
region, for example, one tenth of the length of the combustion 
chamber, the stability behavior of the fundamental and the 
second mode of the high-frequency oscillations can be deter- 
mined by using the simplified model of concentrated combus- 
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FIG. 5 STABILITY DIAGRAM, 7 VS. &, FOR PROPELLANTS HAVING 
n= 1/7 = 0.833 
tion. But the stability behavior of the higher modes as ana- 


lyzed by using the simplified model should not be considered 
too seriously. 


3.3 Two-Step Concentrated Combustion with n = 1/4 


For the study of the nature of the solution of the two-step 
concentrated combustion with one at the injector end. we 
shall consider the representative case of yn = 1 for the pur- 
pose of simplicity. Equation [2.5.8] becomes, 


1 Bexp Ba(l- _ 
1 + B exp [2 a(1 — §)] 

{1 + exp (—az)] — [1 — exp (—az)] exp (-2 ag) 
{1 + exp (—az)] + [1 — % exp (—a7z)] exp (—2 at) 


[3.3.1] 

This equation can be rearranged in the form 
2 exp (—2 at)[1 — B exp (2 a)] — 2 exp (—2 at) [1 + B exp 
(2 a)](%2 + %&) exp (—az) = [1 + B exp [2 a(1 — 2] + 
+ B exp [2 a(1 — SJ}{1 — exp(—2 at)}a exp (—a7)] 
(tz — exp (—a7) 
The coefficient of exp (—2 a7), which is of the order of (@ — 
u,) %, will be considered as negligibly small compared with the 

coefficient of exp (— az). Thus 

exp (a7) = + [1 + B exp (2a)] + — [exp 
(2 at) + B exp {2a(1 — £)}]] [1 — B exp (2a)]~?. . [3.3.2] 


Let a = w for the neutral oscillations and separate the real 
and the imaginary parts of Equation [3.3.2] with the approxi- 
mation of Equations [3.1.3] 


COS WT = —(y 1) t2[!/2( ide + + 
— cos 2 tw] /[1 — cos 2 { 
sin wt = + sin 20 + 


— {sin 2 + sin 2(1 — t)w}]/[1 — cos 2w] | 
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The equation for determining the frequencies of neutral oscil- 
lations of the system is obtained as 


4sintw = + — + 4 sin? weos?w) + 
— — cos? 2 Ew + 4 sin? w cos? (1 — 2 
— 1,2) [(y — cos 2 Ew + 4 sin? w cos w cos (1 — 2 
. [3.3.4] 


It is not easy to solve for w from Equation [3.3.4] directly. 
But a method of successive approximation, analogous to the 
method used in section 3.2 can be used. From Equation 
[3.3.4] we see that if wu. — % issmall compared to u , + %, we 
have as a zero-th approximation w = kr + wu. which is the 
solution given in (4). Introducing »® into the right-hand 
side of Equation [3.3.4], we can calculate w from 
lg — ty 

9 


ii: 3 
4 sint > — ( ) sin? (1 — 2 — 


= {1 — cos cos (1 — 2 4,2. [3.3.5] 


The solution of the frequency of the neutral oscillations can 
thus be written as 


w + Itiz..... .... [8.3.6] 


where / is a function of ¢ and d/t%. From Equation [3.3.5] 
it is obvious that the coefficient / will be less than unity. At 
the nodes of the Ath mode of pressure oscillations, 1.e., 
&¢ = 1/2k, 3/2k,... ete., Equation [3 3.4] gives the minimum 
value of sin w as +(% provided is not toosmall. Con- 
sequently, the frequencies of the neutral oscillations of the 
system are not significantly different from the natural organ- 
pipe frequencies when a small] fraction of the concentrated 
combustion is shifted from the injector end to arbitrary loca- 
tion along the axis 

Having determined w, we can obtain w7 from Equation 
[3.3.3.]. The result of such calculation is given in Fig. 6. It 
is noticed that when the second concentrated combustion 
front is located at the nodes of the kth mode of oscillations, 
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1/2k, 3/2k, ete., cos wt — (y — 1) 


U2t 
2 


or = (2h + 1)r This is identical with 


2 

the result for the case of a single concentrated combustion 
located at the injector end. However, owing to the slight 
decrease of the value of sin w, the ranges of the values of 7 for 
producing unstable oscillations when the second concentrated 
combustion front is located at the nodes are slightly decreased 
as compared with that of the combustion system having ¢ 
single concentrated combustion front at the injector end. 

From Fig. 6 we see that the shift of a small fraction of the 
concentrated combustion from the injector end to any arbi- 
trary locations along the axis is slightly stabilizing. The sta- 
bilizing effect is most significant when the small fraction is 
shifted to the node. Even in this case this shift does not give 
rise to unduly large effect on the stability behavior of the sys- 
tem. It may be inferred that the distribution of a small frac- 
tion of combustion from a given concentrated combustion 
front will not change the qualitative picture of the stability 
behavior of the system as analyzed with the single concen- 


Solutions with Long Nozzle 


4.1 
Exit 


The Boundary Condition at Combustion Chamber 


The proper boundary condition at the combustion cham- 
ber exit is, as mentioned in a previous section, the continuity 
of the disturbances in the combustion chamber and the dis- 
turbances in the de Laval nozzle. When the steady-state 
velocity distribution in the nozzle is linear from the entrance 
to the sonic throat, the ratio J of the fractional variation of 
the velocity to the density disturbance at the nozzle entrance 
is obtained in (6). The real and the imaginary parts of / 
when y = 1.20 are plotted against the reduced frequency 
parameter 8 for different values of the reduced steady-state 
velocity parameter Z = ((y+1)/2) u? (Fig. 2). The reduced 
frequency £ is defined as the angular frequency of the oscil- 
lation divided by the steady-state velocity gradient in the 
subsonic part of the nozzle, i.e., 8 = w-lsu, +1) —w). 
Here /xur, is the length of the subsonic part of the nozzle as a 
fraction of the combustion chamber length. 

For the fundamental mode of high-frequency oscillation, 
The value of 8 corresponding to the 


trated combustion front. we know w is about z. 
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fundamental mode is approximately 2/(y+1) — u]. 
Since u is relatively small, the value of 8 for the fundamental 
mode depends primarily on the value of lsu. In practical 
eases sup is often close to 1/3. If lsu, of the nozzle is very 
small, the value of 8 corresponding to fundamental mode is 


1 — Bexp {2 a(1 — 
1+ Bexp !2 a(1 — 


— tanh aé |- 
u 


(1 — yn) + yn exp (— az) 


iw for the neutral oscillations and B = (1 + J) 


ait very small. Therefore, if the one-dimensional picture could hieally in Fie. 2 f 
. = ‘ S “4 y 4 
lin be extrapolated to such limiting case, the boundary condition ( 
= as or w). stability 
: of constant Mach number would apply satisfactorily to the y -20 as a function of 2(01 u) and B (or w) he stability 
me: case of a nozzle with a very short subsonic part. On the boundary is determined numerically for two particular cases 
a other hand, if lsu, is reasonably long, the stability behavior - = 0.05 and 0.10 
of the fundamental mode can be considerably different from 7 
ci that determined by the constant Mach number boundary corresponding to 
condition. 
A i = 0.213 and 0.301 
We shall consider a practicalexample with [V2/(y+1) 
ced -u| = wsuch that o = 8x, Thus the fundamental mode of For a given value of ¢ and a series of values of w, the left-hand _ 
cil acoustical oscillation corresponds to 6 = 1 and the kth mode side of Equation [4.1.1] is calculated. Call this quantity 
the corresponds to 8 = k. It should be noted that for different 
u) values of u, the condition w» = requires that nozzles of 1 [1 Bexp — 
different lengths are used with the combustion chamber. + Bexp — 
4.2. Solution for the Case of a Single Concentrated Com- 
on bustion Front Then the values of n corresponding to the series of values 
the }-quation [2.5.9.] can be rewritten as of w are given by 
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and the corresponding critical values of 7 by 


[4.1.4] 


The value of sin—! (— }'/yn) is taken in the quadrant consistent 


with cos wr = 1— ((1—X)/yn). 
plotted as shown in Figs. 7, 8, 9, 10, 11. 


The calculated results are 


For the purpose of 


comparison, these curves are also calculated and plotted when 
the boundary condition J (0, u) = (v—1)/2 is used. 


It is noticed in Figs. 7 (a) and 7 (b) 


that curves of n vs. 


8 with the two different boundary conditions are of similar 


shape and vary with the same trend. 


The curve with the 


boundary condition /(8, u) = [(v/u)/(6/p)] for long 
nozzle is shifted toward larger values of n and smaller values 
of 8 as compared to the corresponding curve with the constant 
Mach number boundary condition for very short nozzle. As 
a result, there are many apparent but important results with 
respect to the effect of changing the ratio of the length of the 
subsonic portion of the nozzle to the combustion chamber 


length. 


1 The minimum value of n compatible with any unstable 
pressure oscillations determined by the long nozzle boundary 


condition increases rapidly with increasing 8 or w. 


This is 


shown in Fig. 9. When the constant Mach number or short 
nozzle boundary condition is used, the minimum value of 
n for al] the different modes of oscillation is the same constant 
(y+1) /4 7, slightly less than !/2, which is the minimum value 
of n for the unstable low-frequency oscillation corresponding 


tok = 0. 


2 When ¢ increases, the rate of increase of the minimum 
value of n determined by the long nozzle boundary condition 
is larger than the corresponding rate determined by the short 


nozzle boundary condition. 


Both results are plotted in Fig. 


4 as the dotted and the solid curves, respectively. 

3 For a given value of n, the unstable region of 7 and ¢ as 
determined by the long nozzle boundary condition are smal- 
ler than the corresponding unstable region as determined by 


the short nozzle boundary condition. 


Owing to the rapid 


increase of the minimum value of n for increasing w, there is an 
upper limit of the frequency of possible unstable oscillation 


as determined from Fig. 9 for given value of n. 


different from the result given in section 
there is no such upper limit. 
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= 1/y = 0.833 and uw = 0.213, using the long nozzle bound- 
ary condition, shows that the fundamental mode will be 
unstable when the values of = and é are in certain ranges while 
all the other modes are always stable (Fig. 10). For the case 
of n = 1, both the fundamental and the second mode may be 
unstable but not the higher modes (Fig. 11). 

Consequently, in every respect, the combustion system 
with long nozzle is more stable than the one with very short 
nozzle. 

It should be noticed that the rate of increase of the mini- 
mum value of n compatible with unstable kth mode of oscil- 
lation increases rather fast as k increases. The actual mag- 
nitude of n of the common liquid propellants is not known. 
However, being an exponential index, n is not likely to be 
much larger than unity. Therefore, for practical cases the 
higher modes of pressure oscillations are expected to be stable. 
Furthermore, even if these higher modes become unstable, 
they are not likely to build up rough combustions owing to 
the more effective viscous damping of the high-frequency com- 
ponents. It is the unstable fundamental and the first few 
high-frequency modes that result in rough combustion and 
heuce are of practical interest. 

A comparison of the results in Figs. 7 (a) and 7 (b) and in 
Fig. 9 for the two cases with different wu shows that by increas- 
ing « from 0.213 to 0.301 while the nozzle is so modified as to 
maintain the relation w = 87, the minimum value of n com- 
patible with unstable oscillations are slightly decreased; 
that is, the system becomes a little more unstable. The 
change is, however, rather small as produced by 40 per cent 
increase of the flow Mach number. Whether the effect of 
increasing u is destabilizing or stabilizing is as yet uncertain, 
because this small difference in n can be easily compensated 
by the terms of the order of u? which have been neglected in 
the present analysis. It seems, however, that the stability 
behavior is rather insensitive to the change of flow velocity 
in the combustion chamber when lsu, ‘[V/2/(y+1) — is 
kept unchanged. 

As has been mentioned in section 3, if the combustion is 
reasonably concentrated, the simplified model of concentrated 
combustion front can be used to analyze the stability behav- 
ior of the fundamental and the first few high-frequency modes 
of oscillation, but cannot be used to analyze the stability be- 
havior of the higher modes with the constant Mach number 
boundary condition at = 1. Now with the boundary con- 
dition for long nozzle at x = 1, all that we are interested in is 
the stability behavior of the fundamental and the first few 
high-frequency modes, because all the higher modes are ex- 
pected to be stable. Consequently, the simplified model of 
concentrated combustion is a convenient model for analyzing 
the stability of the high-frequency oscillation if the combus- 
tion zone lies in a region which is small compared to the com- 
bustion chamber length. 

From the calculated results with y = 1.20, we see that if 
the combustion zone lies either between & = 0.32 and ¢ = 
0.45, or between ¢ = 0.57 and é = 0.78, all the high-frequency 
modes are stable for any values of 7 when n is not larger than 
1. Quantitative results as obtained in this calculation will 
be somewhat different when the steady-state velocity profile 
in the subsonic portion of the nozzle is not linear; but the 
qualitative results as shown in the present calculation is 
expected to remain unchanged. 


Summary of Conclusions 
From the previous results, we can draw in general the fol- 
lowing qualitative conclusions: 

| The frequencies of the unstable pressure oscillations in 
the combustion chamber excited by the interaction between 
the combustion process and the pressure oscillations are close 
to the natural frequencies of the combustion system. 

2 For a given configuration of the concentrated combus- 
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tion unstable pressure oscillations are possible when the 
pressure index n of interaction is sufficiently large and when 
the value of the pressure sensitive time lag 7 of the system 
lies in certain ranges which are functions of n. These ranges 
of = for unstable pressure oscillations of a given mode in- 
creases if n increases from its minimum value. 

3 There is an absolute minimum value of n which is 
(y + )/4y. Ifa combustion system has a value of n less than 
(y+1)/4y, all the high-frequency oscillations are stable no 
matter what the value of 7 is, where the combustion is con- 
centrated, and what boundary condition at the exit of the 
combustion chamber is used. 

4 The minimum value of n compatible with unstable pres- 
sure oscillations of a given mode increases when the concen- 
trated combustion front is further away from the nearest 
antinode of the pressure oscillations of that mode. This 
minimum value of n for a given mode becomes “‘infinite”’ 
the concentrated combustion approaches any of the nodes 
of that mode of pressure oscillation. The injector end 
an antinode of all the modes of oscillations and the configu- 
ration with concentrated combustion front at the injector 
end is most liable to become unstable. 

5 When the concentrated combustion front is at the in- 
jector end, the minimum value of n for the kth mode to be 
unstable increases rapidly with increasing k except when the 
nozzle is very short and the constant Mach number boundary 
condition is to be used at the axial exit. The higher modes of 
oscillations are relatively more stable than the lower modes 
in practical cases. 

6 Rockets with very short nozzles are more unstable than 
similar rockets with long nozzle. A rocket with short nozzle 
has lower minimum value of n and larger unstable ranges of 
7. The qualitative stability behavior of the two rockets are, 
however, similar and the results as determined by short 
nozzle boundary condition are a helpful qualitative guide in 
the numerical determination of the stability boundary with 
long nozzle boundary condition. 

7 The stability behavior of a system with concentrated 
combustion is rather insensitive to the change of the flow 
Mach number of the gas in the combustion chamber if the 
velocity gradient in the subsonic portion of the nozzle is 
not changed. 

8 If the combustion is mostly concentrated in a region 
whose width is only a small fraction of the combustion cham- 
ber length, the stability behavior of the fundamental and the 
next few higher modes of oscillations can be satisfactorily ana- 
lyzed by using the simplified model of concentrated combus- 
tion front. The higher modes of oscillations are expected to 
be stable when the long nozzle boundary condition is used at 
the combustion chamber exit. If the combustion is distrib- 
uted so that the combustion zone covers considerable por- 
tion of both the stable region and the unstable region of tof | 
the fundamental mode, there is no obvious position of oa 
concentrated combustion front in the simplified model which 
can be satisfactorily used for the analysis of the stability be- | 
havior of the actual system. 

Since the analysis is made on a one-dimensional basis in — 
the axial direction, all results apply to longitudinal oscilla- — 
tions only. 
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The aerodynamic forces on a cylinder which is small 
_ compared to the mean free path are calculated for the case 
* of an incident gas flow which is in non-Maxwellian equilib- 
rium. This nonequilibrium condition is associated with 
_ the presence of viscous stress and heat flux terms. It is 
shown that the shear stress and normal heat flux give rise 
to lift forces on the cylinder, while the tangential heat flux 
and the deviation of the normal stresses from the hydro- 
_ static pressure give rise to additional drag forces. A modi- 
fied version of the interaction of gas molecules and sur- 
faces is introduced which involves three reflection coef- 
ficients governing the average transfer of energy and nor- 
mal and momentum. 


N 
omenclature 
Wire: 
= element of surface area 
= drag coefficient, Equation [3.4] 


= lift coefficient, Equation [3.5] 


y’ 
= partial drag and lift coefficients proportional 
to term indicated in subscript, Equations 
[3.4], [3.5] 
dE;, dE,, dE,, = energy fluxes incident, reflected, and diffusely 
reflected from dA, Equation [2.7] 
Fp = drag force per unit length of cylinder b frol 
Fy = lift force per unit length of cylinder (=~ 
dF = force on dA 7 
— 
dF; = component of dF due to incident molecules 
— 
dF, = component of dF due to reflected molecules 
f = molecular velocity distribution function, Equa- 
tion [2.5] 
ied = Maxwellian distribution, Equation [2.4] 
Io, 11 = modified Bessel functions of zero and first order 
k = thermal conductivity 
M = Mach number 
m = molecular mass 
Nw = number of molecules incident on dA per unit 
time per unit surface area 
p = pressure 
Di, Dir Peo = pressure on dA due to incident, reflected, and 


diffusely reflected molecules, respectively, 
Equation [2.10] 
normal stress deviation terms, Equations [2. 5}, 
[2.6] 
heat flux terms, Equations [2.5], [2.6] 7 
cylinder radius ye 
gas constant 
molecular speed ratio, Equation [3.3] 
Su = molecular speed ratio referred to cylinder tem- 
peratures, Equation [3.3] (s. assumed equal 
to s in computations) 


Prz, Puy, Pzz 


r = 
R = 


s = 


te = gas temperature 

T; = temperature of incident molecules, Equation 
[2.13] 

= cylinder temperature 
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pete Aerodynamic Forces on a Cylinder for the Free 
Molecule Flow of a Nonuniform Gas 


S. BELL! and S. A. SCHAAPE 


molecular velocity component in 2, y, 2-diree- 


u,v, w == 
tion, Fig. 1 

U,V,W = gas velocity components 

u’, v', w’ = molecular velocity components in z’, y’, 2'- 
directions, Fig. 1 

x; ¥, z = co-ordinates referred to flow direction, Fig. | 

59,2 = co-ordinates referred to surface element, Fig. | 

a = thermal accommodation coefficient, Equation 
[2.7] 

Y a = ratio of specific heats 

= viscosity coefficient 

0 - = angle between surface normal and flow diree- 
tion, Fig. 1 

p = gas density 


a, 0" = reflection coefficients for tangential and normal 
momentum transfer, Equations [2.9], [2.10] 

shear stress due to incident, reflected, and « lif- 
fusely reflected molecules, respectively, 
Equation [2.9] 


Try, Trz Tye = Shear stress terms, [Equations [2.5], [2.6] 


1.0 Introduction 

r — forces exerted upon a body in free molecule flow, i.e., 

on a body whose characteristic dimension is small com- 
pared to the molecular mean free path, have been determined 
theoretically by Tsien (1),* Stalder and Zurick (2), and others 
(10 to 12). Experiments by Stalder, Goodwin, and Creager 
(3), and by Kane and Estermann (4) have provided general 
confirmation of the theoretical developments. 

These investigations have been confined to the case in 
which the gas flow past the body is in Maxwellian equilibrium. 
This restriction is satisfactory for free flight considerations 
but not for what may, perhaps, prove to be one of the most 
important applications of free molecule flow analyses. This 
is wind tunne] measurement of the velocity, and perhaps other 
flow variables, of a gas stream by means of small “free mole- 
cule flow” probes. The possibilities are especially promising 
in connection with hypersonic or supersonic wind tunnels 
operating at low.test section densities (13). These probes 
might be of either the heat transfer or aerodynamic force 
type, or perhaps both. It is clear, however, that for probing 
boundary layers or shock waves or other characteristic parts 
of a flow field, the instrument will often be exposed to a gas 
stream which is not in Maxwellian equilibrium. The oc- 
currence of a large shear stress in a boundary layer, for 
example, produces a departure of the molecular velocity 
distribution from the Maxwellian. As will be seen below, 
this nonuniformity gives rise to additional force terms. 

Approximate analyses of the forces on bodies in such non- 
Maxwellian free molecule flows have been made by Einstein 
(5) and Epstein (6) in connection with the radiometer prob- 
lem. In their discussions the nonuniformity was confined to 
the case of a local heat flux, and in any case the results were 
only of order-of-magnitude validity. 

In the presént report a,more detailed analysis of the free 
molecular flow of a nonuniform gas is made. The results are 


* Numbers in parentheses refer to the References on page 317. 
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genera] enough to include arbitrary (small) stress and heat 
flux. The body geometry has been confined to that of a 
right circular cylinder oriented perpendicular to the gas flow 
yelocity, and it is assumed the body’s thermal conductivity 
is so large that it is at uniform temperature. Only the 
aerodynamic force characteristics are determined. The 
extensions to the determination of heat transfer characteris- 
ties, to other body geometries and to the case of very low 
internal conductivity, will be presented in subsequent reports. 
The method, however, will be apparent from the present 


paper. 
2.0 Method of Analysis Le 


21 The Incident Molecules 


The general method followed is a direct extension of the 
method for the Maxwellian case used by Tsien (1) and Stalder 
and Zurick (2), and was suggested by Goodwin. The force 


per unit area dF/dA on a surface element dA is broken up 
into a part dF;/dA due to the incident molecules, and a 


part dF,/dA due to the reflected or re-emitted molecules. 
Total force characteristics are obtained by integration over the 
surface of the body. 


The force dF;/dA is analyzed in terms of the component 
p; normal to the surface (pressure) and the component 7; 


tangential to the surface (shear stress). See Fig. 1. These 

y 
— i] 

y’ 
x 
FIG. 1 CO-ORDINATE SYSTEM — 


quantities are obtained in terms of the molecular velocity 
distribution function, f(u, v, w), which is the number density 
per unit volume of phase space of molecules with velocity 


components, u,v, w. One has 


mf f- [2.1] 


where m is the molecular mass. The velocity components 
u,v, wand u’, v’, w’ are related (see Fig. 1) by 


Pi 


Ti 


= —u'cos@ + v’ sin 


u 
v = —u’ sin9 — v' cos 
w= w’ \ 


where @ is the angle between the (outward) normal to the sur- 
face and the x (flow)-direction. 
The function f is given in the Maxwellian case by 


fi »—((u — U)? + v? + w?)/2RT 2.4 


where p, 'T, U, R are the gas density, temperature, flow 
velocity, and the gas constant, respectively. For the non- 
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Maxwellian case, one uses the distribution function considered 
by Maxwell (7) and Grad (8) 


= fl] 


— Ul — — Uw — 27,00 


(u — U)? + v? + w? 
2{1- 


where pzz, Try, Gz, etc., are the viscous stress and heat flux 
components. These quantities, together with p, U, and T, 
are treated as constants throughout the region occupied by 
the cylinder. For the equilibrium case, and for small gra- 
dients of velocity and temperature, these terms are related to 
the velocity and temperature by the Navier-Stokes relations 


@ 
Try + ] ... [2.6] 
oy Or 
6 
or 


For rapid processes (e.g., shock waves) or for large gradients 
of temperature and velocity, more complicated relations 
exist; see (8) for details. In the present report, the results 
will be obtained in terms of the variables listed in Equation 
2.5], so that they are independent of whether or not Equa- 
tions [2.6] are valid. It should also be mentioned that for 
sufficiently nonuniform situations the approximation of 
Equation [2.5] itself becomes poor. Additional terms in the 
expression for the molecular velocity distribution would then 
be required and would contribute to the gross force exerted 
upon the body. The relative dependence upon the parame- 
ters Prr, Try, Yr, ete., however, would be unchanged. 


2.2. The Reflected Molecules 


The determination of dF,, the force on dA due to the re- 
emitted molecules, requires a specification of the nature of 
the interaction of the impinging molecules with the surface. 
The treatment of this question will be somewhat generalized 
in this report from the usual treatment of Maxwell (7) which 
has been followed in the subsequent work of Tsien (1), 
Stalder and Zurick (2), and others. Traditionally, the inter- 
action has been described in terms of two parameters, the 
thermal accommodation coefficient a and the specular re- 
flection coefficient ¢. The quantity a is defined by 
dE; — dk, 


where dEH;, dE,, and dE, are the energy fluxes, respectively, 
incident on and re-emitted from the surface, and the flux 
which would be re-emitted if all incident molecules were re-— 
emitted with a Maxwellian distribution corresponding to the 
surface temperature 7... The quantity a then is a measure of | 
the degree to which the incident molecules are “accom-_ 
modated”’ to the surface. Measured values of a are usually 
slightly less than unity (9). It will be observed that all 
quantities on the right side of Equation [2.7] are well-de-— 
fined and can be either calculated or measured. 

The definition of ¢, however, is somewhat less satisfactory. — 
It is supposed that the incident molecules are either reflected 
“specularly”’ (i.e., with simple reversal of normal velocity) — 
or “diffusely”? (which for some applications can be just “ran- — 
domly,’”’ but eventually must be defined more precisely as 
Maxwellian corresponding to a temperature 7, which is not 
necessarily the same as 7',). The coefficient o is then defined — 
as the fraction of diffusely reflected molecules, 1 — o being the — 
fraction which is specularly reflected. The definition in this ss 
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Py,(—u’ sin — v’ cos 0)? + paw” 


For completely specular reflection one has a = ¢ = o’ = 
1 


(—u’ cos + sin @ — U)? + (— u’ sin — v’ cos 6)? + 


form of the quantity « thus becomes meaningless unless the 
actual interaction with the surface is of the type specified. 
As measured values of o have actually been determined, how- 
ever, this difficulty has been avoided because only gross 
tangential forces are dealt with. Thus, according to the 
definition of o, one has 


where 7; and 7, are the incident and reflected fluxes of tan- 
gential momentum (shear stress). For all practical purposes, 


_ Equation [2.8], rather than the specular-diffuse model, has 
4 


been the defining equation for ¢. The two definitions are in 
agreement if the interaction actually is of the specular-dif- 


fuse type (which, of course, it almost certainly is not), other- 


wise Equation [2.8] defines o. Equation [2.8] may be 
written also, for symmetry with Equation [2.7], in the form 


% 


_ where 7,, is defined as the tangential momentum flux which 


would be reflected from the surface if all incident molecules 
were re-emitted with a Maxwellian distribution corresponding 
to the surface temperature 7... Clearly 7. = 0. 

Thus the difficulty inherent in the original definition of 
is avoided, in so far as the determination of tangential forces. 
However, it oceurs again in the determination of normal 
forces, and again recourse has been taken in (1) and (2), 
for example, to the specular-diffuse model. If the actual 


interaction is different, then the determination of heat trans- 
fer, tangential] force, and normal force in terms of just the two 
parameters a and o would yield results which one would 


expect to be mutually inconsistent. This difficulty can be 
avoided if an additional parameter o’, say, is defined by 


where pi, pr, and p.» are, respectively, the incident and re- 


flected fluxes of normal momentum (pressure), and the flux 
which would be re-emitted if all the incident molecules were 
re-emitted with a Maxwellian distribution corresponding to 
the surface temperature T7’,,. tre 
The quantity p.. is given (2) by : 


where N,, is the number of molecules per unit time incident 
 onaunit surface area 


It is given by 


Ny = fr u’fdu'dv'dw’ 


hte “| du'dv'dw’ .-[2.16) 


in terms of the drag force Fp and the lift foree Fz. These are 
given by 

Fo = (—p cos 6 + 7-sin 6)rdé 


With the help of Equations [2.9] and [2.10] these reduce to 


For completely diffuse reflection, a = 1. Forth 
hypothetical partly diffuse, partly specular interaction, an( 


for a surface at rest with respect to the gas, one would have 


oVT: V oo — a)T; + cal y 
= 


—>aasT,—T7; 


13} 
However, in general, for actual physical interaction the three 
parameters would be independent. They are sufficient, how- 
ever, to determine free molecule flow heat transfer and ero. 
dynamic force characteristics.‘ 


‘ Strictly speaking, it would be logical for completeness io in- 
troduce two additional interaction parameters, one for possible 
asymmetry in the momentum transfer in the direction mutually 
perpendicular to the normal and tangential directions, ani! the 
other to allow for possible net mass interaction, i.e., ingassing or 
outgassing. It is tacitly assumed in the present report that such 
additional phenomena are not present. 


2.3. The Combined Expressions 

The total force on a unit of surface area is analyzed int 
a normal component p and a tangential component 7 given by 


r = — Tr) 


[2.14] 


T= 


(2 )pi + 
oT; \ 
In terms of f, using Equations [2.11], [2.12], [2.1], and [ 2.2), 
one has 


(2.15 


The total force on a unit length of the eylinder is obtained 


2x 
(—p sin 6 — 7 cos 6)rdé 


Combining Equations [2.17], [2.16], [2.5], and [2.3], one has 
finally 


cos 6 + ou’v’ sin x 


2RT 


— 27:,(—u’ cos @ + v’ sin @ — U)(—w’ sin @ — v’ cos 8) — 272.(—u’ cos + 


(—u’ cos @ + v’ sin@ — U)? + (—w’ sin @ — v’ cos 6)? + w’ 


Ee u’cos@ + v' sing — U)? + 


1 
( 2pRT 


v’ sin§@ — U)w’ — sin @ — v’ cosé)w’ — 2 


(—u’ cos 6 + v’ sin @ — U)? + (—u’ sin 6 — v’ cos 6)? 


—w cos 6 + v’ sin@d — U) + q,(—u’ sin @ — v’ cos 6) + . [2.18 


+ 


2RT 
Py(—u’ sin @ — v’ cos 0)? + 


v’ sin @ — — 2 7,(—u’ sin — v’ cos — 2 (: 


— 27:,(—u’ cos @ + v’ sine — U)(—u’ sin — v’ cos — 2 
(—u’ cos @ + v’ sing — U)? + (—wu’ sin 6 — v’ cos 0)? + w 


5RT 


2 
E —u’ cos + v’ sin + 


—u' cos? + v’ sind — U) + wu’ sin — v’ cos 6) + we’) rdu'dv'dw’d@. . . [2.19] 


) x 
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Performing the integration indicated in Equations [2.18] 
and [2.19], one obtains 
pr 


(3 + 


Fp 


6S 2 
(n(—)-1 (= = (an (— 2+— 1-—)a(- 
(4 — 20’ + (<) ( 6qy =) (=) 
where Jy and J, are the modified Bessel functions, and s is these viscous stress terms wil] be small compared to the force 
the molecular speed ratio due to the uniform flow (of the order of a few per cent or less). 
: U ; These quantities might thus be neglected in connection with 
[3.3] use of the cylinder as a probe. 


Ye 

It will be observed that the drag force for the nonuniform 
case differs from that for the uniform case by terms propor- 
tions], respectively, to the heat flux in the x-direction and to 
the deviation of the normal stresses from the hydrostatic. 
The lift foree, which is zero for the uniform case, contains 
terms proportional to the shear stress and the heat flux in 
the y-direction. 

Equations [3.1] and [3.2] may be rewritten in the form 


= Cp, + + Coprr + 
Fr Try qu 

= — C 3.5 

I ey + aU’ CL, [3.5] 


The indicated partial lift and drag coefficients, Cx,,,, Cog, 
etc.. correspond to the contribution to the total lift or drag 
arising from the nonuniformity denoted in the subscript. 
These results are presented graphically in Figs. 2 and 3 for 
the case of diffuse reflection. It should also be mentioned 
that the results of Equations [3.1] and [3.2] agree with those 
of (2) for the uniform case. 

Inspection of Figs. 2 and 3 reveals that the partial lift 
and drag coefficients corresponding to the nonuniform terms 
are all small compared to the coefficient for the uniform case 
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PARTIAL DRAG COEFFICIENT FOR NONUNIFORM FREE MOLE- 
CULE FLOW PAST A CYLINDER (DIFFUSE REFLECTION ) 
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+ (1 + 28?)/; 


Copp. Unless a gas flow is very nonuniform, i.e., if the flow 
itself is in the continuum or slip flow region, the quantities 
Pzr, Pyy, and rz, are small compared to p. It follows that the 
contribution to the total force on the cylinder arising from 


4s, 60S 


The relative importance of the forces arising from the 
heat flux terms, however, depends on the magnitude of the 
ratio of g, or g, to pU. At high flow velocities these forces 
wil] thus be negligible compared to the force for the uniform 
case. At very low velocities they will not. In a stagnation 
region such forces constitute the entire force on the cylinder. 


4.1. In the free molecule flow of a nonuniform gas (i.e., 
one in which heat flux or viscous stresses are present) past a 
cylinder, both the lift and drag forces will be affected by the 
nonuniformity. 

4.2. Force contributions due to stress and heat flux terms 
are, in general, small compared to the force for the uniform 
case, except in regions of low flow velocity where the forces 
arising from the heat flux terms become important. 


Conclusions 


4.0 
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1 ‘Superaerodynamics, Mechanics of Rarefied Gases,’ by 
H. S. Tsien, Journal of the Aeronautical Sciences, vol. 13, Decem- 
ber 1946, pp. 643-664. 
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Rockets 


Navy Bureau of Aeronautics announced that 
smokeless JATO units are now in full production with 


-— : of deliveries currently being made to the Navy and the Air 
Force. 

«TT HAS been reported that the Chrysler Corporation has 


reeeived a contract for the production of the ‘Redstone’’ 
_ surface-to-surface missile at Detroit, Mich. The missile 
- was developed at the Army Ordnance Redstone Arsenal at 
—- Huntsville, Ala., by a group of former German engineers 
- from the V-2 group at Peenemuende. The group is headed by 
Wernher von Braun. 


FAIRCHILD’S “Lark’’ guided missile is shown (Fig. 1) 
USS Norton Sound. This photograph is the first cleared for 
public release and shows the firing of this anti-aircraft type 


guided missile from a Navy vessel. The Lark is produced 


Solar Aircraft Compar 


FIG. 1 SHIPBOARD LAUNCHING OF FAIRCHILD “LARK” 


by the Guided Missiles Division of the Fairchild Engine and 
Aviation Corporation at Wyandanch, Long Island, N. Y. 
_ Currently Larks are being flown by all three branches of the 
err Services in their training and evaluation programs. 
- Design of the missile was begun under a Navy project during 
the latter months of World War II as a defense against 
_ Japanese suicide planes. Under the Lark project in the years 
; after the war, Fairchild engineers developed a homing system 
- considered one of the most advanced yet devised. The cloud 
_of smoke partially obscuring the stern of the ship is emitted 
by the booster rockets which push the “bird’’ to flight speed 
and then drop off as the missile continues under the thrust of 
own rocket motors. 
A TOTAL of 100 launchings of the Aerojet-General ‘“Aero- 
bee’’ have been conducted from White Sands Proving Ground, 
Holloman Air Development Center, and the USS Norton 
Sound. An additional quantity is on hand to be launched 
in the near future, and production of new units will continue 
this year. At the present time, a new program is under way 
to further extend the altitude and range of the Aerobee. 
Although the design performance of the new ‘‘Aerobee-Hi”’ 


; made in general utility permitting an even wider use. 


THE Air Force is modifying 6 Martin B-61 Matador pilot- 
less missiles for testing and training programs associated with 
the Hughes “Falcon,” and also for ground-to-air weapons, 
pilotless missiles, and rockets. 


AN all-rocket power plant producing 6000 Ib of thrust and 
weighing only 210 lb has been announced by Reaction Motors, 
Inc. The rocket engine, designated as 1500N4C by RMI, 
burns alcohol and liquid oxygen. Similar units have powered 
jet-fighter aircraft of the Navy and Air Force and single- 
stage missiles. 


PRODUCING 1000 Ib of thrust for 14 sec, the Aerojet- 
General Corporation 14AS-1000G-1 solid propellant JATO 
has been certified by CAA for use in commercial aircraft. 


leaving the launcher on the Navy’s guided missile ship, the — 


Novel Jet Devices 


a ila Nagler Helicopter Company of Westchester County 
Airport, N. Y., has been using rocket motors in autogyro ap- 
plications. The rocket motors, located at the hub of the 
blades, are discharged through a nozzle at the tip of the blades. 
Two models have been displayed. One is called a Converti- 
craft, an Aeronea light plane fitted with rocket-powered auto- 
gyro blades. Another is called the Heliglider, a one-man lift- 
ing unit. The strap-on unit, weighing 65 lb, will lift a man 
off the ground at 2700 fpm climb rate to a ceiling of 2700 ft. 
Range is about 5 miles. 


ANOTHER novel use for reaction engines involves the use 
of the pulsejet engine which was used on the German \-1 
“buzz-bomb.”’ Devenco of New York uses the pulsejet 
engine to atomize insecticides and fumigants, to produce 
smokes and fogs, and for use as a blowtorch to kill weeds or 


to remove snow and ice. 
'Furbojet Engines 


A JET thrust reverser, developed by the French national 
aircraft engine facility at Snecma, has been flight-tested in a 
de Havilland Vampire in France. Reverse thrust is said to 
be about half the normal positive thrust of the engine. A 
compressor bleed of air is injected into the center of the en- 
gine exhaust flow. The exhaust flow is thus forced into a 
series of 20 annular deflector vanes placed concentrically 
along the thrust line at the tailpipe. The flow is thus turned. 
A requirement of about 2 per cent of the engine air flow is 
needed to produce the negative thrust which acts as a brake. 
However, as much as 10 per cent thrust reduction has been ex- 
perienced when the device is not in operation due to air bleed- 
ing through the annular vanes. 


A DEVICE to make the operation of high-speed jet: :ir- 
craft easier and safer has been announced by Solar Aircraft 
Company of San Diego, Calif. Called the Microjet, 
it directly senses engine pressures and replaces complex 
electronic equipment. Although most details are under 
security restrictions, the device is stated to be pneumatically 
operated. It computes what the turbine discharge pressure 
should be under all flight conditions, and at the same time 
notes any error between the actual engine pressure and what 
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it should be. If there is an error, electrical signals are auto- 
matically sent out to other engine controls which correct the 
pressure conditions, 


THE Westinghouse Electric Corporation and Rolls-Royce, 
Ltd., of Derby, England, have signed an agreement provid- 
ing for technical co-operation between the two companies for 
a period of ten years. The agreement includes exchange of 
information in the design, development, and production of 
gas turbine and aero-engines. The first application will be 
in the collaboration of Rolls-Royce with Westinghouse in 
future development of Westinghouse jet engines. 


Courtesy El Al Israel Airlines 


EL AL ISRAEL AIRLINES’ CURTISS C-46 SHOWING TWIN JETS 
INSTALLED UNDERNEATH THE FUSELAGE 


FIG. 2 


DR. ERICH SCHATZKI, aeronautical engineer and de- 
signer, suggested that the installation of jet engines as thrust 
boosters on Curtiss C-46 aircraft, developed by him for El Al 
Israel Airlines, might offer the most promising solution in 
improving C-46 performance in compliance with CAA 
regulations. Curtiss C-46 aircraft would be equipped with 
two Marbore II engines each. Both engines would be housed 
in individual, interchangeable pods (Fig. 2) underneath the 
fuselage. No structural modifications are necessary. 


IT IS understood that Allison is preparing to flight-test its 
J-71 turbojet engine in a B-45 flying test bed. The J-71 isa 
9600-lb thrust engine having a 16-stage axial-flow compressor 
and a three-stage turbine. This engine has recently com- 
pleted a 200-hr test run in an altitude test chamber of the 
NACA. Reports in the press state the J-71 will be used in 
the Douglas B-66, the Republic F-105, the Northrop F-89, 
and the McDonnell F-3H. 


THE de Havilland Aircraft Company is reported to have 
developed a new lightweight jet engine, ‘“‘Gyron,’’ said to be 
Britain’s first in the 15,000-lb to 20,000-lb thrust class. 
This engine employing a seven-stage, single-spool compressor 
was developed for around $2.5 million under a joint de Havil- 
land-General Electric private venture. 


ROLLS-ROYCE has installed ‘“teeth’’ around the perime- 
ter of a jet engine tailpipe and has reduced low-frequency 
noise normally encountered in engine tests, although there 
has been some increase reported in high-frequency noise. 


Aircraft, United States 


THE Strategic Air Command has disclosed that the B-47B 
Stratojet recently flew more than 12,000 miles nonstop in 24 
hours. Aerial refueling techniques were employed. 


PRATT & Whitney’s J-57 axial-flow gas turbines will 
power Boeing’s new prototype jet transport. The engines 
are similar to the J-57’s used on the B-52 long-range bomber. 
Boeing expects its jet transport to be in the air by mid-1954. 
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Called Project X by Boeing, the four-engine, swept-wing jet 
liner is now called Model 707. The aircraft is expected to 
have performance characteristics which will give it a definite 
speed and power edge over the de Havilland Comet. The 
707 is to carry 100 passengers at a speed of nearly 600 mph. 


THE first flight of the North American YF-100 has been re- 
ported to have taken place at Edwards AFB, Calif. Sweep- 
back is said to be 45 deg. Other features are a streamline ex- 
ternal fuel tank of new design, canopy design for high speed, 
and a lowered horizontal tail fin. The F-100 has been 
ordered into quantity production by the USAF. The weight 
is expected to be around 24,000 lb. Power is supplied by a 
Pratt & Whitney J-57 engine delivering 15,000-lb thrust with 
an afterburner. 


FIG. 3 shows a view of the XA3D, the Navy’s new swept- 
wing jet attack aircraft, now being built by the Douglas Air- 
craft Company. It is powered by two jets, each slung in a 
pod under the wing outboard of the fuselage. It will be in 
the 600 to 700-mph class and will carry a crew of three. As 
in other carrier aircraft, the wing of the A3D will fold to per- 
mit easier handling and more compact storage aboard air- 
craft carriers. Detailed performance data are secret, but no 
known airplane of comparable size, now in service or con- 
templated for early service introduction, can carry an equiva- 
lent bomb load as high and as fast as the A3D. 


Courtesy U.S. Navy 


THE XA3D, THE NAVY’S NEW SWEPT-WING JET ATTACK 
AIRCRAFT 


FIG. 3 


A NEW world speed record made on July 16 is claimed for 
the F-86D which flew 715.7 mph. A General Electric J47-17 
turbojet engine powered the Air Force’s North American 
Sabre Jet. The J47-17 engine is equipped with electronic 
controls and an afterburner which utilizes additional fuel for 
extra thrust. 


FIG. 4 shows the new F-94C Starfire, all rocket-armed in- 
terceptor now rolling off Lockheed Aircraft Corporation pro- 
duction lines. Besides the standard 24 rockets housed in its 
sleek nose, the USAF plane now carries a dozen 2.75-in. 


4 


"Courtesy Lockheed Aircraft Corporation 


FIG. LOCKHEED’S NEW F-94C PRODUCTION LINE 
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rockets on each wing in pods visible on the plane at left, 
where a workman is ready to install the rocket launcher’s 
white plastic breakaway cover. The stovepipe apparatus 


i behind each F-94C is a silencer to muffle engine sound during 


ground tests. 


FIG. 5 shows standard HVAR 5-in. rockets mounted in the 
launching position on the rotary bomb door developed by 
The new preloaded 


- Martin for the XB-51 tactical bomber. 


Courtesy Glenn L. Martin Company 


FIG. 5 HVAR 5-IN. ROCKETS MOUNTED IN LAUNCHING POSITION 
ON ROTARY BOMB DOOR 


door rotates 180 deg just prior to bomb or rocket release. 
Through use of this door, a jet bomber traveling at high speed 
is just as stable a platform as one flying at World War II 
speeds, and U. S. military aircraft will not have to slow down 
to make @ successful bomb or rocket release as they have in 
the past. 


REPUBLIC Aircraft Corporation has disclosed details on 
two of its aircraft. The first of the new aircraft is the 
Thunderstreak, a sleek new high-speed swept-wing F-84F 
jet fighter, which is being produced in quantity for the U. 8. 
Air Force and NATO countries. Although much of the F- 
84’s performance is classified, it is listed as a plane of more 
maneuverability and higher speed than the Thunderjet. 
General data: span 33 ft, 6in.; length 43 ft, 4 in.; height 14 
ft, 4 in. at rudder tip; speed 600 mph plus; combat radius: 
not disclosed but exceeds F-84G, which with external tanks is 
more than 1000 miles; service ceiling 45,000 ft plus; arma- 


Courtesy Republic Aviation Corporation 


FIG. 6 REPUBLIC XF-91, HIGH-ALTITUDE INTERCEPTOR FIGHTER 
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ment six 0.50-caliber machine guns plus externally mounted 
rockets, bombs, and napalm; power Wright-Sapphire J-6) 
turbojet engine developing 7200-lb thrust. The Republic 
XF-91 (Fig. 6) is a high-speed, high-altitude intereeptor 
fighter which has made several flights at speeds faster than 
sound to become the first combat-type aircraft with this 
achievement. Incorporating several radical features, the XF- 
91 is powered by a General Electric J-47 5200-lb thrust turbojet 
engine plus afterburner, and emergency power is supplied by a 
6000-lb thrust rocket engine built by Reaction Motors. An 
unusual wing design features inverse taper, sweptback wings, 
also a variable incidence of the wing permitting changes in 


angle of attack. ; 


THE RAF recently flew a Supermarine Swift fighter, the 
Swift F.4, which boasts more power and more armaments t)ian 
the early Swift F.1. The Rolls-Royce Avon engine with 
afterburner has a thrust of about 9000 lb. Carried are tour 
30-mm cannon. 


Aircraft, Foreign 


A COMET jet airliner, the first to visit the North American 
continent, arrived at Ottawa on May 29, having flown from 
London across the Atlantic via Iceland and Labrador. ‘The 
total flying time for the journey of 3550 miles was 10 hr, 
48 min. 


RECENT movies show launchings of the British-made 
Fairey vertical take-off aircraft at the Woomera missile range 
in Australia. This delta wing craft is powered by a Fairey 
Beta I rocket engine having two 900-lb thrust motors. A 
booster assisted take-off is produced by two 600-lb thrust 
rockets. Additional information on the rocket propulsion 
system or the flight test data is not available at this time. 


THE Jindivik pilotless aircraft, which is radio-controlled 
and is used as a target drone, has been placed into operation 
at the Woomera rocket range. Take-off is from a tricycle 


dolly, landing with a retractable skid. The aircraft is 23 ft 
long and has a wing span of 19 ft; powered by an Armstrong- 
Siddeley turbojet engine, two models are disclosed, one hav- 
ing a thrust of 1000 lb, another 1600 Ib. Use is as a target 
for other weapons now in development at Woomera. 


Courtesy Aviation Week Magazine 


FIG. 7 THE SUD-OUEST 8S.0. 4050 VAUTOUR 


THE French Sud-Ouest 8.0. 4050 Vautour is shown in 
Fig. 7. Two 5500-lb thrust Atar jets mounted in underwing 
nacelles give it sonie-speed capabilities in level flight. 


FLIGHT tests are reported to be under way with a new 
Ledue ramjet aircraft, the Leduc 021 which is being fostered 
by the French Air Ministry. The craft is said to weigh 
about 11,000 lb and capable of a speed of Mach 0.95 in level 
flight. New features include retractable outrigger wheels 
with tandem main wheels, tiptanks, and a jettisonable cockpit 
located upstream from the air inlet. 
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THE French Snease 8.E. 5000 Baroudeur jet ground-attack 
plane is capable of taking off from a small field without a run- 
way by means of a rocket-propelled sled which is jettisoned 
after the craft is airborne. Take-off space is thus claimed to 
be under 700 yards. Landing is with retractable skids. 
The wings of the plane are sweptback 35 deg, which is ex- 
pected to allow sonic flight capabilities with its Atar 101 
power plant of 6000-Ib thrust. 


High Altitude and Extraterrestrial 
Research 


TO determine the effects of cosmic radiation at altitudes of 
50,000 to 100,009 ft, the Air Research and Development Com- 
mand is sending fruit flies, sealed in pressurized tubes, aloft 
in controlled balloons. Launchings will be from the West 
Coast and drifts will be eastward across the U.S. Data from 
these tests are expected to aid our knowledge of the effects of 
cosmic radiation for high-altitude aircraft. 


THE U.S. Navy operating in Baffin Bay, Northern Green- 
land, has launched rockets from balloons to secure cosmic 
radiation data. ‘Deacon’ type solid propellant rockets 
weighing about 300 lb are launched from ‘‘Skyhook’”’ balloons 
and have attained altitudes of about 40 miles. 


ANOTHER high-altitude project using the Viking rocket 
calls for the study of the Aurora Borealis from the USS Norton 
Sound in co-operation with Norwegian scientists at Tromso. 
Dr. H. Newell is in charge of the program. = 


Courtesy Dou craft Company 


FIG. 8 “HELMET DESIGNED FOR EMERGENCY ESCAPE FROM AIR- 
CRAFT FLYING AT SUPERSONIC SPEEDS 


DEVELOPMENT of a new slotted helmet designed to 
protect pilots forced to bail out from their aircraft while 
traveling at supersonic speeds, was announced by the U. S. 
Air Force Research and Development Command (see Fig. 
8). Engineers of the Douglas Aircraft Company, Santa 
Monica Division, who designed and built the helmet, said 
the slots or vents which are cut into the forward crown section 
of the headgear greatly reduce wind shock and air lift and 
keep the helmet secure on the pilot’s head. The helmet has 
been successfully tested in an outdoor wind tunnel at simu- 
lated speeds up to Mach 1.04. 


SAINT Louis University will become one of the first insti- 
tutions in the world to offer a program of graduate study 
dealing with the problems of space travel. The program will 
be offered through Saint Louis University’s Parks College of 
Aeronautical Technology. The first phase of the program will 
be a seminar offered to senior students in aeronautical en- 
gineering and will deal with rocket stability and design. 
Faculty members of Saint Louis University w]li conduct the 
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seminar, but it is anticipated that leading authorities in the 
field will be brought in as guest lecturers. 

THE Naval Aero Medical Laboratory at Philadelphia has 
tested its new “space suit’’ and will soon conduct tests aloft. 
Made of rubber with a large plexiglas helmet fastened to the 
neck, the suit has been successfully tested at simulated alti- 
tudes up to 70,000 ft. 


Facilities and Equipment 


HUGHES Aircraft Company is reported to have avionics 
(electronics in aviation) contracts totaling $200 million per 
year with a $600 million backlog. Research and develop- 
ment on the Hughes air-to-air Falcon missile, F-98 (carried by 
the Convair F-102 interceptor), led to development of 
avionic control and guidance systems such as the radar fire- 
control system used with the F-86D, F-89D, and F-94C, 
the F2H-4, and the Canadian CF-100. 


NACA has disclosed details of its new 8-ft wind tunnel at 
Langley Field, Va. Capabilities are reported to be Mach 1.4 
with variable speed, temperature, and air pressure. It is ex- 
pected that the tunnel will be in normal operation in about a 
year. Also at Langley is the new Gas Dynamics Laboratory 
used to gain insight into the physics of high-speed flight. 
Flight speeds up to Mach 9.0 or at altitudes of about 200,000 
ft can be simulated. 


G. M. GIANNINI & Co., designers and manufacturers of 
precision instruments used in guided missiles and radar control, 
plans to build a new 24,000-sq ft assembly plant in Pasadena, 
Calif. 


ENGINEERING Research Associates has been made a 
division of Remington Rand, Inc. The country’s largest 
manufacturers of electronic computers have facilities at St. 
Paul, Minn., and Arlington, Va. 


CONSOLIDATED Engineering Corporation of Pasadena, 
Calif., is planning a 20-acre ‘Instrument Park” in the Hast- 
ings Ranch area of Pasadena. The Corporation is a manu- 
facturer of electronic analytical instruments. 


REACTION Motors of Rockaway, N. J., has announced 
eround-breaking for a new plant site. In addition, it was re- 
vealed that a backlog of over $5.5 million exists. 


CONSTRUCTION proceeded from the roof down on Boe- 
ing’s new B-52 hangar at Seattle. The roof structure, 90 
ft above ground level, rested on wood towers while supporting 
steel columns were built beneath it. Another feature is its 
785 ft long, 65 ft high door—the largest unobstructed door in 
the country. Four of the all-jet B-52’s ean be housed in the 


hangar. 


PRODUCTION of 8 to 12-ton capacity liquid oxygen and 
nitrogen transfer, storage, and transport units has been an- 
nounced by Hofman Laboratories, Newark, N. J. A semi- 
trailer consists primarily of a liquid storage vessel and an 
outer shell which also functions as part of the chassis. All 
components are attached to the container for accessibility 
and easy operation. A vacuum space between the container 
walls maintains high thermal efficiency. While the trailer is 
at rest the evaporation loss will not exceed 250 Ib per day. 
An electric-motor-driven centrifugal liquid-oxygen pump, 
‘apable of transferring 100-150 gal of liquid oxygen per min, 
is installed at the rear of the trailer with all other controls. 


THE Components Development Section of the Aircraft 
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Courtesy General Elec 
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SHOCK TUBE FOR STUDYING JET ENGINE COMBUSTION 
PROCESSES 


FIG. 9 


Gas Turbine Division of the General Electric Corporation 
(Evendale, Ohio) is using shock waves to study combustion 
processes (Fig. 9). The effect of such waves is studied in a 
miniature shock tube. Development on jet engine com- 
bustion and the effect of shock waves are becoming increas- 
ingly important, G. E. engineers explain, since the combus- 
tion chamber is one of the most important sections of a jet 
engine. 


MATHIESON Chemieal Corporation has announced it is 
producing hydrazine at its new three-million-dollar plant at 
Lake Charles, La. The principal use of the hydrazine will 
be by the armed forces in various rocket programs. Hy- 
drazine will be prepared by a modified Raschig Process 
which utilizes the reaction of ammonia with sodium hypo- 
chlorite to form monochloramine; the latter is reacted with 
ammonia to form hydrazine. Hydrazine, a powerful reduc- 
ing agent, now sells for about $3 per lb; the price is expected 
to fall to about $0.50 per lb with large-scale production and 
improved techniques. 


THE Armour Research Foundation of the Illinois Institute 
of Technology, Chicago, IIl., is currently engaged in several 
jet projects. These include: effects of explosions in jet- 
engine test cells, analysis of forces and surface pressures in jet- 
engine fuel pumps, instrumentation and resistance probe 
for true free-air thermometer, guided missile warheads, 
rocket-assist artillery ammunition, physical constants of 
propellant materials, and a bent-nozzle isotope separator. 


_ High-Frequency Combustion Instability in Rocket Motor with 


Concentrated Combustion 
4 “Aspects of Combustion Stability in Liquid Propellant 


(Continued from page 313) 
Rocket Motors, Parts I, II,” by L. Croeco, JouRNAL OF THE 
AMERICAN Rocket Society, vol. 21, November 1951, pp. 163- 
178; vol. 22, January-February 1953, pp. 7-16. 

5 “High Frequency Combustion Instability in Rockets With 
Distributed Combustion,”’ by L. Crocco and Sin-I Cheng, paper 
presented at the Fourth Symposium on Combustion, September 
1952. 

6 “Supercritical Gaseous Discharge with High Frequency 
Oscillation,’ by L. Crocco, paper presented at the Eighth In- 
ternational Congress of Applied Mathematics and Mechanics, 
August 1952. 
7 “The Transfer Functions of Rocket Nozzles,” by H. S. 
Tsien, JOURNAL OF THE AMERICAN Rocket Society, vol. 22, 
May-June 1952, pp. 139-143. 
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on Miniature 
Ball Bearings 


_ Ultra small precision ball bearings supply accurate 
alignment, long life and attention-free operation to 
aviation instruments, railway and marine indicating and 
recording devices, medical appliances, and similar fine 
mechanisms. 


Complete specifications are given in new Catalog 
52B. Engineering data sheets, speed load charts, and 
measuring unit conversion tables are available. For these 
and for further information write or phone H. D. Gilbert 
at the plant. 


Miniature 
Incorporated 


(Precision. Bearings 
AS Keene, New Hampshire 


Pioneer Precisionists to World's 
Foremost Instrument Manufacturers 


SAVE SPACE, 
WEIGHT, FRICTION 


Aerodynamic Forces on a Cylinder for the Free Molecule Flow of a 


Nonuniform Gas 


(Continued from page 317) 

“Theoretical Aerodynamic Characteristics of Bodies in a 
ies Molecule Flow Field,”’ by J. R. Stalder and V. J. Zurick, 
NACA Technical Note 2423, 1951. 

3 “A Comparison of Theory and Experiment for High Speed 
Free Molecule Flow,” by J. R. Stalder, G. Goodwin, and M. 0. 
Creager, NACA Technical Note 2244, 1950. 

4 “A Torsion Balance for Measuring Forces in Low Density 
Gas Flows,” by I. Estermann and E. D. Kane, Journal of Applied 
Physics, vol. 20, June 1949, pp. 608-610. 

5 “Zur Theorie der Radiometerkrifte,” by 
Z. Physik, vol. 27, 1924, p. 1. 

6 “Zur Theorie des Radiometers,”’ 
vol. 54, 1929, p. 537. 

7 “On Stresses in Rarefied Gases Arising from Inequalities 
of Temperature,’’ by J. C. Maxwell, Transactions of the Royal 
Society of London, vol. 170, 1879, p. 231. 

8 “On the Kinetic Theory of Rarefied Gases,” by H. Grad, 
Communications on Pure and Applied Mathematics, vol. 2, 1949, 
pp. 331-407. 

9 ‘Kinetic Theory of Gases,” by L. B. Loeb, McGraw-Hill 
Book Co., Inc., New York, N.Y., 1934. 

10 “Theory of Drag in Highly Rarefied Gases,” by M. Hei ‘¢ - 
man, Communications on Pure and Applied Mathematics, vol. 
1948, pp. 259-273. 

11 “Applications of the Theory of Free Molecule Flow to 
Aeronautics,’ by H. Ashley, Journal of the Aeronautical Sciences, 
vol. 16, 1949, pp. 95-104. 

12 “Transport Phenomena in Very Dilute Gases,”’ 
Wang Chang, NORD 7924-U17H-3-5, 1950. 

13. “Wind Tunnel Testing Problems in Superaerodynamics, ”’ 
by H. 8S. Tsien, Journal of the Aeronautical Sciences, vol. 15, Oct. 
1948, p. 577. 


A. Einstein, 


by P. Epstein, Z. Physik, 


by C.8 
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H. K. WILGUS, Associate Editor 


Successful ARS Summer Meeting at Los Angeles, 


me July 2, Draws Record Attendance = 


Southern California Section Acts as Host 


ie conjunction with the ASME Summer 

Meeting held at the Statler Hotel, Los 
Angeles, Calif., June 30-July 2, 1953, the 
American Rocket Society sponsored two 
technical sessions and a luncheon on the 
last day of the meeting. More than 150 
ARS members registered, with the Sou- 
thern California Section acting as unoffi- 
cial host and participating in the program. 


Morning Session 


The morning technical session, at which 
four papers were presented, attracted 
more than 100 members. The chairman 
and co-chairman for this session were 
R. D. Geckler, chief technical specialist, 
Aerojet-General Corporation, Azusa, 
Calif., and Leon Green, Jr., senior engi- 
neer, Solid Engine and Chemical Division, 
Aecrojet-General Corporation. 

The following are brief summaries of 
the papers presented: 

“Rocket-Powered Wind Tunnel” by 
F. Kreith, P. B. Stewart, and E. S. 
Starkman, University of California, Ber- 
keley, Calif. Aerodynamic performance 
testing sometimes requires large capacity, 
high velocity, and subatmospheric pres- 
sures. Not many large-scale installations 
are available which satisfy such require- 
ments, primarily because of the cost of the 
blowers and evacuating pumps required. 
A low-cost design is proposed which uses 
the exhaust from rocket engine thrust 
chambers as motive fluid for jet pumps 
built integrally into the wind tunnel. 
Configurations, estimated performance, 
and cost analysis were presented for a 
large, high subsonic wind tunnel capable 
of operating at altitude conditions from 
5000 ft to about 60,000 ft. 

“Correlation of Experimental Data on 
the Disintegration of Liquid Jets’ by 
C. C. Miesse, Aerojet-General Corpora- 
tion, Azusa, Calif. The paper summarized 
the various theories of jet disintegration, 
from Rayleigh’s first analysis to contem- 
porary theories, evaluating them accord- 
ing to their ability to effect a correlation 
of the wave length, maximum drop size, 
and breakup-length data obtained from 
over 60 photographs of jet disintegration. 
These photographs included tests with 
two different liquids (water and liquid 
nitrogen), two different types of nozzles, 
and a wide range of nozzle diameters and 
jet. velocities. 

“New Instruments for Rocket Motor 
Testing’ by R. B. Bowersox, R. A. 
Buchanan, and J. R. Zweizig, Jet Propul- 
sion Laboratory, California Institute of 
Technology, Pasadena, Calif. The paper 
described three instruments developed at 
the Jet Propulsion Laboratory to meet 
special operating requirements in rocket 
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motor operation and _ research tests. 
These instruments are: (1) a ten-channel 
recording potentiometer which makes 20 
measurements per sec with an accuracy of 
+0.3 per cent, and has input-voltage 
ranges of 10, 50, and 100 millivolts; (2) 
a chopper-stabilized d-c amplifier having 
negligible drift and fixed gains of 10, 25, 
50, 100, or 200; (3) a rocket motor per- 
formance computer which computes mix- 
ture ratio and characteristic velocity con- 
tinuously during a motor test with an 
accuracy of +0.5 per cent. 

“Effect of Wave Propagation in Feed 
Lines on Low Frequency Rocket Insta- 
bility” by R. H. Sabersky, California In- 
stitute of Technology, Pasadena, Calif. 
This paper analyzed a liquid monopro- 
pellant rocket system for low frequency 
instability, taking into account the com- 
pressibility of the fluid in the propellant 
line. The propellant is assumed to turn 
into gas at a fixed time interval after 
entering the chamber, this time interval 
being taken to be independent of pressure 
or temperature. A numerical example 
was given, illustrating that the effect of 
compressibility may become important 
for practical cases. 


Major Yeager Speaks at ARS 
Luncheon 

More than 250 members and guests, 
including many representatives of the 
West Coast aircraft industries, attended 
the ARS Luncheon, held in the Golden 
State Room of the Statler Hotel, to hear 
Major C. E. (“‘Chuck’’) Yeager of the AF 
Experimental Flight Test Center, Ed- 
wards AF Base, Calif. 

B. L. Dorman of the Southern California 
Section, manager of the Test and Field 
services Division, Aerojet-General Cor- 
poration, was the presiding chairman. 


J. H. Kindelberger, chairman of the Board, 
North American Aviation, introduced 
Major Yeager who spoke on the topic, 
“Experiences in Rocket-Propelled Air- 
craft.” Major Yeager was the first man 
to penetrate the sound barrier (to fly 
faster than sound). 


Afternoon Session 


The afternoon technical session, at 
which three papers were presented, was 
presided over by 8. K. Hoffman, president 
of the Southern California Section, and 
chief of the Propulsion Section, Aero- 
physics Department, North American 
Aviation. The co-chairman for the ses- 
sion was L. KE. Stocking of Douglas Air- 
craft. 

The first paper presented was “An Air- 
Transportable Liquid-Oxygen Generator 
—Its Operation and Application” by C. A. 
Bleyle, R. B. Hinckley, and C. L. Jewett, 
Arthur B. Little, Inc., Cambridge, Mass. 
It discussed the design and operation of a 
lightweight air-transportable liquid-oxy- 
gen generator, constructed for the U. 8. 
Air Force. This plant, which operates on 
a low-pressure cycle, separates oxygen 
from air at the rate of 10 tons per day. 
Unique design features include a special 
type of heat exchanger, aluminum con- 
struction, gas-turbine-driven air compres- 
sor, and skid-mounted sections for loading 
into standard cargo aircraft. 

The second paper was “Solid Propellant 
Rockets—Basic Concepts, Present Status, 
and Trends in Development” by Charles 
KE. Bartley, Grand Central Aircraft Com- 
pany, Pacoima, Calif. Only the unclassi- 
fied sections of this paper were discussed. 

“Take-Off from Satellite Orbit,” by H. 
S. Tsien, California Institute of Tech- 
nology, Pasadena, Calif., appeared in{the 
JOURNAL, July-August 1953,pp. 233-236. 


Note: Due to lack of space, news ‘of 
ARS Sections will appear in the next%issue 
of the JouRNAL. 


THE LUNCHEON OF THE ARS SUMMER MEETING, HELD IN THE GOLDEN STATE ROOM, STATLER 
HOTEL, LOS ANGELES, CALIF., JULY 2, 1953 
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H. S. SEIFERT, California Institute of Technology, Associate Editor 


RAUMFAHRTFORSCHUNG, H. Gartmann, 
editor, R. Oldenbourg Co., Munich, 
1952, 199 pp. $3. 
Reviewed by Frep WoLF 
California Institute of Technology 
Jet Propulsion Laboratory 


With the co-operation of the Stuttgart 
Society for Space Research the editor has 
arranged for a number of specialists to 
write @ series of rather loosely connected 
essays in German on arbitrarily selected 
topics in the field of space travel. 

Willy Ley, author of the first section, 
starts with the successful flight of the two- 
stage Bumper rocket, which in his opinion 
represents a major historical event. Then 
he reaches into mythological antiquity to 
trace in 20 pages the idea of space 
travel throughout history. He surveys the 
scientific thought immediately preceding 
the modern development of rockets, which 
he dates as beginning with Oberth’s book 
on space travel in 1923. He then roams 
backward through time and _ fictional 
literature searching for the beginning of 
the idea of space travel. 

In the second chapter Werner Schaub 
gives a good 80-page summary of celestial 
mechanics from an astronautical view- 
point. Building up from fundamental 
dynamics, Professor Schaub develops the 
possible orbits and the equations of the 
two- and three-body problem in simple 
steps, interpreting the concise mathema- 
tical formulation in clear uncramped Ger- 
man. Discussion of periodic orbits of the 
restricted three-body problem leads to 
direct applications in astronautics. 

In the third chapter Rolf Engel, U. T. 
Bédewadt, and Kurt Hanisch describe a 
multistep rocket for the purpose of supply- 
ing a satellite space station. The design 
considers two possible orbits for the space 
station, one at 346-mile altitude and the 
other at 1040 miles. The authors briefly 
discuss the effects of altitude and the 
inclination of the orbital plane with the 
ecliptic on the field of view and illumination 
of the station. Apart from this, the main 
emphasis is laid on the supply vehicle 
itself, its trajectory of ascent, structural 
vehicle efficiency, and optimization under 
the numerous conditions imposed on the 
flight. Conservative’ assumptions have 
been made for the propellant performance, 
and the data of the aircraft industry on 
structures were used to arrive at the gross 
weight. The authors also consider the 
recovery of the first stage by means of 
human pilot-controlled descent with a 
rocket-motor brake and give some thought 
to the return trip from the satellite station, 
based on Dr. Sanger’s ideas, on a large- 
wing glider in a “skipping’’ trajectory, 
circling the earth several times to brake the 
velocity gradually. In a final cost esti- 
mate the authors envision a satellite sta- 
tion for experimental and research pur- 
poses with a structural weight of 180 tons 
and 330 tons of equipment. This total 


could be transported in three hundred 
supply vehicles of which every fifth would 
be a return glider. Three years would be 
required to build the satellite at a total 
cost of 500 million dollars, including the 
ground-station organization and its per- 
sonnel of 1600. 

The plan for the space station itself and 
its possibilities for use are treated in the 
fourth chapter by Herman Oberth. The 
author shows his own tentative schematic 
layout and a picture of Ross and Smith’s 
proposal, and lets his fancy run the whole 
gamut from conservative research to 
speculative conjecture. However, the 
student of astronautics, engaged in the 
tedious footwork of slow and sober pre- 
paratory techniques, may welcome an 
occasional boost of morale with such a 
quick glance into the future! 

In the fifth and final chapter by Heinz 
von Diringshofen, medical problems of 
extraterrestrial travel are classified in two 
groups—the well-known effects of high 
acceleration on the human body and other 
recent experiences related to high-speed 
flight; and the totally unknown reaction of 
man to exposure of his delicately balanced 
physical and mental equilibrium to removal 
of gravity. Only the most rudimentary 
glance at this vast unexplored field is 
given, together with an optimistic belief in 
ultimate success which is based on very 
scant evidence. 

The book closes with an enumeration of 
various societies for space research, the 
International Astronautical Federation, 
anditsmembers. Also a historical bibliog- 
raphy on space travel is added which goes 
back to Goddard’s work in 1919. 


Rerractory Harp Merats, by P. 
Schwarzkopf, R. Kieffer, and Associates, 
The Macmillan Co., N. Y., 1953, 447 
pp., 92 illus. $10. 

Reviewed by Pot Duwez 
California Institute of Technology 
Jet Propulsion Laboratory 


This book deals with the structure, 
preparation, and properties of refractory 
carbides, nitrides, silicides, and borides. 
The authors, who are leading authorities in 
the field, present a critical evaluation of 
our knowledge of these extremely impor- 
tant high-temperature materials. The book 
will be particularly welcomed by those 
interested in fundamental research. Since 
the work of Becker published in 1933 
(Hochschmelzende Hartstoffe), there has 
not been a book covering the field, and 
many of the original papers published in 
the intervening period (especially those on 
borides) have not been easy to find. The 
really thorough literature coverage pre- 
sented in this book will be greatly appreci- 
ated and will help in planning future re- 
search programs. 

The applications are dealt with in each 


chapter covering each individual material. 


In addition, however, the last part of the 
book is a comprehensive summary of the 
present position of the carbides, nitrides, 
silicides, and borides in the field of high- 
temperature materials, especially in view 
of their application to gas turbines. The 
comparison of the physical properties of 
these materials with the more convention:! 
high-temperature alloys is particularly 
useful. 

Erratum Note: The editor of this section 
regrets that the price of this book was 
incorrectly listed as $7.50 in the May-June 
issue of the JoURNAL. 


Quantum Cuemistry, by K. 8. Pitzer, 
Prentice-Hall, Inc., N. Y., 1953, 529 pp. 
$10. 

Reviewed by Davin ALTMAN 
California Institute of Technology 
Jet Propulsion Laboratory 


The book is an outgrowth of man) 
original papers by the author on various 
phases of quantum mechanics and of 
several years of teaching a course on this 
subject matter at the University of Cali- 
fornia, which this reviewer was fortunate 
enough to attend. Briefly, the book pre- 
sents the essentials of quantum theory of 
interest to chemists and treats many 
classical problems in sufficient detail to 
provide a working knowledge of the prin- 
ciples. In keeping with the author’s 
opinion that it is not necessary to employ 
advanced mathematics to gain a_ basic 
understanding of the fundamentals of 
quantum mechanics, this book has been 
written mainly for exposition of basic 
principles rather than formulation of 
mathematical techniques for detailed solu- 
tion of problems. There are 24 appendixes 
for treatment of specific topics of more 
mathematical interest. The underlying 
approach of the presentation is contained 
in the following statement by the author: 
‘““Many persons, the writer included, tend 
to think first in terms of qualitative physi- 
cal phenomena and to follow with the 
mathematics appropriate to the model so 
conceived. Such individuals should not 
discard models but rather remind them- 
selves frequently of the distinction be- 
tween the concepts established by experi- 
ment and the parts which are only figments 
of imagination.” 

The book contains thirteen chapters of 
which the first seven are devoted mainly to 
the development of the principles of the 
quantum theory, and the second six chap- 
ters to a more or less detailed discussion o! 
many practical problems of interest to the 
chemist. The first part of the book which 
treats the general theory is adequate, al- 
though it adds very little to the treat- 
ments which have already been given in 
the many other books on this subject 
The second half of the book, however, con- 
tains much valuable information on such 
topics as chemical bonds and valences, 


spectra of polyatomic molecules and statis- 
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tical thermodynamic calculations concern- 
ing them, properties of crystalline solids, 
imperfect gasses, and electrical and nuclear 
phenomena. This discussion is inter- 
spersed with numerous practical examples 
which are very instructive for indicating 
the method of treating the problem. In 
addition the book has problem sets at the 
end of each chapter which make it well- 
suited for use in a teaching course on 
quantum chemistry. 

In general, the book has a sufficient 
volume of new and interesting material 
which cannot be found in any other single 
published source but which would recom- 
mend it to the chemist or practical physi- 
cist working the field. Although the price 
of $10 is not inconsistent with the subject 
matter contained in the book, it appears 
that the quality of the binding and the 
paper could be improved at this price. 


MecHANICAL ENGINEERING THERMODY- 
namics, by D. A. Mooney, Prentice- 
Hall, Inc., N. Y., 1953, 540 pp. $9.35. 

Reviewed by R. S. Wick 
California Institute of Technology 
Jet Propulsion Laboratory 


This textbook is written for junior or 
senior engineering students and consists of 
essentially two parts. The first deals with 
general thermodynamics, and the second 
with engineering applications. Although 
the title implies that this textbook was 
written primarily for mechanical engineer- 
ing students, the first part of the text is 
sufficiently general for a first course in 
thermodynamics in other curricula. The 
first twelve chapters deal with the first 
and second laws and their applications, the 
properties of pure substances, steam, gases 
and gaseous mixtures. The thirteenth 
chapter treats process calculations for sta- 
tionary systems. These chapters could 
very well form a first-semester course in 
thermodynamics. The final twelve chap- 
ters cover the usual engineering applica- 
tions such as steam and gas cycles, tur- 
bine, compressor, refrigeration, air con- 
ditioning, internal-combustion-engine cy- 
cles, and combustion processes. 

This well-written text is designed for the 
student who wants to know not only 
“how” but “why” when it comes to engi- 
neering thermodynamics. Because of its 
clearness and conciseness, it is also recom- 
mended to those interested in a self-study 
or review course in thermodynamics. 


INTERNAL Bauuistics, F. R. W. Hunt, 
editor, Philosophical Library, Inc., 
N. Y., 1951, 311 pp. $12. 

Reviewed by A. J. ZAEHRINGER 
Thiokol Corporation 


This book, which is a collective work on 
internal ballistics and its application to 
guns, was sponsored by the Scientific Advi- 
sory Council of the British Ministry of 
Supply. Although primarily a treatise on 
gun ballistics, there is much which is 
applicable to solid rocket propellants. 
The first five chapters of the book are de- 
voted to the chemical, thermal, and ballis- 
tic properties of propellants, and two bal- 
listie equations are derived relating to the 
mode and rate of burning. Numerical 
tables are given for thermochemical com- 


putations. The gun is introduced in 
Chapters VI and VII, and the energy equa- 
tion and the equation of motion of the shot 
are deduced. 

In Chapter VIII a solution of the ballis- 
tic equations for a linear law of burning is 
given, and Chapter IX is devoted to other 
solutions of the same problem. Solutions 
for nonlinear law of burning are given in 
Chapter X. Chapter XI is devoted to a 
number of approximations which are fre- 
quently used. Chapters XII and XIII 
deal with the measurement of velocity and 
pressure, and the subject of cordite proof 
is treated in Chapter XIV. Chapter XV 
gives an outline of the application of statis- 
tical methods to cordite production and 
proof; Chapter XVI reviews the results of 
recent experimental work in America, Ger- 
many, and England. 

Appendix I deals extensively with the 
theory of leaking guns into which cate- 
gory falls the recoilless gun (which may be 
considered as a rocket and gun built to- 
gether). The book closed with a bibliog- 
raphy of theoretical works on internal 
ballistics. 

Although the greater part of the book 
deals with guns, the early sections are of 
general theoretical interest to those in the 
solid rocket propellant field. For ex- 
ample, thermochemical computations are 
treated, and together with the tables fur- 
nished, basic parameters for single- and 
double-base rocket propellants may be cal- 
culated. Also the sections on pressure 
measurements and statistical methods are 
of value to those engaged in the testing and 
analysis of solid-propellant rocket motors. 
One should not, however, expect to find 
much carry-over of quantitative informa- 
tion from one field to the other, as the 
ranges of variables involved differ by 
orders of magnitude. 


THEORY OF THE INTERIOR BALLISTICS OF 
Guns, by J. Corner, John Wiley & Sons, 
N. Y., 1950, 443 pp. $8. 

Reviewed by Jack LORELL 
California Institute of Technology 
Jet Propulsion Laboratory 


This is a textbook on theoretical interior 
ballistics of guns, written by one of the 
foremost British authorities on the sub- 
ject. The emphasis is on methods rather 
than magnitudes. Thus examples of typi- 
cal data are included, but there are no 
extensive tables of ballistic magnitudes. 
The subject matter is very definitely 
limited to interior ballistics, with little of 
any tie-in either to gun design techniques 
or to exterior ballistics. 

In theoretical interior ballistics, the 
problem of determining muzzle velocity 
and peak pressure when the characteristics 
of shot, charge, and gun are known moti- 
vates the entire subject. This problem, 
known as the central problem of interior 
ballistics, forms the theme for the book, to 
which are related discussions of topics 
ranging from ballistic methods proper to 
the chemistry of propellants. 

In Chapter I, the scope of the book is 
laid out, and a short discussion of experi- 
mental ballistics is included. Chapters II 
and III are devoted to gun propellants. 
Topics such as flame propagation in one 
dimension, bulk burning rates in guns, 

F propellant erosion, and the thermochem- 


istry of propellants are discussed in some 
detail. Chapters IV, V, and VI describe 
the various ballistic theories that are used 
in attacking the central problem of interior 
ballistics. Chapters VII and VIII involve 
the ballistics of special types of guns, for 
which the methods of the previous chap- 
ters must be modified. Of particular 
interest to the rocket engineer, perhaps, is 
the theory of the recoilless gun, which uses 
the rearward thrust of part of the gases 
from the combustion chamber released 
through a nozzle to balance the recoil. 
Chapter IX discusses hydrodynamic prob- 
lems in guns. The main difference be- 
tween the gas flow in rockets and guns is, of 
course, that in the former the flow is 
steady-state, whereas in the latter it is 
nonsteady. As a consequence, supersonic 
velocity can be reached in a straight gun 
tube, whereas a throat is required in order 
to obtain supersonic velocity in a rocket 
motor. Chapter X is a discussion of heat 
transfer in guns. In addition to the text 
there are three appendixes providing aids 
for numerical ballistic calculations, a name 
index, and a subject index. 

As a text in the theory of the interior 
ballistics of guns, this book is excellent. 
However, its use to the rocket engineer 
may be limited since the theories of rocket 
and gun ballistics, though superficially 
closely related, diverge widely when it 
comes to quantitative calculations. For 
example, in guns, pressures are of the 
order of 50,000 psi, which is perhaps two 
orders of magnitude higher than those used 
in rockets. As a result, the thermody- 
namics of the powder gases, for instance, 
are entirely different in the two cases. 


Books 

Structure of Metals. Crystallographic 
Methods, Principles, and Data (Metal- 
lurgy and Metallurgical Engineering 
Series), by C. S. Barrett, McGraw-Hill 
Book Co., Inc., N. Y., 1952 (second edi- 
tion), 661 pp. $10. 

Titanium—Its Occurrence, Chemistry 
and Technology, by J. Barksdale, Ronald 
Press Co., N. Y., 1949, 591 pp. $10. 

Introduction to the Theory of Proba- 
bility and Statistics, by N. Arley and K. R. 
Buch, John Wiley & Sons, Inc., N. Y., 
1950, 236 pp. $4. 

Small Particle Statistics, by G. Herdan 
(with a Guide to the Experimental Design 
of Particle Size Determinations, by M. L. 
Smith), Elsevier Press, Houston, Texas 
1953, 520 pp. $12. 

Detonation in Condensed Explosives, 
by J. Taylor, Oxford University Press, 
N. Y.; Oxford, at the Clarendon Press, 
1952, 196 pp. $5. 

Organic Chemistry, by R. Q. Brewster, 
Prentice-Hall, Ine., N. Y., 1953, 855 pp. 
$9.35. 

Advances in Geophysics, Vol. I, H. E. 
Landsberg, editor, Academic Press, Inc., 
N. Y., 1952, 362 pp. $7.80. 

Elements of Cartography, by A. H. 
Robinson, John Wiley & Sons, Inc., N. Y., 
1953, 254 pp. $7. 

Vision Through the Atmosphere, by 
W. E. K. Middleton, University of To- 
ronto Press, Toronto, Canada, 1952, 250 
pp. $8.50. 

The Comets and Their Origin, by R. A. 
Lyttleton, Cambridge University Press, 
N. Y., 1953, 173 pp. $5. 
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Jet Propulsion Engines 


A Note on the Performance of Ducted 
Fans, by Bryan Thwaites, Aeron. Q., vol. 
4, part IT, Feb. 1953, pp. 179-185. 

Theory and Practice of Gas Turbine 
Power Plants for Helicopters, by Igor B. 
Bensen, Amer. Helicopter, vol. 30, May 
1953, pp. 6-10, 13. 

Hydraulic Regulation of Feeding of 
Turbomachines, by Basile Demtchenko 
(in French), France, Ministére de I’ Air. 
Pub. Sci. Tech., no. 277, 1953, 182 pp. 

Analysis of the Performance Charac- 
teristics of Turbojet Engines as a Function 
of Operating Conditions, by Ludovico 
Pascuci (in Italian), Aerotecnica, vol. 33, 
Feb. 1953, pp. 94-99. 

Possibilities of Turbojets and Ramjets 
at Supersonic Speeds, by Gaspare Santan- 
gelo (in Italian), Aerotecnica, vol. 33, Feb. 
1953, pp. 114-117. 

Research and Development Progress, 
Propulsion, by Abe Silverstein, Aviation 
Age, vol. 19, June 1953, pp. 182-189. 

Flying Tests on Supersonic Ram Jet 
Tubes at Low Subsonic Speed, by Irene 
Sanger-Bredt (in German), VDI-For- 
schungsheft, no. 437, 1953, pp. 40-60. . 

Physical Fundamentals of Jet Propul- 
sion, by Eugen Sanger (in German), 
VDI-Forschungsheft, no. 437, 1953, pp. 
5-25, 46 references. 

Specialization in Turbojet Engine De- 
sign, by Bill Krase and Bob Neitzel, Aero 
Digest, vol. 66, June 1953, pp. 74, 76, 
78-80, 82, 84, 86, 88-95. 

Research on Aerodynamic Noise from 
Jets and Associated Problems, by E. J. 
Richards, J. Roy. Aeron. Soc., vol. 57 
May 1953, pp. 318-342, 25 references. 

Some NACA Research on Centrifugal 
Compressors, by I. A. Johnsen, and 
Ambrose Ginsburg, Trans. ASME, vol. 75, 
July 1953, pp. 805-817. 

Periodic Flow Regenerator—a Summary 
of Design Theory, by J. E. Coppage and 
A. L. London, Trans. ASME, vol. 75, July 
1953, pp. 779-787. 

The Rotary Regenerative Air Preheater 
for Gas Turbines, by A. T. Bowden and 
W. Hryniszak, Trans. ASME, vol. 75, July 
1953, pp. 767-777. 

Ramjet Possibilities, Flight, vol. 114, 
July 3, 1953, pp. 7-12. 

Aircraft Powerplant Seals, by Norman F. 
Rooke, SAE VJ/., vol.61, July 1953, pp. 23-29. 

Exhaust Augmented Engine Cooling, 
Shell Aviation News, no. 179, May 1953, 
pp. 19-20. 

Theory and Practice of Gas Turbine 
Power Plants for Helicopters, part II, by 
Igor B. Bensen, Amer. Helicopter, vol. 31, 
June 1953, pp. 10-13, 16. 

Tailpipe ‘“Teeth’’ Cut Jet Engine Noise, 
eg Week, vol. 59, July 6, 1953, pp. 
28-29. 


Heat Transfer and Fluid 
Flow 


On the Noise Emanating from a Two- 
Dimensional Jet Above the Critical Pres- 


sure, by Alan Powell, Aeron. Q., vol. 4, 
part II, Feb. 1953, pp. 103-122, 12 refer- 
ences. 

Heat Transfer to Bodies in a High-Speed 
Rarefied-Gas Stream, by Jackson R. 
Stalder, Glen Goodwin, and Marcus O. 
Creagor, NACA Rep., no. 1093 (supersedes 
NACA TN 2438), 1952, 9 pp. 

Physical Phenomena in Gas Jets. 
Final Report, by Maryland University, 
Institute for Fluid Dynamics and Applied 
Mechanics, Oct. 1952, 53 pp. Contents: 
Optical Study of Compressible Jets, by 
D. Bershader and B. Cary, pp. 1—20; On 
Supersonic Flow of a Two-Dimensional 
Jet in Uniform Stream, by 8. I. Pai, pp. 
21-32 (J. Aero. Sci., vol. 19, Jan. 1952, pp. 
61-65); Problems in the Study of Turbu- 
lence, by R. Betchov, pp. 44-53; Axially 
Symmetrical Jet “sys of a Compressible 
Fluid, by S. I. Pai, 33-43 (Q. Appl. 
Math., vol. 10, July 1952, pp. 141-148). 

Shape and Thickness of Shock Fronts in 
Argon, Hydrogen, Nitrogen, and Oxygen, 
by E. F. Greene, and D. F. Hornig, J. 

Chem. Phys., vol. 21, April 1953, pp. 617- 
623. 

Flow Changes in Gases in Which Mass 
and Impulse are Conserved, by L. G. 
Dawson, Aeron. Q., vol. 4, part II, Feb. 
1953, pp. 193-204. 

A Note on the Theory of the Constant- 
Area Mixing of Compressible Flows as 
Applied to High-Speed Wind Tunnel De- 
sign, by F. G. Irving, Aeron. Q., vol. 4, 
part II, Feb. 1953, pp. 164-178. 

On the Stability of the Flow of Gas out 
of a Compressor, by J. M. Stephenson, 
J. Roy. Aeron. Soc., vol. 57, May 1953, pp. 
345-346. 

Schlieren Analysis of Supersonic Three- 
Dimensional Flows, by Claude Picard (in 
French), Recherche Aeron., no. 32, March- 
April 1953, pp. 15-19. 

Calculation of a Nozzle and of a Profile 
in the Hodograph Approximation, by 
A. Gilles (in French), Recherche Aéron., no. 
32, March-April 1953, pp. 1—10. 

Calculation of Cascade Flows by a 
Field of a Series of Whirls, by A. Betz (in 
German), Z. Angew. Math. Mech., vol. 33, 
no. 4, April 1953, pp. 113-116. 

A Visualization Study of Secondary 
Flows in Cascades, by Arthur G. Hanson, 
NACA TN no. 2947 (formerly RM no. 
E52F19), May 1953, 93 pp. 

Measurement of Heat Transfer and Skin 
Friction at Supersonic Speeds; Prelimi- 
nary Results of Measurements on Flat 
Plate at Mach Number of 2.5, by J. E. 
Johnson and R. J. Monoghan, Gt. Brit. 
Aeron. Res. Council, Current Paper, no. 59 
(formerly ARC TR 12,483, RAE TN Aero 
1994), 1951, 18 pp., 17 figs. 

One-Dimensional High-Speed Flows, by 
J. Kestin and 8. K. Zaremba, Aircr. Engng., 
vol. 25, June 1953, pp. 172-175, 179; 19 
references. 

Flows Through Nozzles and Related 
Problems of Cylindrical and Spherical 
Waves, by Yu Why Chen, Comm. Pure 
Appl. Math., vol. 6, May 1953, pp. 179- 
228; 22 references. 


Pulsed Surface Heating of a Semi-Finite 
Solid, by J. C. Jaeger, Q. Appl. Math., vol. 
11, April 1953, pp. 132-137. 

Subsonic Flow of Air Through a Single- 
Stage and Seven-Stage Compressor, by 
Chung-Hua Wu, NACA TN, no. 2961, June 
1953, 32 pp. 

Experimental Investigation of the Mis- 
ing Loss Behind the Trailing Edge of a 
Cascade of Three 90° Supersonic Turning 
Passages, by Luke L. Luccini, NACA RM 
L50F2la, Aug. 1950 (declassified from 
Confidential 4/13/53), 31 pp. 

A Method for Stabilizing Shock Waves 
in Channel Flow by Means of a Surge 
Chamber, by Stanford E. Neice, NACA 
TN no. 2694, June 1953, 46 pp. 

Theoretical and Experimental Investig:- 
tion of the Flow in a Duct of Varying 
Cross-Section, with Particular Application 
to the Design of Ducts for Free Flight 
Ground-Launched Model Tests, by C. H. 
E. Warren, R. E. Dudley, and P. J. 
Herbert, Gt. Brit. Aeron. Res. Council, 
Current Paper no. 60 (formerly ARC TR 
no. 13,726, RAE TN Aero 2064), 1951, 25 
pp., 12 figs. 

Radial-Axial Heat Flow in Regions 
Bounded Internally by Circular Cylinders, 
by J. H. Blackwell, Can. J. Phys., vol. 31, 
May 1953, pp. 472-479, 

E ‘ficiency and Drag Axial-Flow 
Compressor Stage, by J. M. Stephenson, 
Airer. Engng., vol. 25, co 1953, pp. 
158-160. 

Preliminary Analysis of Axial-Flow 
Compressors Having Supersonic Velocity 
at the Entrance of the Stator, by Antonio 
Ferri, NACA RM no. L9GO6, Sept. 1949 
(declemied from Confidential 4/13/53), 
36 p 

Fheoretical and Experimental Analysis 
of Low-Drag Supersonic Inlets Having a 
Circular Cross Section and a Central 
Body at Mach Numbers of 3.30, 2.75, and 
2.45, by Antonio Ferri and _ Louis 
Nucci, NACA RM L8H13, Nov. 1948 (de- 
classified from Confidential 4/13/53), 89 


pp. 
Combustion 


Factors Controlling the Combustion of 
Zirconium Powders, by Holger C. Ander- 
son, Lawrence H. Belz, J. Electrochem. 
Soc., vol. 100, May 1953, pp. 240-249. 

A Theoretical Model of a Gaseous Com- 
bustion Wave Governed by a First-Order 
Reaction, by Raymond Friedman and 
Edward Burke, J. Chem. Phys., vol. 21, 
April 1953, pp. 710-714. 

Aerothermodynamics and Combustion 
Theory, by Theodore von Karman, Aero- 
tecnica, vol. 33, Feb. 1953, pp. 80-86. 

Radical Reactions of Nitric Acid in 
Flames, by G. K. Adams, W. G. Parker, 
and H. G. Wolfhard, Discussions Faraday 
Soc., no. 14, 1953, pp. 97-103. 

Free Radicals in Explosions Studied by 
Flash Photolysis, by R. G. W. Norrish, 
— Faraday Soc., no. 14, 1953, pp. 
16-22 

The Mode of Action of Lead Tetraethy! 


Eprtor’s Note: This collection of references is not intended to be comprehensive, but is rather a selection of the most significant and 


stimulating papers which have come to the attention of the contributors. 
literature is unavailable because of security restrictions. 
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as an Inhibitor of Combustion Processes, 
by G. H. N. Chamberlain, D. E. Hoare, 
and A. D. Walsh, Discussions Faraday 
Soc., no. 14, 1953, pp. 89-96. 

A Slot Burner Method for Studying 
Combustion Wave Instability, by G. H. 
Markstein and L. H. Somers, J. Chem. 
Phys., vol. 21, May 1953, p. 941. 

Burning Velocity Measurement, by 
Dorothy M. Simon and Edgar L. Wong, 
J. Chem. Phys., vol. 21, May 1953, p. 936. 

Reaction Processes Leading to Spontane- 
ous Ignition of Hydrocarbons, by Charles 
E. Frank and Angus U. Blackham, NACA 
TN no. 2958, June 1953, 27 pp. 

Flame Propagation Rates the 
Chemical Nature of the Attachment Sur- 
face, by L. Lapidus, J. B. Rosen, and 
R. N. Wilhelm, Princeton Univ., Forrestal 
Res. Center, Chemical Kinetics Project, 
TR no. 1, April 1953, 23 pp., 12 figs. 

The Thermodynamics of Combustion 
Gases: Temperatures of Methane-Air, 
Propane-Air, and Ethylene-Air Flames, by 
Robert W. Smith, Jr., John Manton, and 
Stuart R. Brinkley, Jr., Bureau of Mines. 
Report of Investigations 4983, June 1953, 21 


). 

Studies on the Spontaneous Ignition of 
Fuels Injected into a Hot Air Stream. V. 
Ignition Delay Measurements of Hydro- 
carbons, by B. P. Mullins, Fuel, vol. 22, 
July 1953, pp. 363-379. 

Studies on the Spontaneous Ignition of 
Fuels Injected into a Hot Air Stream. 
IV. Ignition Delay Measurements on 
Some Gaseous Fuels at Atmospheric and 
Reduced Static Pressures, by B. P. 
Mullins, Fuel, vol. 32, July 1953, pp. 343- 
362. 

Spontaneous Ignition of Fuels Injected 
into a Hot Air Stream. III. Effect of 
Chemical Factors upon the Ignition Delay 


of Kerosene-Air Mixtures, by B. P. 
-— Fuel, vol. 32, July 1953, pp. 327— 


Analytical (Power Series) Solutions to 
the Equations of Flame Propagation, by 
J. O. Hirschfelder and Edwin S. Campbell. 
Wisconsin Univ. Naval Res. Lab. CM-784, 
May 1953, 117 pp. 


Fuels, Propellants and 
Materials 


Heat Capacities, Latent Heats and 
Entropies of Fluorine from 13° to 85°K: 
Heats of Transition, Fusion, Vaporization, 
and Vapor Pressures of the Liquid, by Jih- 
Heng Hu, David White, and Herrick L. 
Johnston, Ohio State Univ. Res. Founda- 
tion, TR no. 283-23, April 1953, 11 pp. 

An Investigation of Boron Carbide, by 
Frank W. Glaser, David Moskowitz, and 
Benjamin Post, J. Appl. Phys., vol. 24, 
June 1953, pp. 731-733. 

Determination of Excessive Sludge 
Formation in Mixed Acid Storage Drums, 
by Arthur E. Lenehan, Naval Air Rocket 
Test Sta. Rep. no. 29, Feb. 1953, 8 pp. 

Optical Absorbance of the Ternary Sys- 
tem of Nitric Acid-Nitrogen Dioxide- 
Water, by Scott Lynn, David M. Mason, 
and B. H. Sage, Calif. Inst. Tech. Jet Pro- 
pulsion Lab. Progress Rep. no. 20-187, 
Jan. 1953, 22 pp. 

Storage and Handling of Nitrogen 
Oxides, by F. T. Selleck, H. H. Reamer, 
and B. H. Sage, Indust. Engng. Chem., vol. 
45, April 1953, pp. 814-819. 

Studies on Hydrazine. Technical Re- 
port No.,1. The Photolysis of Ammonia 
at 1849A in a Flow System, by C. C. 
McDonald, A. Kahn, and H. E. Gunning, 
Illinois Inst. Tech. Department of Chem- 


istry, Feb. 1953 (declassified from Re- 
stricted June 23, 1953), 20 pp. 

Effect of Prestraining on Recrystalliza- 
tion Temperature and Mechanical Proper- 
ties of Commercial, Sintered, Wrought 
Molybdenum, by Kenneth C. Dike and 
Roger A. Long, NACA TN 2973, July 
1953, 25 pp. 

Investigation of Effects of Grain Size 
upon Engine Life of Cast AMS 5385 Gas 
Turbine Blades, by Charles A. Hoffman 
and Charles A. Gyorgak, NACA RM no. 
£53D06, July 1953, 21 pp. 

Behavior of Materials Under Conditions 
of Thermal Stress, by Samuel S. Manson, 
NACA TN no. 2933, July 1953, 105 pp. 


Physical-Chemical 
Topics 


Analytical Procedures for Rocket Pro- 
pellants. VII. Mixed Acid, by John D. 
Clark, Naval Air Rocket Test Sta. Rep. no 
24, Sept. 1952, 18 pp. 

Kinetics of Two Exchange Reactions 
Involving Diborane, by P. Calvin May- 
bury and W. S. Koski, J. Chem. Phys., 
vol. 21, April 1953, pp. 742-747. 

The Thermal Decomposition of Nitrous 
Oxide, by Lewis Friedman and Jacob 
Bigeleisen, J. Amer. Chem. Soc., vol. 75, 
Mav 5, 1953, pp. 2215-2217. 

The Reactivity of Free Radicals, Dis- 
cussions Faraday Soc. no. 14, 1953, 256 pp. 

Free Radicals by Mass Spectrometry. 
Part II. The Thermal Decomposition of 
Ethylene Oxide, Propylene Oxide, Di- 
methyl Ether, and Dioxane, by F. P. 
Lossing, K. U. Ingold, and A. W. Tickner, 
Discussions Faraday Soc. no. 14, 1953, pp. 
34-44. 

Studies of Free Radical Reactivity by 


Reynolds Electrical 


Engineering Co., Inc. 


Electrical and Construction Engineers 


week 


EL PASO 
HOUSTON 


LAS VEGAS 
AN ORGANIZATION OF ELECTRICAL 
ENGINEERS TRAINED TO THE INTRICATE 
NEEDS OF THE CONSTRUCTION INDUSTRY 


We are the leading manufacturer of precision 
orifices in the world. We manufacture Spin- 
nerettes used to produce synthetic fibers. Each 
Spinnerette requires from 2 to 25,000 preci- 
sion orifices with tolerances of plus or minus 
0.00004 inches. Let this background help your 
production of precision fuel injector nozzles. 
Write or call Department O for information. 


BAKER & CO., INC. 


113 ASTOR STREET, NEWARK 5, NEW JERSEY 
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the Methods of Flash Photolysis; the 
Photochemical Reaction Between Chlorine 
and Oxygen, by George Porter and Frank- 
lin J. Wright, Discussions Faraday Soc. no. 
14, 1953, pp. 23-34. 

The Absorption — of Free NH: 
Radicals, by G. Herzberg and D. A. 
Ramsay, ” Discussions Faraday Soc. no. 14, 
1953, pp. 11-16. 

Infra-Red Spectrum of Nitric and 
Deuteronitric Acid Vapour. Completion 
of the Identification of the Fundamental 
Frequencies. Entropy of Nitric Acid. 
Barrier Resisting Rotation of the Hy- 
droxy] Group, by Henry 
Ingold, and i. G G. Poole, . Chem. Soc., 
Nov. 1952, pp. 

Molecular Sound Velocity and Molecu- 
lar Compressibility of Liquid Mixtures, 
by Otohiko Nomoto, J. Chem. Phys., vol. 
21, May 1953, pp. 950-951. 

Kinetics of Reaction Between Ammonia 
and Oxygen in a Quartz Vessel, by Henry 
Wise and Maurice F. Frech, J. Chem. 
Phys., vol. 21, May 1953, pp. 948-949. 

The Normal Boiling Points of Oxygen 
on the Thermodynamic Scale, by J. G. 
Aston and G. W. Moessen, J. Chem. Phys., 
vol. 21, May 1953, p. 948, 

Monometric Determination of the Den- 
sity of Liquid Ozone, by Callaway Brown 
and Karl D. Franson, J. Chem. Phys., vol. 
21, May 1953, pp. 917-919. 

Applied Chemical Kinetics. VI. De- 
termination of Kinetic Constants. c. 
Parallel reactions, by J. C. Jungers and 
A. Giraud (in French), Rev. Inst. frangais 
Pet. Am. Combust. Liquides, vol. 8, April 
1953, pp. 152-168. 

Contribution to the Theory of Fast 
Reaction Rates, by R. DeVogelaere and 
M. Boudart, Princeton Univ. James Forres- 
tal Res. Center. Chemical Kinetics Proj- 
ect, TR no. 3, July 1953, 40 pp.; 6 figs. 

Kinetics and Mechanisms of the 
‘‘Water-Gas” Reactions, by W. M. 
Graven and F. J. Long, Princeton Univ. 
James Forrestal Res. Center. Chemical 
Kinetics Project, TR no. 2, May 1953, 27 
pp., 10 figs 

Frictional Behavior of Polyethylene, 
Polytetrafluoroethylene and Halogenated 
Derivatives, by R. C. Bowers, W. D. 
Clinton, and W. A. Zisman. Naval Res. 
Lab. Rep. no. 4167, May 1953, 21 pp. 

Transport Properties of Gases Obeying a 
Modified Buckingham (Exp-Six) Poten- 
tial, by Edward A. Mason, Wisconsin 
Univ. Naval Res. Lab. Rep. ONR-1, June 
1953, 64 pp. 


Instrumentation and Ex- 
perimental Techniques 


Meter Measures Fuel Flow in Pounds 
(Gaveo Corp., Mass Flowmeter), Aviation 
Week, vol. 58, April 6, 1953, pp. 62, 64, 66. 

Temperature Measurement; the Use 
of Heat-Sensitive Colour-Changing Pig- 
ments, Aircr. Prod., vol. 15, May 1953, p. 
184. 

Apparatus for the Measurement of the 
Particle Content of Air (in High Speed 
Wind Tunnels), by P. Durrenberger (in 
French), Recherche Aéron., no. 32, March- 
April 1953, p. 14. 

Microsecond Photography of Rockets in 
Flight, by E. Barkofsky, R. Hopkins, and 
8. Dorsey, Electronics, vol. 26, June 1953, 
pp. 142-147. 

Frequency Modulated Magnetic Tape 
Recording and Playback Instrumentation 
System, by Joseph Petes, Naval Ordn. 
Lab. Rep. no. 1167, pp. 75-136; Navord 
Rep. no. 2713, Feb. 1953. 

Development of an Inductance Type 
Accelerometer for Use in a Frequency 
Modulated Recording System, by Benja- 


min Sussholz, Naval Ordn. Lab. Rep. no. 
1167, pp. 21-67; Navord Rep. no. 2712, 
Feb. 1953. 

Three-Dimensional Liquid Analog for 
the Determination of Temperature Dis- 
tribution, by G. A. Sterbatzel and J. L. 
Beal, Project nie TM no. CAL-44, Oct. 
1952, 22 pp., 14 figs. 

The eopler Schlieren Apparatus, by 
D. W. Holder and R. J. North, Gt. Brit. 
Aeron. Res. Council, Reps. and "Mem. no. 
2780 (formerly ARC TR no. 13,068), 1953, 
13 pp., 6 figs. 

A Variable High Frequency Light Syn- 
chronized with a High-Speed Motion-Pic- 
ture Camera to Provide Very Short Expo- 
sure Times, by Walter Frank Lindsey and 
Joseph Burlock, NACA TN no. 2949, May 
1953, 17 pp. 

Possibility of Simulation of a Flying 
Structure With Free Control Surfaces by a 
Servomechanism, by Coupry and 
R. Valid (in Frene h), Recherche Aéron., 
no. 32, March-April 1953, pp. 29-34. 

A Small Pirani Gage for Measurements 
of Non-Steady Flow Pressures, by Miro- 
slav John Pilney, NACA TN no. 2946, 
June 1953, 36 pp. 

Experiments with Static Tubes in a 
ae Air-Stream. Parts I and II, 

by D. W. Holder, R. J. North, and A. 
Chinneck, Gt. Brit. Aeron. Res. Council 
Reps. and Mem. no. 2782 (formerly AR 
TR no. 13,269 and 13,268), 1953, 14 pp., 
13 figs. 

Measurement and Registration of Short 
Time Intervals, by Jan Groeneveld (in 
German), Z. Phys., vol. 134, April 17, 
1953, pp. 645-647. 

An Instrument for the Direct Measure- 
ment of Intense Thermal Radiation, by 
Robert Gardon, Rev. Sci. Instrum., vol. 24. 
May 1953, pp. 366-370. 

A Method of Control of a Predeter- 
mined Flow Rate, by H. H. Reamer and 
B. H. Sage, Rev. Sci. Instrum., vol. 24, no. 
5, May 1953, pp. 362-366. 

A High Temperature Thermo-Regulator, 
by F. Rosebury, Rev. Sci. Instrum., vol. 
24, no. 5, May 1953, pp. 398-399. 

Continuous Recording of Pressures for 
Supersonic Wind-Tunne! Calibration, by 
L. P. Gieseler, Navord Rep. no. 2744, 
Jan. 1953, 14 pp., 7 figs. 

Measurement by Times Spark Shadow- 
graphs of Shock Velocities in the Shock 
Tube, by Herbert L. Hoover, Lehigh Univ. 
Inst. of Res. TR no. 3, July 1953, 24 pp. 

A Modern Jet Engine Test Plant (Stand- 
ard Motors, Coventry), Aeroplane, vol. 
84, May 15, 1953, p. 643. 


Terrestrial Flight, Ballis- 
tics, and Vehicle Design 


Ghost Jets in the Sky (United Air Lines 
research on simulated jet flight service), 
Aviation Age, vol. 19, May 1953, pp. 6-9. 

The Air Force and Guided Missiles, by 
E. P. Mechling, Ordnance, vol. 37, March- 
April 1953, pp. 789-790. 

The Case for Pod-Mounted Jet Engines, 
by David A. Anderton, Aviation Week, vol. 
58, May 18, 1953, pp. 29-30, 32-34, 37-38. 

The Case for the Underslung Nacelle on 
the Jet Transport, by B. T. Salmon, Avia- 
tion Age, vol. 19, May 1953, pp. 84-91. 

Properties and Tables of Generalized 
Rocket Functions for Use in the Theory of 
Rockets with a Constant Slow Spin, by 
J. Barkley Rosser and R. J. Walker, 
Cornell Univ., 1953, 111 pp. 

An Experimental Solution to the 
Lagrange Ballistic Problem, by A. E. 
Seigel, Navord Rep. 2693 (Aeroballistic 
Research Rep. 80), Sept. 1953, 9 pp. 

Guided Missiles—Comments on Surface 


Based Antiaircraft Missiles, by Henry H, 
Porter, Aero. Engng. Rev., vol. 12, July 
1953, pp. 24-29. 

Parachute Recovery of Ryan “Firebee” 
Jet Propelled Target Aircraft, Shell Avia- 
tion News, no. 179, May 1958, pp. 21. 

More about Jindivik, Aeroplane, vol. 84, 
June 1953, pp. 790. 

Recent ‘Developments in Rockets and 
Guided Missiles in the United States, by 
Norman J. Bowman, J. Space Flight, vol. 
5, May 1953, pp. 1-9. 

8. Turboprop Transport Future, by 
C. R. Smith, Aero Digest, vol. 66, a 
1953, pp. 17-19. 


Space Flight 


The Food and Atmosphere Control 
Problem on Space Vessels. Part I. 
Chemical Purification of Air, by Norman 
J. Bowman, J. Brit. Interplan. Soc., vol. 
12, May 1953, pp. 118-123. 

Optimum Conditions for Multi-Stage 
Rockets Rising Vertically in a Gravita- 
tional Field, by Irene Sanger-Bredt (in 
German), VDI-Forschungsheft, no. 437, 
1953, pp. 26-39. 


Astrophysics, Aerophysics, 
and Atomic Physics 


Table of Isotopes, by J. M. Hollander, 
I. Perlman, and G. T. Seaborg, Revs. Mod. 
Phys., vol. 25, April 1953, pp. 469-651. 

The Energy Levels and the Structure of 
Light Nuclei, by D. R. Inglis, Revs. Mod. 
Phys., vol. 25, April 1953, pp. 390-450. 

Light Pulses from the Night Sky, by 
J. V. Jelley and W. Galbraith, Phil. Mag., 
vol. 44, June 1953, pp. 619-622. 

Cloud-Chamber Observations of the 
Heavy Charged Unstable Particles in 
Cosmic Rays, by H. 8. Bridge, C. Peyron, 
B. Rossi, and R. Safford, Phys. Rev., vol. 
90, June 1, 1953, pp. 921-933. 

Ratio of Neutral to Charged Particles in 
the Nuclear Interacting Component of 
Cosmic Rays, by K. Greisen and W. D. 
Walker, Phys. Rev., vol., 90, June 1, 1953, 
pp. 915-920 

Cosmic Radiation Intensity-Time Varia- 
tions and their Origin. I. Neutron 
Intensity Variation Method and Mete- 
orological Factors, by J. A. Simpson, 
W. Fonger, and 8S. B. Treiman, Phys. Rev., 
vol. 90, June 1, 1953, pp. 934-950. 

The Instability of a Layer of Fluid 
Heated Below and Subject to Coriolis 
Forces, by S. Chandrasekhar, Proc. Roy. 
Soc., vol. 217A, May 7, 1953, pp. 306-327. 

The Coriolis Force of a Sphere in a 
Flow, by Friederich Bayer-Helms and 
Heinz Pick (in German), Z. Phys., vol. 
134, April 17, 1952, pp. 582-595 

Fluctuations and Latitude “Effect of 
Cosmic Rays at High Altitudes and Lati- 
tudes, by H. V. Neher, V. Z. Peterson, and 
E. A. Stern, Phys. Rev., vol. 90, May 15, 
1953, pp. 655-674. 

The Penetrating Component of Cosmic 
Radiation in the Upper Atmosphere, by 
J. D. Pullar and E. G. Dymond, Phil. 
Mag., vol. 44, June 1953, pp. 565-577. 


General Topics 


From Dirigibles to Missiles, by G. A. 
Crocco (in Italian), Aerotecnica, vol. 33, 
Feb. 15, 1953, pp. 6-9. 

The Birth of an Engine; British Prac- 
tice in Aviation Engine Procurement and 
Development, by F. R. Banks, Aero. 
Engng. Rev., vol. 12, June 1953, pp. 31-43. 

Aeronautical Acoustics, Flight, vol. 63, 
June 12, 1953, pp. 755-758. 
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Wright Brothers pioneers in tas _— in 190 


The Institute of Aero- 


ed Marquardt engineers . 

to prelese the vegies . powered craft in 1903, pioneered by the Wrig 
Brothers, to ramjet speeds of mach 4. Marquarc 


of the original Wright 
Brothers engine for per- vhas been the pioneer in the research and develo: 


manent display at their ment of ramjets, afterburners, air turbine. 
Durand Aeronautical Museum. 
“sory drives and thrust-controls of adv 
design. These units are @urrently in produ 
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As an Oxidant 


NITROGEN 
TETROXIDE 


offers outstanding advantages ci 
to of rocket motors } 


a _ SP ECIFIC IMPULSE: Nitrogen Tetroxide exceeds many other well- 
& | ape known oxidizers in pounds of thrust developed per pound 
Pate aie of fuel consumed per second. 


EASY T0 HANDLE: Nitrogen Tetroxide may be shipped, piped and 
stored in ordinary carbon steel equipment. It possesses 
high chemical stability, high density, low freezing point, 
reasonably low vapor pressure. 


heed 


ae Nitrogen wees “s Xide is available at tow cost in 125-pound I.C.C. approved 


bic steel cylinders and 1-ton containers. 

Address your inquiry to the Product Development Department 


ALLIED a & DYE CORPORATION 


age bk 40 RECTOR STREET, NEW YORK 6, N. Y. 
Technical service and development on Nitrogen Tetroxide—formerly handled by the Product Development Department, 
Solvay Process Division—are now handled by Nitrogen Division, Allied Chemical & Dye Corporation. 
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LIQUID OXYGEN 
GENERATORS 


An air-transportable 

plant for the separa- 
4 tion of high-purity oxy- 
44 be gen from the air. 


ih 


LIQUEFIED-GAS 
STORAGE CONTAINERS 
PROTOTYPE 
liquefied gases such as 
liquid hydrogen or oxy- 
DEVELOPMENT 


Containers for the stor- 


—_ _ Arthur D. Little, Inc. has blended 
its scientific and engineering skills in 
the Mechanical Division to provide industry 


LIQUID OXYGEN with a unique service... 
PUMPS 
» NY , A compact, noncontam- Scientists in the fields of ay 
inating unit for supply- chemistry, physics, metallurgy, mathematics, 
ing high-pressure oxy- biology, electronics and technical economics 


gen gas from low-pres- 
sure liquid oxygen stor- 
age containers. 


regularly engaged in diversified 
research and development projects 
join with our engineers in prototype 
development of equipment requiring 
a high level of engineering skills. 


Our staff is experienced in interpreting the 
ideas of industry and following 
through with the ee of 
HEAT EXCHANGERS specialized equipment. 
A gas-to-gas heat ex- 
changer with excep- 
tional characteristics, 


HELIUM oh 
REFRIGERATORS 

\\ A practical system for 
the prevention of evap- 
oration loss in stored 
liquefied gases. 


WRITE FOR BROCHURE RJ28-5 


MECHANICAL DIVISION 
AIR COOLERS 


A unit which provides Arthur D. Little, Juc. 


cooled, compressed air 30 MEMORIAL DRIVE — CAMBRIDGE 42, MASS. 
which is free of dirt, 


oil, or entrained water. Offices in New York, St. Louis, and Mexico City 


CREATIVE TECHNOLOGY SINCE 1886 


THERMODYNAMICS © HEAT TRANSFER ¢ REFRIGERATION TO MINUS 456°F » GAS LIQUEFACTION 
e VACUUM ENGINEERING + ELECTROMAGNETISM + MECHANICAL DESIGN + VIBRATION 
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parts like these! 


MODEL A-7 


CAPACITY: 20 TONS 


Put the Model A-7 to work forming 


Incorporating the latest engineer- 
ing developments, this new Model 
A-7 Hufford provides semi-auto- 
matic operation in a low tonnage 
machine. It insures considerably 
higher output than older, manual 
models of a similar tonnage range 
and offers, in addition, increased 
forming versatility with independ- 
ent arm actuation. 


FEATURES 


SEMI-AUTOMATIC OPERATION-—Stretch pres- 
sures, wrapping speed and arm positions are 
pre-set for any given workpiece. These factors 
are uniformly duplicated on each work cycle. 
GREATER UNIFORMITY OF PARTS—with semi- 
automatic operation all human variations, nor- 
mally affecting part uniformity, are eliminated. 


NO SPECIAL OPERATING SKILL—Once machine 
functions are properly established, no special 
skill is required to produce excellent parts. 
FEWER LOSSES—Work breakage from over- 
stretching is virtually eliminated on production 
runs, 
INDEPENDENT ARM ACTUATION—Each arm is 
under independent control, simplifying setups 
on non-symmetrical dies and adding to form- 
ing versatility. 
The Model A-7 is the ideal machine 

for modernizing your stretch-press 
department. Write for quotation! 


FORD 


4h Yorks D 


GUNDO, CALIFORNIA 
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CURRAN ENGINEERING CO. 


Manufacturer of 


MECHANICAL COMPONENTS 
from 
METALS, CERAMICS, AND PHENOLICS 


Consultants and Specialists of 
ROCKET IGNITER ASSEMBLIES 
and 
LONGITUDINAL SHAPED CHARGE 
CUTTERS 


A 


“CENGO" Process for 
HIGH ‘TEMPERATURE- HIGH DIELECTRIC 
INSULATING OF METALLIC ASSENBEES 


4423 W. Jefferson Blvd. 


Los Angeles 16 California. 


as 


CONTROL 


rocke 


SOUND 


INDUSTRIAL 


CONTROL 


ing testing of the 


how” 
experience — can 
instant service, 

We welcome the 
problem. For full 


SOUND CONTROL 


INDUSTRIAL 


an electronic instrument for 
the measurement of static dnd 


dynamic pressures 


Versatile, positive, dependable — 
under the most adverse conditions 


m The water-cooled capacitance pick-up permits 


sures at extremely high combustion temperatures. Remote location of the 
pressure indicator is provided by cable connection, with lengths up to 100 feet 
performing satisfactorily. The pressure indicator can be installed with the dia- 
phragm flush with the wall of the pressure chamber, eliminating surge effects. 


INDUSTRIAL 


<4 silenced by 
ISC mufflers 


@ In the laboratory, in the test cell, on the 
airfield, hundreds of INDUSTRIAL SOUND 
installations are subduing the © 
noise, heat and gas velocities generated dur- > 


Whatever your noise problem, ISC’s skilled 
engineering, design, and installation 
— gained through years of practical 


WRITE ISC TODAY 
Foreign Licensees: 
Cementation (Muffeiite) Limited, London 
Les Travaux Souterrains, Paris : 


Industrial Sound Control, Inc. 


45 Granby Street, Hartford, Conn. 
2119 So. Sepulveda Bivd., Los Angeles, Calif. 


SOUND O 


‘ 


big jets. 


“know 


4 
help you, and is at your 


challenge of the unusual © 
information 


qNnos 


+ > 


* 
oO UMD: CONT ® 


the measurement of pres- 


Write or telephone for information 


PHOTOCON researcH 


421 N. Foothill Blvd., Pasadena 8, Calif. 


Phone SYcamore 2-4131 


the rocket industry . . 
combustion pressures 


instrument of 
. for 


DYNA-GAGE DG-101 


Output: + 10 volts at 


25,000 ohms. 


Power requirements: 
PS-102. 

Size: 

83%” x 1114” x 


CAPACITANCE PICKUP 


TYPE 


wA TER £0 PRE SSYRE 
Rocker Gace 


WATER PASSAGE 
CONTACT STUQ 
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multiply your 
ao engineering manpower! 


EQUATION ‘SOLVING #¢¢¢« 


SYSTEM SIMULATION, 229286 


. SOLVE IGH- ORDER EQU 
complete solution time 

equations up to the 20th orde 
to two hours with the EASE com pind 


2. REDUCE RESEARCH AND DEVELOPMENT TIME... 
functioning of an entire dynamic system 
can quickly be simulated on the EASE. 
Effects of changing loads, forces, condi- 
tions and other variables can be deter- 
mined rapidly and accurately—at a frac- 
tion of the time and cost of operational 
testing. 


3. SPEED PRODUCTION TESTING... coupled with 
voltage transducers, the EASE computer 
can be channeled to physical systems, sub- 
systems or components. By substituting 
for, and simulating, associated equipment 
it provides substantial savings in test time. 


LEAKPROOF? 


CIRCLE-SEAL provides —, 
dead tight sealing at 
pressures @ as low as 
1” HO « as high as 
proof pressure. 

Dead tight means NO 
bubbles, NO dripping. 


EASE 30-channel linear computer 


INVESTIGATE THESE IMPORTANT EASE FEATU 


200 SERIES combines low cracking pres- 
sure, high strength. Well suited to both 
high pressure pneumatic systems and low 
pressure pneumatic and vacuum service, 


800 SERIES provides physical and func- 
tional characteristics comparable to swing 
check valves—insures zero leakage even 
with fuel or fuel vapor—can be mounted 


SIMPLE TO OPERATE AND MAINTAIN... 
any engineer or mathematician who 
can write the equation can solve it on 
the EASE computer with only a few 
hour’s training. 


COMPACT, COMPLETELY SELF-CONTAINED 
... the entire unit requires less than 8 
sq. ft. of floor space, is complete with 
its own regulated power supply... 
you simply plug it in to a 20 ampere, 
110 v. a.c. line! 


LOW COST ... the EASE is the world’s. 


first high- quality computer to be 
mass-produced in practical commer- 


in any position. 


200 Series 800 Series | 
Pressures: | 0 to 3000 psi | 0 to 600 psi 
Cracking | normally 
Pressure: | less than 1 psi 6° to 8 H20 
tol’ 
Sizes: Tube Tube 
air, carbon air, carbon di- 
dioxide, oxide, alcohol, 
Suitability: | nitrogen, water, aircraft 
helium, and jet fuel, 
oxygen, etc. | fuel vapor, et. 


cial form. The result is low cost with- 
out sacrifice in utility or quality. 


2200 WRIGHT AVE... 


For complete data, please request Bulletin F-9 


RICHMOND, 


BECKMAN INSTRUMENTS 


CIRCLE 4 SEAL 


precision check valves 


JAMES*DOND* CLARK 


| 2181 E. Foothill Bivd., Pasadena 8, Calif. 
Engineering data sent free of charge on request? 
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THROUGH THREE GENERATIONS... 


__ Recognized as outstanding producers of small parts and 


- close tolerance work .. . with plenty of Yankee ingenuity 
a plus a wide open mind for the best ideas of others. 


In our modern plant we also have the automatic screw 


machines for producing quality parts up to 1%” diameter 


8%” long. 


Through constant research and development work 

eG 4 in the interest of our customers we have maintained our 

leadership as specialists in fine instrument screws for over 


sixty years... with a variety of products which cover 


_ every department of the screw machine products 


MEMBER OF 
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SCREW MACHINE 
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Specialists in Screw Machine Products Since 1893 
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SUPERIOR PERFORMANCE 


= ~ On 3 Cycles to 20 Megacycles 


The 
LA-239 C 


Oscilloscope 


1. Wider Bandwidth: Complex 
waves from 5 Cycles to 15 Megacycles. 
Sine waves from 3 Cycles to 20 Mega- 
cycles. 


2. Extended Sweep Frequencies: 
Linear from 10 Cycles to 20 Megacycles 
internally synchronized. Triggered 
sweep, from single random impulses to 
irregular pulse-intervals up to as high as 
6 Megacycles. 


3. Square Wave Response: Rise 
time 0.042 Microseconds; only 5% 
droop on flat-topped pulses as long as 
30,000 Microseconds duration. 


4. Greater Stability: Electronically regu- 
lated power supplies throughout to maintain 
accuracy and constant operation under varying 
line conditions or line surges. You can display 
surges on the line from which Model LA-239C 
is being powered without distortion of the 
trace! 


5. Higher Signal Sensitivity: Maximum 
sensitivity without Probe: 10.4 millivolts. With 
Probe: 100 millivolts. (Maximum signals, 125 
V. Peak and 450 V. Peak respectively.) 


6. Timing Markers: Interval Markers of 
0.2; 1; 5; 20; 100; 500; or 2,000 Microseconds 
may be superimposed on the trace for the 
accurate measurement of the time base. 


7. Voltage Calibration: Signal amplitude 
is compared against a 1,000 cycle square 
wave (generated internally) the amplitude of 
which is controlled by a step-and-slide at- 
tenuvator calibrated in peak volts. (A jack is 
provided to deliver 40V Peak for use in 
calibrating other instruments.) 


8. Sweep Delay: Any portion of the sweep 
longer than a 10 Microsecond section may be 
expanded by 10:1 for detailed study of that 
portion of the signal. 


9. Power Source: 110 to 130 V AC; from 
50 to 1,000 cycles. 295 Watts. (Fused at 4 
Amperes.) 


10. Dimensions: In Bench Cabinet: 19 in. 
Wide; 15% in. High; 16% in. Deep. In Rack 
Mounting (With cabinet removed to fit stand- 
ard relay rack): 192 in. Wide; 14 in. High. 
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THE LAVOIE MODEL LA-239C has been designed to surpass 
the high performance of the TS-239A/UP, which has been the 
standard test oscilloscope for the Armed Services since its 
introduction. Model LA-239C is the result of a long period of 
research and development which has included the study of 
new tubes, new circuits, and new techniques. Rugged design 
has been combined with functional simplicity to produce an 
instrument as attractive as it is efficient. 


To create a circuit that will produce a certain complex 
wave form, or study transients and pulse phenomena, no 
better precision instrument is available today. 


Lavoie Laboratories take pride in offering this precision 
oscilloscope as the combination of engineering perfection 
and manufacturing skill. 


MORGANVILLE, NEW JERSEY 
DESIGNERS AND MANUFACTURERS OF ELECTRONIC EQUIPMENT 
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